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SUBSONIC AERODYNAMIC CHARACTERISTICS OF INTERACTING LIFTING 
SURFACES WITH SEPARATED FLOW AROUND SHARP EDGES 
PREDICTED BY A VORTEX -LATTICE METHOD 

John E. Lamar and Blair B. Gloss 
Langley Research Center 

SUMMARY 

Because the potential flow suction along the leading and side edges of a planform 
can be used to determine both leading- and side -edge vortex lift, the present investiga- 
tion was undertaken to apply the vortex-lattice method to computing side-edge suction 
force for both Isolated and interacting planforms. Although there is a small effect of 
boxmd vortex sweep on the computation of the side -edge suction force, the results 
obtained for a number of different isolated planforms produced acceptable agreement 
with results obtained from an application of the suction analogy to a method employing 
continuous induced-velocity distributions. The computed side -edge suction results gen- 
erally remained within 1 percent of the converged ones for 6 singularities chordwise at 
each of 20 spanwise locations. 

The lift characteristics obtained from the present method for several isolated 
wii^s agreed as closely with experimental data as did those obtained from the method 
employing continuous induced -velocity distributions. In addition, by using the method 
outlined in this report, better agreement between theory and experiment was noted for a 
wing in the presence of a canard than had previously been obtained. 

INTRODUCTION 

The development of methods for predicting the aerodynamic characteristics of air- 
craft, which during portions of their flight envelopes, develop separated flow with reat- 
tachment around the leading and side edges of the wing, has been the subject of studies 
for many years and has had varying degrees of success. Polhamus in references 1, 2, 
and 3 provided a method by which the effects of separated flow around a sharp leading 
edge on lift and drag could be estimated by an analogy that relates these forces to the 
attached flow leading -edge suction force. Hence, current lifting -surface computer pro- 
grams which estimate leading -ec^e suction are useful in predicting the leading -edge sep- 
aration effects on the forces. 


The effects of side-edge separation on the aerodynamic characteristics have been 
estimated by determining the attached flow side force and then employing a "suction 
analogy" similar to that used at the leading edge. This concept was developed and cor- 
related with experiment and other methods in reference 4. The procedure employed in 
reference 4 is based on a modified Multhopp method which is outlined in reference 5. 

The technique of reference 4 is well suited to single, isolated planforms. However, 
for lifting planforms in combinations or for flat body -wing configurations, the vortex - 
lattice method (i.e., ref. 6) is better suited because its elemental panel representation 
provides a straightforward extension to the more complex configurations. 

This paper is concerned with the extension of the vortex -lattice method of refer- 
ence 6 to the computation of the attached flow side force. Others have published similar 
work recently (refs. 7, 8, and 9, for example) but they do not provide comparisons of 
alternate paneling arrangements or convergence studies with their analyses. Thus, the 
purposes of this paper are (1) to provide comparisons of alternate paneling arrangements, 
(2) to provide convergence studies, (3) to provide the details of side -force and pitching- 
moment computation, and (4) to present comparisons between experimental and theoreti- 
cal results for sii^le planform configurations and interacting planform configurations. 

The program changes made in Langley computer program A2794 since the publica- 
tion of reference 6 are discussed in appendix A. These changes include corrections, 
improvements, and the additions relating the the side -edge -suction computations. In 
addition, input and output data for a sample case are presented in appendix B, and a list- 
ing of the program is provided in appendix C. 

SYMBOLS 


A aspect ratio 

b wing span 

Cq drag coefficient, 

^oo®ref 

^D,o experimental value of drag coefficient at Cl = 0 


Cl lift coefficient, 

^oo^ref 


pitching -moment coefficient about Y axis, Patching moment 

9oo^ref^ref 


2 



ACm, i contribution to pitching -moment coefficient from vortex system operating on 
an elemental panel 


Cn 


normal-force coefficient, 


Normal force 
QooSref 


Cs 


total leading -edge suction-force coefficient, 
2(Leadi ng-edge suction force on one panel) 
Qoo^ref 


Cy, 


twice total side -edge -suction -force coefficient of one wing panel, 

S6 N/2 

2(Side force along side edge on one wing panel) = 2 ^ se)j 

^oo®r6f 


^ACy se) i contribution to side -edge -force coefficient from ith elemental panel 


Cref 


reference chord 


ct 


tip chord 


AFy,i 

Kp 


contribution to side force from ith elemental panel 
a(sin a cos a) 


Kv,Ze 


a sin^ ot 


I 

M 

M 


Q^y,se 
a sin^ or 

length of trailing filament between adjacent chordwise horseshoe vortices 
bound vortex filament length in chord direction 
Mach number 



N 


Nc 

Ns 

Qoo 

Sref 

U 

w 

X,Y,Z 

X,Y,Z 


x,y 


total number of horseshoe vortices that contribute to the side -edge suction 
force 

number of elemental panels in chordwise row 
number of chordwise rows on wing semispan 
free -stream dynamic pressure 
reference area 

free -stream velocity 
induced downwash velocity 

axis system of a given horseshoe vortex (see fig, 1) 
body -axis system for planform input (see fig. 1) 

distance along X- and Y-axis, respectively 


Xj-ef moment reference point (taken to be zero herein) 

‘distance along tip chord to centroid of side-edge force 

Ax’ chordwise distance from midpoint of particular vortex filament to moment 

reference point 

« angle of attack, deg 

/3 = \jl - m 2 


vortex strength 

chordwise sum of vortex strengths to a particular elemental panel, sr (see 

fig. 1) 
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leading -edge sweep angle, positive for sweepback, deg 



taper ratio 


\ 

p density 

ij/ sweep angle of bound vortex, deg 

Subscripts: 

B bound vortex 

c centroid 

i particular horseshoe vortex 

j particular item of location 

L left 

le leading edge 

p potential or attached flow 

R right 

se side edge 

tot total 

vie vortex effect at leading edge 

vse vortex effect at side edge 

THEORETICAL DEVELOPMENT 

The attached flow side force is developed in accordance with the Kutta-Joukowski 
law for forces generated by a vortex filament. Figure 1 shows vortex filaments which 
have a streamwise component interacting with the net downwash at the filament midpoint 
to produce an elemental side force. The net side force on an elemental panel due to a 
swept horseshoe vortex system on the left wing panel is 
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AFy^i = p|^r- |^(wl - Va)lj^ (wr - Ua)ZRj + ^ (wg - Ua)Azj 


( 1 ) 


and the contribution to the side-force coefficient is 


(^CY,se)j Sref|u 


U 


a I 


L - 


w 


U 


a\l 


R 


tan 1^1 rf^B ' 
tan i// U\ U ^ 



( 2 ) 


The side force is of order or 2^ which is appropriate since it is associated with 
edge suction. If the trigonometric terms were retained, the side force is actually a 
function of sin^ a since the a term is really a sin a term and T/U, T'/U, and 
w/U are all proportional to sin a for the wind and body axes coincident. Hence, 
Rv,se can be formulates as 


^v,se 



(3) 


For small « the sin « - a which, for numerical purposes, is taken to be 1 radian in 

equations (2) and (3) and leads to 


N/2 


*^v,se “ 2 ^ pCy^gg^_ 

i=l 


(4) 


This provides an additional contribution to the total lift, as indicated in the following 
equation: 


^L,p 

^L,tot = Kp sin Q! cos" 



Cl 


,vse 


sin O' sin Q! cos a 




(5) 
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For planforms having sharp edges, the drag coefficient can be written as 
Cd = Cd,o + Cl, tot tan a 



In the numerical determination of the side force, it is realized that computational 
time savings could be made with the utilization of a swept horseshoe vortex system. The 
savings are due to the vortex filament length and the net downwash associated with the 
right trailing filament of a swept horseshoe vortex being the same as those for the left 
trailing filament on the adjoining inboard swept horseshoe vortex. However, the swept 
bound vortex may lead to local and overall errors in the side force, just as it did for 
leading-edge thrusfl (ref. 10). This potential problem will be investigated although it 
should be less serious than that for thrust because in the side -force computation, the 
bound vortex interaction with the net downwash at its midpoint will only contribute a por- 
tion to the total side force, rather than the entire result. 

To assess the importance of the inclusion of the swept bound vortex, numerical 
studies are presented in the next section. They are based on paneling the wing in vari- 
ous ways to emphasize the influence of the bound vortex differently. 

The pitching-moment contribution about the Y-axis associated with the side -edge - 
suction force is determined from each elemental horseshoe vortex by 


AC 


m.i 


_2__JL 
Sref ^ref 




Ax^ 





|tan \l/\ p/WB 
tan i// U\ U 



The sign of each term is chosen with the realization that the overall rotation of the 
trailing flow field on the left wing panel is clockwise as viewed from the rear. This 
rotation causes the vortex elements behind the moment reference point to contribute a 
noseup moment if associated with the left trailing leg or a sweptback bound vortex and a 
nosedown moment for the vortex elements associated with the right trailing leg (fig. 1). 
The total pitching moment is obtained by using the following expression. 


'm,p 




^m,vse " 



m,i 


sin or sin a 


^c,le 




Cm.tot = Kp sin » cos « ^ + K^,le Isin a | sin a -5^ + 




'‘C,se 


v,se a I sin a 


when the particular Xc -terms equal Xref ' 

Only those horseshoe vortices or portions thereof which would intersect the side 
edge if they were projected laterally to the local spanwise extent of the planform (those 


iThe leading -edge thrust problem and the program changes made to correct it are 
described in appendix A. 
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that do so are said to directly oppose the side edge) are included in the summation for 
the side force and pitching moment. This procedure is the same as that used for com- 
puting the leading -edge thrust, where all the distributed thrust along the chord is pro- 
jected forward and assumed to act at the leading edge. One reason for computing the 
side -edge suction force in this manner, rather than with the method presented in refer- 
ence 8, is that, in the application of the method of reference 8 to a cropped diamond wing, 
the entire side -edge suction force would be calculated over a wing panel with the only 
reduction coming from the removal of the contribution from the leading -edge suction. 
This retains the contribution to the side force from the aft part of the cropped delta wing 
from which no edge force is expected. 

NUMERICAL STUDIES 
Panel Arrangements 

Table I presents a comparison of Ky ge and the side -edge load centroid obtained 
by four different methods for the three wing planforms presented in figure 2. Method 1 
is the continuous loading method of reference 4 and the results of this method are taken 
to be the standard. Method 2 is the present method, which was described previously. 
Method 3 is the same as method 2 except each planform is considered as two wings (the 
dashed lines in fig. 2 show break lines) and the side-edge suction force is computed only 
on the aft wing. Method 4 is the vortex -lattice method described in reference 7 with the 
results being supplied by R. G. Bradley of General Dynamics Corp. 

The layout of the bound vortices directly opposing the side edge of a wing tip will 
have less sweep for method 3 than for method 2 because method 3 panels the wing as two 
planforms. In fact the cropped delta wing is a special case for method 3 since the bound 
vortices inboard of the wing tip have no sweep. By comparison, the results of method 3 
agreed more closely with the continuous induced-velocity approach of reference 4 
(method 1) than the results of method 2. The results from method 3 for the cropped delta 
wing agreed closest with method 1. This agreement indicates that there is an effect of 
bound vortex sweep on Ky^ge; however, the maximum difference between the results of 
method 2 and those of method 1 for the configurations shown in table I was only about 
4 percent for Ky^ge and about 4 percent for the centroid location. The table also shows 
that the results of method 4 (ref. 7) are somewhat higher than those of method 1. 

In an effort to study further the effect of bound vortex sweep angle on Ky ge and 
its chordwise centroid, results were obtained for a rectangular wing sheared to various 
sweep angles; these results are presented in table II. The reason for selecting this type 
of planform was to provide a critical evaluation of method 2, since the bound vortex sweep 
angles will all be (1) the same, (2) maximum for the planform, and (3) equal to that of the 
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leading ec^e. For any other simple planform, discounting those with reversed taper, the 
bound vortex sweep angles would become less positive the closer the vortices are located 
to the trailing edge. Increasing the leading -edge sweep angle in this manner leads to a 
reduced number of horseshoe vortices that directly oppose the side edge and that can 
make a contribution to side force or its moment. Table II shows that, in general, the 
Kv gg values and chordwise centroid locations of method 2 are smaller and more for- 
ward, respectively, than those of method 1. The maximum percent errors are 9.4 for 
Ky se Ax/ct in terms of ct. The highest sweep angle reported in this 

study was 75° because for higher sweep angles method 1 was unable to determine a suit- 
able control point pattern to insure a valid logarithmic singular correction. 

Method 2 is chosen in this paper as the method to be used in subsequent calculations 
and will be designated as present method because it allows the analysis of two wings in 
the presence of each other (for example, a canard -wing configuration), whereas method 3 
does not. Another reason is that the values of Ky^ge presented in table I, as deter- 
mined by the four methods, have only small differences and the effect on the total lift 
answer would amount to less than a 4 -percent maximum error for aisles of attack up to 
30°. 

Table HI presents a comparison of Kp, Ky^^e, and Ky gg as computed by the 
present method and method 1 (ref. 4) for several different aspect -ratio rectangular wings 
and the three wing planforms shown in figure 2. There is seen to be good agreement 
between the two methods. As the aspect ratio for the rectangular wings approaches zero, 
Ky gg is less than n for the present method and greater than tt for method 1. Ref- 
erence 4 shows for rectangular wings that as |3A approaches zero, the theoretical value 
of Ky gg should approach tt. 

Values of Ng and Ng ^ and 25, respectively, were used in obtaining the results 
in table HI because a preliminary investigation indicated this combination to be adequate. 
Subsequently, a convergence study was undertaken to determine the minimum require- 
ments of Ng and Ng- The results of this study are discussed in a subsequent section. 

Effect of Ng and Ng 

The effect of varying Ng and Ng on Ky^gg and Ax/ct is examined herein. 

The solutions are from the present method and are presented in figures 3 and 4 for iso- 
lated planforms. For most of the wings considered in figures 3 and 4, there are many 
combinations of Ng and Ng which will yield results within 1 percent of what a^ears 
to be the converged values of Ky gg and Ax/c^. In particular the pattern for Ng - 6 
and Ng = 20 gives good agreement with the converged result except for the wing with 
A = 3.5 and A = 75°. The vortex -lattice representation for this wing does not provide 
enough vortices that oppose the side edge for a converged result to be determined with 
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method 2. This pattern is seen to be larger than the pattern for Ng = 4 and Ns = 20, 
which was determined in reference 6 to be generally adequate to yield acceptable 
9Cm/9CL solutions. 


CORRELATION WITH EXPERIMENT 

Figures 5 to 8 present the comparison of the theoretical results obtained by the 
present method with experimental data obtained from reference 4 for flat rectangular 
wings with sharp leading and side edges. These wings had aspect ratios of 0.20, 0.40, 
1.00, and 3,00. Also, the comparison of theoretical results (present method) with exper- 
imental results (ref, 4) for the three swept wings presented in figure 2 is shown in fig- 
ures 9 to 11. Since the values of and its chordwise centroid obtained by 

methods 1 and 2 and reported in table HI closely agree, the theoretical results obtained 
by using the present method should agree as well with the experimental data as did the 
results of method 1 (ref. 4). As pointed out in reference 4, the reason for the disagree- 
ment between experiment and theory shown in figures 9 and 11 is that sweptback wings 
which have large amounts of area behind the point of maximum span develop lift values 
in excess of those predicted because of the additional induced effects associated with the 
actual shed -vortex system. These additional induced effects are the leading -edge vortex 
acting along the side edge and over the trailing triangular portion of a cropped diamond 
wing. 

Figures 12 to 14 present a comparison of experimental lift (ref. 9) with the present 
theoretical lift on a wing in the presence of a canard for three different canard configu- 
rations: (1) in the wing chord plane, (2) above the wing chord plane, and (3) above the 
wing chord plane with 18. 6o of anhedral. (Only the lift on the portion of the model drawn 
With solid lines in the sketches in figs. 12 to 14 is plotted.) The addition of the side-edge 
vortex lift gives better agreement between theory and experiment as compared with that 
shown in reference 11 for these particular models at angles of attack up to wing stall. 

The theoretical results presented in reference 11 are replotted in figures 12 to 14 as 
short dashed lines. It should be noted that the canard data are not presented in this report 
since the experimental data of reference 11 indicated the canard did not develop full 
leading-edge vortex lift. 


CONCLUDING REMARKS 

Because the potential flow suction along the leading and side edges of a planform 
can be used to determine both leading- and side -edge vortex lift, the present investiga- 
tion was undertaken to apply the vortex -lattice method to computing side -edge suction 
force for isolated or interacting planforms. Although there is a small effect of bound 
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vortex sweep on the computation of the side -edge suction force, the results obtained for a 
number of different isolated planforms produced acceptable agreement with results 
obtained from a method employing continuous induced -velocity distributions. The com- 
puted side -edge suction results generally remained within 1 percent of the converged ones 
for 6 singularities chordwise at each of 20 spanwise locations. 

The lift characteristics obtained from the present method for several isolated 
wings agreed as closely with experimental data as did those obtained from the method 
employing continuous induced -velocity distributions. In addition, by using the method 
outlined in this report, better agreement between theory and experiment was noted for a 
wing in the presence of a canard than had previously been obtained. 

Lar^ley Research Center 

National Aeronautics and Space Administration 
Hampton, Va., 2366 5 
June 19, 1975 
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APPENDIX A 


CHANGES, IMPROVEMENTS, AND ADDITIONS TO 
LANGLEY COMPUTER PROGRAM A2794 

The purpose of this appendix is to describe in some detail the changes, improve- 
ments, and additions to Langley computer prc^ram A2794 of reference 6, 

Changes to the input cards of the computer program presented in reference 6 are 
as follows: 

(1) On the configuration card an additional field of 5 has been added (columns 66 
to 70) and has a specification of F5.1. This field is used to obtain entry into the tip- 
suction overlay. By putting a 1. in this field, a tip -suction computation will be made. 

If this option is not required, the field can be left blank or a 0. put in it. 

(2) With the 1. specified in change (1), it is necessary to provide another input data 
card with a format of 4F10.5. This card contains the limits of "y-integration" over 
which the leading -edge suction distribution is to be integrated. Normally these limits 
would be the plane of symmetry (0.) and the left wing tip (-b/2); however, others could be 
used. Four fields are provided for the two planforms since each planform would need a 
beginning (inboard) and ending (outboard) y-location. The order is for the first planform 
beginning and ending y-limits followed by the second planform y-limits. This card then 
b6conies the l3.st card of the input deck for a configuration, 

(3) With the 1. specified in change (1) and SCW specified as 0. the numbers in 
TBLSCW(I) must be the same on a given planform but can change with planform in order 
for the program to function properly. This restriction was placed on TBLSCW(I) to save 
computer execution time. 

(4) With the 1. specified in change (1), SCW or TBLSCW(I) must be larger than 1. 

Program Improvements 

Improvements to the vortex -lattice computer program presented in reference 6 are 
detailed below. A listing of the complete revised computer program is presented in 
appendix B. 

Since the leading-edge thrust and its distribution are obtained by the difference 
between Lift x a and induced drag on an overall and local basis, improvements in the 
accuracy of these terms would yield necessarily more accurate thrust results. Refer- 
ence 6 has determined ranges of Ng and Ng required for the convergence of lift for a 
wide assortment of planforms. A similar study for the near-field induced drag found 
convergence but not to the far -field values. Reference 10 relates this problem to the 
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violation of the Munk stagger theorem (ref. 12) in two regards: The induced drag asso- 
ciated with the induced velocity from the bound vortices and from the near -field limit of 
the trailing vortices should sum to zero, and the bound vortices should be of uniform 
length. 

To implement those ideas, a new overlay was developed (OVERLAY 4) to take the 
set of r/U already determined in OVERLAY 2 along with the original vortex lattice in 
OVERLAY 1 and to repanel the wing with a network of unswept horseshoe vortices of 
equal spacing whose values of P/U were determined by Interpolation. It is with this 
setup that the near-field induced drag is now determined. 


OVERLAY 5 

The side -force computation is performed in OVERLAY 5 as outlined in the section 
of this paper entitled "Theoretical Development." Coplanar wings which have an unswept 
leading edge on the second wing require special attention in the computation of side force 
and pitching moment at the leading -edge region of the second wing. The side force act- 
ing on the trailing vortex filaments of the first wing which intersect the leading edge of 
the second wing is computed between lines A-A and B-B (shown in sketch (a)) of the sec- 
ond wing. 

Wing planform 



Sketch (a). 
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OVERLAY 5 computes integrating the local leading -edge suction over 

the desired portion of the configuration to obtain the total leading -edge suction coeffi- 
cient Cg. This was done to allow the program user the flexibility of choosing the 
leading-edge region of the wing over which vortex lift is assumed to exist. By making 
the small angle of attack approximation, which is done throughout the potential flow part 
of this computer program, is computed by using the following expression: 


The values of Kp and computed by this program in the manner outlined are 

appropriate only for untwisted or uncambered liftii^ surfaces. The lift and pitching- 
moment coefficients are then computed by using the expressions found in reference 4. 
The side -force computation on planforms with dihedral is performed in a manner sim- 
ilar to that for planforms with no dihedral. For reliable side -edge loading results the 
program should be restricted to planforms which do not have swept forward leading 
edges. 
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SAMPLE CASE 

Input data, the sketch, and output data for a sample case with a canard -body -wing 
combination are presented in this appendix. The canard shown has 18.6° anhedral, and 
the leading -edge suction is integrated from the body -wing and body-canard intersections 
to the respective tips. 

The following list contains the output variable names not defined in reference 6: 


KP 



KV LE 

Ky,Ze 


KV SE 

I^v,se 


ALPHA 

a 


CN 

^N,tot 


CLP 

Cl,p 


CLVLE 

^L,vZe 


CLVSE 

Ky^se |sin Q!l sin a cos a 


CMp2 

pitching -moment coefficient due to 

^L,p 

CMVLe2 

pitching -moment coefficient due to 

^L,vZe 

CMVSe2 

pitching -moment coefficient due to 

^L,vse 

cm2 

total pitching moment 


CD 

f^L,tot “ 


CL**2/(PI*AR) (Cl, tot) 



^Reference point is the origin of the X,Y,Z axis system. 
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Sample Input Data and Sketch 


COLUMN NUMBERS FOR INPUT DATA 

000000000 1 1 1 1 1 1 1 1 1 12^222?2a2^33333333314A44.A<4^^^,aS55555S5556b666f>666677777777778 
1 2345670901 234S6709O1 234567800 12345678R0 1234567890 1234567890 1234567890 1234567890 


GROUP ONE DATA 


2. 

1. 

9.1756 

10. 

0. 

0. 

0.0 

0,0 

0.0 

-1.2 

-.55 

0. 

-2.2 

-.9 

0. 

-3.2 

-1.2 

0. 

-4.2 

-1.42 

0. 

-5. 

-1.5 

0. 

-11.65 

-1.5 

fVJ 

. 

GO 

1 

-18.60 

-6.51 

0. 

-20.00 

-6.51 

-18.62 

-18.85 

-1.5 

0. 

-18.85 

0.0 


7. 

0. 

0. 

-18.05 

0.0 

0. 

-18.85 

-1.5 

0. 

-20,4 

-1.5 

0. 

-28.7 

-10. 

0. 

-30.9 

-10, 

0, 

-32.2 

-1.5 

0. 

-38. 

-1.5 

0. 

-38. 

0,0 



GROUP TWO 

DATA 

1. 

6. 13. .30 

1 . 0, 

-1.5 

-15. 

-1.5 


159.99696 

-1.69 

1 . 

1 . 

1 . 

1 • 

1 • 

1 . 

1 . 

1 . 

1 • 

1 . 

0.0 

1 . 
u 
1 . 

1 . 

1 • 

1 . 

1 . 


0 . 0 , 0 . 

-15. 


1 . 
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Sample Output Data 

GEOMETRY DATA 


FIRST REFERENCE PLANFORM HAS 10 CURVES 

VARIABLE S^^EEP PIVOT POSITION XtSI * 

BREAK POINTS FOR THE REFERENCE PLANFORM 

SWEEP DIHEDRAL 

ANGLE ANGLE 


ROOT CHORD HEIGHT = 


POINT 


1 

2 

3 

4 

6 

7 

8 
9 

10 

ll 


-1.69000 


X 

Y 

REF 

REF 

0.00000 

0.00000 

-1.20000 

-.55000 

-2.20000 

-.90000 

-3.20000 

-1.20000 

-4.20000 

-1.42000 

-5.00000 

-1.50000 

-11.65000 

-1.50000 

-ia.60300 

-6.51000 

-20.00000 

-6.51000 

-18.85000 

-1. 50000 

-18.85000 

0.00000 


65.37644 

0.00000 

70.70995 

0. 00000 

73.30076 

0.00000 

77.59258 

0. 00000 

84.28941 

0.00000 

90. 00000 

0. 00000 

54.21355 

-18.62000 

90.00000 

0. 00000 

12.92778 

-18.62000 

0.00000 

0.00000 


0.00000 YiSi * 


MOVE 

CODE 

1 

1 

I 

I 

1 

1 

I 

1 

I 

I 




SECOND REFERENCE PLANFORM HAS 

7 CURVES 


ROOT CHORD HEIGHT = 

O.dlDOO 

VARIABLE SWEEP 

PIVOT POSITION 

X<S) = 

0.00000 



BREAK POINTS FOR 

THE REFERENCE 

PLANFORM 


POINT 

X 

REF 

Y 

REF 

SWEEP 

ANGLE 

DIHEDRAL 

ANGLE 

MOVE 

CODE 

1 

2 

3 

4 

5 

6 

7 

8 

-18.85000 
-18.85000 
-20.41003 
-28. 70000 
-30.90000 
-32.20000 
-38. 30000 
-38.00000 

0.00000 
-1.50000 
-1.50000 
-10.00000 
-10.00030 
-1.50000 
-1. 50000 
0.00000 

0.00000 

90.00000 

44.31794 

90.00000 

-8.69550 

90.00000 

0.00000 

0.00000 
0.00000 
0. 00000 
0.00000 
0. 00000 
0.00000 
0.00000 

1 

1 

1 

1 

1 

1 

1 


Y(SJ * 


0.00000 


0«00000 
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00 


CONFIGURATION NO. 1 


CURVE 

i IS 

SHEPT 

65.37644 

DEGREES. ON 

PLANFORM 

1 

CURVE 

I IS 

SWEPT 

0.00000 

DEGREES ON 

PLANFORH 

2 




BREAK points FOR 

THIS CONFIGURATION 


POINT 

X 

Y 

1 

SWEEP 

DIHEDRAL 

HOVE 





ANGLE 

ANGIE 

CODE 

I 

0.00000 

0.00000 

“1.69000 

65.37644 

0. 00000 

1 

2 

“1.20000 

“.55000 

“1.69000 

70.70995 

0.00000 

1 

3 

“2.20000 

“.90000 

“1.69000 

73.30076 

0.00000 

1 

4 

“3.20000 

“1.20000 

“1.69000 

77.59258 

0.00000 

1 

5 

“4.20000 

“1.42000 

“ 1.69000 

84.23941 

0.00000 

1 

6 

“5.00000 

“1.50000 

“1,69000 

90.00000 

0.00000 

1 

7 

“11.65300 

“1.53300 

“ 1.69000 

54.21355 

“18.62000 

1 

8 

“la. 60000 

“6.51000 

-.00200 

90.00000 

0.00000 

1 

9 

“20.00000 

“6.51300 

-.00200 

12.92778 

-18.62000 

1 

10 

“18.85000 

“1.50000 

“1.69000 

0.00000 

0.00000 

1 

1 1 

“18.85000 

3.00000 

“ 1.69000 







SECOND PLANFORH BREAK POINTS 



1 

“18.85000 

0.00000 

0.00000 

0.00000 

0.00000 

1 

2 

“18.85000 

“.55000 

0,00000 

0. 00000 

0.00000 

1 

3 

“ 18 . 35000 

-.90000 

0.00000 

0.00000 

0.00000 

1 

4 

“18.85300 

“1.23000 

0.00000 

0.00000 

0.00000 

1 

5 

“18.85000 

-1.42000 

0.00000 

0.00000 

0.00000 

1 

6 

“18.85000 

“1.50000 

0.03000 

90.00000 

0. 00000 

1 

7 

“20,40000 

“1.50000 

0.00000 

44.31794 

0.00000 

1 

8 

-25.29212 

-6. 51000 

0. 0 3000 

44.31794 

0.00000 

i 

9 

“28. 70000 

“10.00000 

0.00000 

90.00000 

0. 00000 

1 

10 

“30.90000 

-10.00000 

0., 00000 

“8.69550 

0.00000 

1 

11 

“32.20000 

“L. 50000 

0.00000 

90.00000 

0. 00000 

1 

12 

“38.00000 

-1.50000 

0.00000 

0.00000 

0. 00000 


13 

-38.00000 

0.00000 

0.00000 





174 HORSESHOE VORTICES USED ON THE LEFT HALF OF THE CONFIGURATION 


PLANFORH 


TOTAL 


SPANWISE 
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AERODYNAMIC DATA 


CONFIGURATION NO. 1 


STATIC LONGITUDINAL AERODYNAMIC COEFFICIENTS ARE COMPUTED 


X X 

C/4 3C/4 


-18. 17029 

-18.32212 

-6.1^552 

-18.47395 

-18.62578 

-6.14552 

-18.77761 

-18*92944 

-6.14552 

-19.08127 

-19.23310 

-6.14552 

-19.38493 

-19.53676 

-6. 14552 

-19.68859 

-19.84042 

—6. 14552 

-17.19422 

-17.41637 

-5.41655 

-17.63853( 

-17.86068 

-5.41655 

-18.08284 

-18.30500 

-5.41655 

-18-52715 

-18.74931 

-5. 41655 

-18.97146 

-19.19362 

-5.41655 

-19.41577 

-19.63793 

-5.41655 

-16.21814 

-16.51062 

-4.68758 

-16.80310 

-17.09558 

-4.68758 

-17.33807 

-17.68055 

-4.63758 

-17.97303 

-13.26551 

-4.68758 

-13.55799 

-18.85048 

-4.68758 

-19.14296 

-19.43544 

-4.68758 

-15.24206 

-15.60487 

-3.95862 

-15.96768 

-16.33048 

-3.95862 

-16.69329 

-17.05610 

-3.95862 

-17.41891 

-17.78172 

-3.95862 

-18.14452 

-18.50733 

-3.95862 

-18.87014 

-19.23295 

-3.95862 

-14.26598 

-14.69911 

-3.22965 

-15.13225 

-15.56538 

-3.22965 

-15.99852 

-16.43165 

-3.22965 

-16.86479 

-17.29792 

-3.22965 

-17.73105 

-18.16419 

-3.22965 

-18.59732 

-19.03046 

-3.22965 

-13.28990 

-13.79336 

-2.50068 

-14.29682 

-14.80028 

-2.50068 

-15.30374 

-15.80720 

-2.50068 


I 


-.1Z481 
-.12481 
-.12481 
-.12481 
-.12481 
12481 
-.37041 
-.37041 
-.37041 
-.37041 
-.37041 
-.37041 
-.61602 
-.61602 
-.61602 
-.61602 
-.61602 
-.61602 
-.86163 
-.86163 
-.86163 
-.86163 
-.86163 
-.86163 
-1.10724 
-1.10724 
-1.10724 
-1.10724 
-1. 10724 
-1.10724 
-1.35284 
-1.35284 
-1.35284 


S C/4 SWEEP 



ANGLE 

.38462 

53.24644 

.38462 

48.39306 

.38462 

43.62423 

.38462 

37.24015 

.38462 

29-56186 

.38462 

20.51844 

.38462 

53.24644 

.38462 

48.89306 

.38462 

43.62423 

.38462 

37.24015 

.38462 

29.56136 

.38462 

20.51344 

.38462 

53.24644 

.38462 

48.39306 

.38462 

43.62423 

.38462 

37.24015 

.38462 

29.56186 

.38462 

20.51844 

.38462 

53.24644 

.38462 

48.39306 

.38462 

43.62423 

.38462 

37,24015 

.38462 

29. 561 86 

.38462 

20.51344 

.38462 

53.24644 

.38462 

43.39306 

.38462 

43.62423 

.38462 

37,24015 

.38462 

29.56136 

.38462 

20.51344 

.38462 

53.24644 

.38462 

43.39306 

.38462 

43.62423 


DIHEDRAL 

LOCAL ALPHA 

ANGLE 

IN RADIANS 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-13.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18-62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18 .62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 

-18.62000 

0.00000 


DELTA CP AT DESIRED 
CL * 1. 00000 

3.72658 

1.69346 

1.04361 

.66053 

.41117 

.22533 

3-41634 

1.60920 

1.06980 

.74750 

.50397 

.29033 

3.06694 

1.45406 

,99044 

.71586 

.50163 

.29981 

2.75671 

1.31328 

.90859 

.66947 

.47964 

.29326 

2.48005 

1.18951 

,83349 

.62596 

.45656 

.28069 

2.22473 

1.08091 

.76725 


O 
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CO 

o 


“ 16.31066 

- 16.81413 

- 17.31759 

- 17.82105 

- 18.32451 

- 18.82797 

- 12.37593 

- 12.94524 

- 13.51455 

- 14.08387 

- 14.65318 

- 15.22249 

- 15.79180 

- 16.36111 

- 16.93043 

- 17.49974 

- 18.06905 

- 18.63836 

- 5.19375 

- 6.38125 

- 7.56875 

- 8. 75625 

- 9.94375 

- 11.13125 

- 12.31875 

- 13.50625 

- 14.69375 

- 15.88125 

- 17.06875 

- 18.2562 5 

- 4.33125 

- 5.59375 

- 6.85625 

- 8.11875 

- 9.38125 

- 10.64375 

- 11.90625 

- 13.16875 

' - 14.43125 

- 15.69375 

- 16.95625 

- 18.21875 

- 3.37292 

- 4.71875 

- 6.06458 

- 7.41042 

- 8.75625 

- 10.10208 

- 11.44792 

- 12.79375 

- 14. 13950 

- 15.48542 

- 16.83125 

- 18.17708 

- 2.41458 

- 3 . 843 75 

; - 5.27292 

- 6.70208 

' - 8.13125 

- 9.56042 

^- 10.98958 

- 12.41875 

- 13.84792 

- 15.27708 

- 16.70625 

- 18.13542 

- 1.36042 

- 2.88125 

- 4.40208 

- 5.92292 

- 7.44375 

- 8.96458 

- 10.48542 

- 12.00625 

- 13.52708 

- 15.04792 

- 16.56875 

- 18.08958 

SECOND PLANFORM HORSESHOE 

- 28.43420 

- 28.65373 

- 28.87327 

- 29.09280 

- 29.31233 

- 29.53186 

- 29.75140 

- 29.97093 

- 30.19046 

- 30.40999 

- 30.62952 

- 30.84906 

-2 7.71927 

- 28.01120 

- 28.30313 

- 28. 59506 

- 28.88699 

- 29.17892 


- 2.50068 

- 1.35284 

- 2.50068 

- 1.35284 

- 2.50068 

- 1.35284 

- 1.81810 

- 1.58282 

- 1.81810 

- 1.58282 

- 1.81810 

- 1.58282 

- 1.81810 

- 1.58282 

- 1.81810 

- 1.58282 

- 1.81810 

- 1. 58282 

- 1.46000 

- 1.69000 

- 1.46000 

- 1.69100 

- 1.46000 

-1 .69000 

- 1.46000 

- 1.69000 

- 1.46000 

- 1.69000 

- 1.46000 

- 1.69000 

- 1.31000 

- 1.69000 

- 1.31000 

- 1 .69000 

- 1.31000 

- 1.69000 

- 1.31000 

- 1.69000 

- 1.31000 

- 1.69000 

- 1 .31000 

- 1.69000 

- 1.05000 

- 1.69000 

- 1.05000 

- 1.69000 

- 1.05000 

- 1.69000 

- 1.05000 

- 1.69000 

- 1.05000 

- 1.69000 

- 1.05000 

- 1.69000 

-. 72500 

- 1.69000 

-.72500 

- 1.69000 

-.72500 

- 1.69000 

-.72500 

- 1.69000 

-.72500 

- 1.69000 

-.72500 

- 1 . 69 J 00 

-.27500 

- 1.69000 

-.27500 

- 1.69000 

-.27500 

- 1.69000 

-.27500 

- 1.69000 

-.27500 

- 1.69000 

-.27500 

- 1.69000 

VORTEX DESCRIPTIONS 

- 9.61538 

0.00000 

- 9.61538 

0 . 00000 

- 9.61538 

0.00000 

- 9.61538 

0.01300 

- 9.61538 

0.00000 

- 9.61538 

0.00000 

- 8.84615 

0.00000 

- 8.84615 

0.00000 

- 3.84615 

0.00000 


.38462 

37.24015 

.38462 

29.56186 

.38462 

20.51844 

.33567 

53.24644 

.33567 

48.89306 

.33567 

43.62423 

.33567 

37.24015 

.33567 

29.56186 

.33567 

20.51844 

.04000 

84.04287 

.04000 

82.80077 

•04000 

80.909 72 

.04000 

77.69198 

.04000 

71.07536 

.04000 

51.34019 

.11000 

77.0 7090 

.11000 

74-46967 

.11000 

70 . 6 0793 

.11000 

64.35899 

.11000 

52.97327 

.11000 

29.60445 

.15000 

72.61761 

.15000 

69.24593 

.15000 

64.35899 

.15000 

56.79343 

.15000 

44.19307 

.15000 

22.61987 

.17500 

69.93693 

.17500 

66. 14953 

.17500 

60.75117 

.17500 

52.63333 

.17500 

39.80557 

.17500 

19.65382 

.27500 

64.44004 

.27500 

59.93142 

.27500 

53.74616 

.27500 

45.00000 

.27500 

32.47119 

.27500 

15.25512 


«38462 

42.90475 

.38462 

36.54497 

• 38462 

28.94001 

.38462 

20.03721 

.38462 

10.00798 

.38462 

-.67404 

.38462 

42.90475 

.38462 

36.54497 

.38462 

28.94001 


- 18.62000 

0.00000 

- 18.62000 

0.00000 

- 18.62000 

0.00000 

- 18.62000 

0.00000 

- 18.62000 

0.00000 

- 18.62000 

0.00000 

- 18.62000 

0.00000 

- 18.62000 

0.00000 

-18 .62000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0.00000 

o.ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


• 5«092 
. 4^657 
• 26 X 38 
2-03972 
.87773 
•80396 
•^6991 
.49827 
.22120 
.26730 
. 039 X 6 
.05710 
1.05621 
.69435 
-40845 
.40797 
.06902 
.08068 
.72423 
.69553 
.41693 
.48129 
.10972 
.09115 
.52299 
.65413 
.41141 
.50002 
.14697 
.09444 
.40121 
.61283 
.40307 
.48819 
.17141 
.09460 
.32584 
.56958 
.39267 


4.03015 

1.54306 

.79538 

.45779 

.27680 

.15319 

4.06426 

1.72615 

1.01160 
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- 29.47085 

- 29.76278 

- 30.05471 

- 30.34664 

- 30.63857 

- 30.93051 

- 27.00434 

- 27.36867 

- 2/. 73299 

- 28.09732 

- 28.46165 

- 28.82596 

- 29. 19031 

- 29.55464 

- 29.91897 

- 30.28330 

- 30.64762 

- 31.01195 

- 26.28940 

- 26.72613 

- 27.16286 

- 27.59959 

- 28.03631 

- 28.47304 

- 28.90977 

- 29.34649 

- 29. 78322 

- 30.21995 

- 30.65667 

- 31.09340 

- 25.73998 

- 26.23234 

- 26.72471 

- 27.21707 

- 27. 70944 

- 28.20180 

- 23.69417 

- 29.18653 

- 29.67890 

- 30.17126 

- 30.66363 

- 31.15599 

- 25.19055 

- 25.73856 

- 26.28656 

- 26.83456 

- 27.38257 

- 27.93057 

- 28.47857 

- 29.02657 

- 29.57458 

- 30.12258 

- 30.67058 

- 31.21859 

- 24.47562 

- 25.09602 

- 25.71642 

- 26.33683 

- 26.95723 

- 27.57763 

- 28.19803 

- 28.81843 

- 29.43883 

- 30.05923 

- 30.67963 

- 31.30003 

- 23.76069 

- 24.45349 

- 25. 14629 

- 25.83909 

- 26.53189 

- 27.22469 

- 2/. 91749 

- 28.61029 

- 29.30308 

- 29.99588 

- 30.68868 

- 31.38148 

- 23.04576 

- 23.81095 

- 24.57615 

- 25.34135 

- 26.10655 

- 26.87175 

- 27.63694 

- 28.40214 

- 29.16734 

- 29.93254 

- 30.69773 

- 31.46293 

- 22.33083 

- 23.16842 

- 24.00602 

- 24.84361 

- 25,68121 

- 26.51880 

- 27.35640 

- 28.19400 

- 29.03159 

- 29.86919 

- 30.70678 

- 31.54438 


- 8.84615 

0.00000 

- 8.84615 

o.onoo 

- 8.84615 

0.00000 

- 8 . 0/692 

0.00000 

- 8.07692 

0.00000 

- 8.07692 

0.00000 

- 8.07692 

0.00000 

- 8.07692 

0.00000 

- 8.07692 

0.00000 

- 7.30769 

0.00000 

- 7.30769 

0.00000 

- 7.30769 

0.00000 

- 7.30 769 

0.00000 

- 7.30769 

0.00000 

- 7.30769 

0.00000 

- 6.71654 

0.00000 

- 6.71654 

0.00000 

- 6.71654 

0.00000 

- 6.71654 

0.00000 

- 6.71654 

0.00000 

- 6.71654 

0.00000 

-^.12538 

o.ouooo 

- 6.12538 

0.00000 

- 6.12538 

0.00000 

- 6.12538 

0.00000 

- 6.12538 

0.00000 

- 6.12538 

0.00000 

- 5.35615 

0.00000 

- 5.35615 

0.00000 

- 5.35615 

0.00000 

- 5.35615 

0.00000 

- 5.35615 

0.00000 

- 5.35615 

0.00000 

- 4.58692 

0.00000 

- 4.58692 

0.00000 

- 4.58692 

0.00000 

- 4.58692 

3.00000 

- 4.58692 

0.00000 

- 4.58692 

0.00300 

- 3.81769 

0.00000 

- 3.81769 

0.00000 

- 3.81769 

0.00000 

- 3.81769 

0.00000 

- 3.81769 

0.00000 

- 3.81769 

0.00000 

- 3.04846 

3.00000 

- 3.04846 

0.00000 

- 3.04846 

0.00000 

- 3.04846 

0.00000 

- 3.04846 

0.00000 

- 3.04846 

0.00000 


.38462 

20.03721 

.38462 

10.00798 

.38462 

-.67404 

•38462 

42.90475 

.38462 

36.54497 

.38462 

28.94001 

.38462 

20.03721 

.38462 

10.00798 

.38462 

-.67404 

.38462 

42.90475 

.38462 

36. 54497 

•38462 

28.94001 

.38462 

20.03721 

.38462 

10.00798 

•38462 

-.67404 

.20654 

42.90475 

.20654 

36.54497 

.20654 

28.94001 

.20654 

20.03721 

.20654 

10.00798 

.20654 

-.67404 

.38462 

42.90475 

.38462 

36.54497 

.38462 

28.94001 

.38462 

20.03721 

.38462 

10.00798 

•38462 

-.67404 

.38462 

42.90475 

*38462 

36.54497 

.38462 

28.94001 

.38462 

20.03721 

.38462 

10.00798 

•38462 

-.67404 

.38462 

42.90475 

.38462 

36.54497 

• 38462 

28.94001 

.38462 

20.03721 

.38462 

10.00798 

.38462 

-.67404 

.38462 

42.90475 

.38462 

36.54497 

•38462 

28.94001 

.38462 

20.03721 

•38462 

10.00798 

.38462 

-.67404 

•38462 

42.90475 

.38462 

36.54497 

.38462 

28.94001 

.38462 

20.03721 

•38462 

10.00798 

•38462 

-.67404 


o.ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

O.ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0.00000 


IS3 


.62038 

.37816 

.20487 

3.90372 

1.67207 

1.01718 

.64912 

.40451 

.21969 

3.75602 

1.55480 

.95388 

.62259 

.39546 

.21653 

3.75572 

1.38093 

•87260 

•58406 

•37780 

.20889 

1.55688 

1.12328 

.78078 

.54240 

.35835 

•20024 

1.20203 

•80400 

•64936 

•48163 

.32869 

•18653 

1.07099 

.64833 

.54302 

.42692 

.30159 

.17359 

•97541 

•56518 

.47066 

.38345 

.27953 

.16170 

.88966 

.51438 

.42166 

.35200 

.26443 

.14966 
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to 

to 


- 21.61539 

- 22.52589 

- 2.27923 

0.00000 

- 23.43538 

- 24.34537 

- 2.27923 

0.00000 

- 2>. 25537 

- 26.16586 

- 2.27923 

0.00000 

- 2 ^. 0^536 

- 27.98535 

- 2.27923 

3 . 000 DO 

- 23.39534 

- 29.80534 

- 2.27923 

0.00000 

-33 . a 533 

- 31.62583 

- 2 . 2/923 

0.00000 

- 21.07505 

- 22.03981 

- 1.69731 

0.00000 

- 23.00457 

- 23.96934 

- 1.69731 

O.OOJOO 

- 24.93410 

- 25.89886 

- 1.69731 

0.00000 

- 26. 86363 

- 27 . 82839 

- 1.69731 

0.00000 

- 23.79315 

- 29.75792 

- 1.69731 

0.00000 

- 30. 72268 

- 31.68744 

- 1.69731 

0.00000 

- 19,64792 

- 21.24375 

- 1.46000 

0.00000 

- 22.33953 

- 24.43542 

- 1.46000 

0.00000 

- 26.03125 

- 27.62708 

- 1.46000 

0.00000 

- 29. 22292 

- 30.31875 

- 1.46000 

0.00000 

- 32-41458 

- 34.01042 

- 1.46000 

0.00000 

-3 5.6062 5 

- 37.20203 

- 1.46000 

0.00000 

- 19.64792 

- 21.24375 

- 1.31000 

0.00000 

- 22.83953 

- 24.43542 

- 1.31000 

0.00000 

- 26.03125 

- 27.62708 

- 1.31000 

0.00000 

- 29.22292 

- 30.31875 

- 1.31000 

0.00000 

- 32.41458 

- 34.01042 

- 1.31000 

0.00000 

- 35.60625 

- 37.20208 

- 1.31000 

0.00000 

- 19.64 792 

- 21.24375 

- 1.05000 

0.00000 

- 22.83958 

- 24.43542 

- 1.05000 

0.00000 

- 26.03125 

- 27.62708 

- 1.05000 

0.00000 

- 29.22292 

- 30.81875 

- 1.05000 

0.00000 

- 32.41458 

- 34.01042 

- 1. 05000 

0.00000 

- 35.60625 

- 37.20208 

- 1.05000 

0.00000 

- 19,64792 

- 21.24375 

-.72500 

0.00000 

- 22.83958 

- 24.43542 

-. 72500 

0.00000 

- 26.03125 

- 27.62708 

72500 

0.00000 

- 29.22292 

- 30.81875 

-.72500 

0.00000 

- 32.41453 

- 34.01042 

-.72500 

0.00000 

- 35.60625 

- 37.20208 

-.72 500 

0.00000 

- 19.64792 

- 21.24375 

-. 2/500 

0.00000 

- 22. 83958 

- 24.43542 

-. 2/500 

0.00000 

- 26.03125 

- 27.62708 

-.27500 

0-00000 

- 29.22292 

- 30.81875 

-.27500 

0.00000 

- 32.41458 

- 34.01042 

-.27500 

0.00000 

- 35.60625 

- 37.20208 

-.27500 

0.00000 


.38462 

42.90475 

0.00000 

0.00000 

.38462 

36.54497 

0.00000 

0.00000 

•36462 

28.94001 

0.00000 

0. 00000 

.38462 

20.03721 

0.00000 

0.00000 

.38462 

10.00798 

0.00000 

0.00000 

.38462 

-.67404 

0.00000 

0.00000 

.19731 

42.90475 

0.00000 

0.00000 

.19731 

36,54497 

0.00000 

0.00000 

.19731 

28.94001 

0. 00000 

0.00000 

.19731 

20.03721 

0.00000 

0.00000 

. 19731 

10.00798 

0.00000 

0.00000 

.19731 

-.67404 

0.00000 

0.00000 

. 04000 

0.00000 

0.00000 

0.00000 

•04000 

0.00000 

0.00000 

0.00000 

.04000 

0. 00000 

0.00000 

0.00000 

.04000 

0.00000 

0.00000 

0.00000 

.04000 

0.00000 

0.00000 

0.00000 

.04000 

0.00000 

0.00000 

0.00000 

.11000 

0.00000 

0.00000 

0.00000 

•11000 

0.00000 

0.00000 

0.00000 

•11000 

0.00000 

0.00000 

0.00000 

.11000 

0.00000 

0.00000 

0.00000 

.11000 

0.00000 

0.00000 

0.00000 

.11000 

0. 00000 

0.00000 

0.00000 

.15000 

0.00000 

0.00000 

0.00000 

•15000 

0.00000 

0.00000 

0.00000 

.15000 

0.00000 

0.00000 

0.00000 

.15000 

0.00000 

0.00000 

0.00000 

.15000 

0.00000 

0.00000 

0.00000 

•15000 

0.00000 

0.00000 

0.00000 

.17500 

0.00000 

0.00000 

0.00000 

.17500 

0.00000 

0.00000 

0.00000 

.17500 

0.00000 

0.00000 

0.00000 

.17500 

0.00000 

0.00000 

0.00000 

.17500 

0.00000 

0.00000 

0.00000 

.17500 

0.00000 

0.00000 

0.00000 

.27500 

0.00000 

0.00000 

0.00000 

-27500 

0.00000 

0.00000 

0.00000 

.27500 

0.00000 

0.00000 

0.00000 

.27500 

0.00000 

0.00000 

0.00000 

.27500 

0.00000 

0.00000 

0.00000 

•27500 

0.00000 

0.00000 

0.00000 


REF . CHORD 

C AVERAGE 

TRUE AREA 

REFERENCE AREA 

8/2 

REF . AR 

TRUE AR 

9.17560 

13.51260 

270.25200 

159,99696 

10.00000 

2.50005 

1.48010 


•80015 

•49U8 

.37956 

•33268 

.25925 

.13259 

.71316 

•52206 

•30060 

•35297 

.27534 

.09532 

.39401 

.37683 

•36901 

.21472 

•03357 

.00860 

•38667 

.40078 

.35211 

•21240 

.05136 

•01443 

•40181 

.40029 

•33914 

.21200 

•06333 

•01950 

.42072 

.39237 

•33066 

.21129 

•07126 

.02322 

.43575 

.38501 

.32522 

.21049 

•07638 

•02574 


NUMBER 

.30000 
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CJMPLFTe CONFIGURATION 

DESIRED CL COMPUTED ALPHA 
1,00000 17,90535 


WING-BODT CHARACTERISTICS 
l_IFT INDUCED DRAG (FAR FIELD SOLUTIGNJ 

ri (WBJ COl AT CL(^B) CO I / (CL ( H8 J **2 ) 

(1/(PI^AR) - .12732 I 

.5972^ .06733 ,ld876 


STATION 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


13 

14 

15 

16 
17 

la 

19 

20 
21 
22 

23 

24 

25 

26 


2Y/8 

-, 614^>5 
54165 
— , 463 1 6 
39586 
-.32296 
-.25007 
-.18181 
-.14600 
-.13100 
-.10500 
-.07250 
-.02^50 


-.96154 
-.88462 
-.80769 
-.73077 
-.67165 
-.61254 
-.53562 
-.45809 
-.381 77 
-.30485 
-.22792 
-.16973 
-. 14600 
-.13100 


complete CONFIGURATION CHARACTERISTICS 


CL 

PER RADIAN 
3. 19992 


ALPHA 

PER DEGREE 
.05585 


CLITWISTJ 

0.00000 


alpha at CL=0 T CP 
-0.00000 -.41849 


CM/CL 

-2.33969 


CMO 

0.000 00 


AODITIONAL LOADING 
UITH CL BASED ON S(TRUE) 

SL COGF CL RATIO C RATIO 


LOAD DUE ADD- LOAD AT BASIC LOAD 
TO TrtlST CL* 0.00000 AT CL*0 


SPAN LOAD AT 
DESIRED CL 


.29459 
.4241 7 
.51395 
.58240 
.63523 
.67490 
.70180 
.74890 
.75574 
.76401 
.77124 
.77651 


2. 18482 
2. 14999 
1.9/871 
1.80761 
1.65146 
1.50949 
1.38810 
.71015 
.67406 
.63924 
.60767 
.57494 


. 13483 
. 19729 
.25974 
• 32220 
.38465 
.44710 
. 50558 
1.05457 
1.12118 
1. 19518 
1.26919 
1.35059 


0.00000 
o.ooouo 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
Q. 00000 
0.00000 
O.OOOQO 


0.00000 

O-OOUOO 

0.00000 

0-00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


0.00000 
o.uoooo 
0. 00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 


,17440 

.25112 

.30428 

.34^80 

.37608 

.39956 

.41549 

.44337 

,44742 

-45231 

.45660 

.45972 


contribution of the second PLANFORM TO SPAN LOAD DISTRIBUTION 


39826 

04280 

.19496 

5842 7 

2.25367 

.25925 

71649 

2.21450 

.32355 

81880 

2.11118 

,39704 

83381 

2,02130 

.43725 

62500 

1.28426 

.48666 

56651 

1.02823 

.55095 

,54809 

.89085 

.61525 

.542 52 

♦79836 

.6/954 

,54273 

.72963 

,74384 

,54496 

.67435 

.80813 

,54492 

.63602 

.85677 

,55/25 

,39321 

1.41720 

.56563 

.39912 

1-41720 


0.00000 

0-00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 -00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0-00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0-00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


.23578 
.34590 
.42418 
.48475 
.52324 
.37002 
.33 539 
.32449 
-32119 
.32131 
.32263 
.32261 
.32991 
.33487 


-AT CL OES- 
X LOCATON OF 
LOCAL CENT PR 


-18.50420 
-17.74451 
-16.97449 
-16.20531 
-15.44203 
-14.6893 7 
-13.99U0 
-12.85910 
-11.99776 
-11.11893 
-10. 36793 
-9.6/063 


-28.82026 
-28.31349 
-27. 77968 
-27.22058 
-26. 75954 
-26.81144 
-26.44443 
-25.9891 7 
-25.51920 
-25.05822 
-24.61574 
-24.28605 
-24.07250 
-24. 19489 


CO 

09 
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NP 


27 -.10500 
20 -.07250 
29 -.02750 


• 57294 
.57830 
.58193 


.40428 

.40806 

.41062 


1.4I72J 

1.41720 

1.41720 


0.00000 

0.00000 

0.00000 


0.00000 

0.00030 

0.00000 


0.00000 

0.00000 

0.00000 


.33920 

-24 

.3423/ 

-24 

.34452 

-24 


INDUCED DRAG, LEADING EDGE THRUST AND SUCTION COEFFICIENT CHARACTERISTICS 
COMPUTED AT THE DESIRED CL FROM A NEAR FIELD SOLUTION 


STATION 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


SECTION COEFFICIENTS 


L. E. SWEEP 
2Y/B ANGLE 

-.61455 54.21355 

-.54165 54.21355 

-.46876 54.21355 

-.39586 54.21355 

-.32296 54.21355 

-.25007 54.21355 

-. 10101 54. 21355 

-.14600 84.28941 

-.13100 77. 59258 

-.10500 73.30076 

-.07250 70.70995 

-.02750 65.37644 


con C/28 

CT C/2B 

-.00351 

.02261 

-.00^50 

.03094 

-.00070 

.03272 

•00393 

.03234 

.00881 

.03082 

.01540 

.02672 

.02575 

•01890 

.02279 

.02402 

.02687 

.02025 

.03394 

.01370 

.04300 

.00512 

.05588 

-.00735 


CS C/26 
. 03866 
.05292 
.05595 
.05530 
.05271 
.04569 
.10663 
.23456 
.18165 
. 08992 
.01812 
-.01808 


CONTRIBUTION OF THE SECOND PLANFORM TO THE CHORD OR DRAG FORCE 


-.96154 
-.88462 
-.80769 
-.73077 
-.67165 
-.61254 
-.53562 
-.45869 
-.38177 
-.30485 
-.22792 
-. 16973 
-.14600 
-.13100 
-.10500 
-.07250 
-.02750 


44.31794 
44.31794 
44.31794 
44.31794 
44.31794 
44.31794 
44.31794 
44.31794 
44.31794 
44.31794 
44.31794 
44.31794 
0.00000 
0. 00000 
0.00000 
0.00000 
0.00000 


-.01214 
-.00759 
-.00567 
.00017 
.03605 
.06006 
.03478 
.02943 
.02890 
.02974 
.03063 
.033 76 
.03551 
.03483 
•03365 
.03295 
.03312 


.03563 
.04391 
.05030 
.05097 
.01768 
-.01999 
.00107 
.00496 
.00501 
• 00420 
«00341 
.00075 
-.00070 
.00029 
•00200 
.00314 
.00325 


.04979 
.06137 
.07031 
.07123 
.02471 
-.02 793 
.00149 
.00693 
.00700 
.00587 
.00476 
.00124 
-.00070 
.00029 
•00200 
.00314 
-00325 


TOTAL COEFFICIENTS 

C0II/CL**2 » .15956 CT= .15098 CS= .29453 

111 ' ' I 


23822 

24694 

24557 
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KP , KV AND RESPECTIVE CHOROWISE CENTROIDS FOR EACH PLANFORM 


KP- 
KV LE== 
KV $E- 


PLANFORM 
1 .2ttB79 
1.55287 
.22241 


NO. I 

CENTROID AT 
CENTROID AT 
CENTROID AT 


-14.39074 

-14.70269 

-19.26927 


KP = 
KV LE = 
KV SE= 


PLANFORM NO 
1.91113 
1.04260 
.45948 


. 2 

CENTROID 

CENTROID 

CENTROID 


AT -26.24069 
AT -26.77147 
AT -29.83309 


ISD 

CJI 
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■S3 


ALPHA 

0.0000 

2.0000 

4.0000 

6.0000 
8.0000 

10.0000 

12.0000 

14.0000 

16.0000 
18.0000 
20.0000 
22.0000 

24.0000 

26.0000 
28.0000 

30.0000 

32.0000 

34.0000 

36.0000 

38.0000 

40.0000 

42.0000 

44.0000 

46.0000 

48.0000 

50.0000 


PERFORMANCE CHARACTERISTICS FOR PLANFORM 1 


CN 

0.0000 
.0471 
.0983 
.1534 
.2120 
.2739 
.3388 
.4064 
.4764 
.5483 
.6219 
. 6 968 
.7 726 
.8489 
.9255 
1-0019 
1.0777 
1.1526 
1.2262 
1.2982 
1.3681 
1 .4357 
1.5007 
1.5626 
1.6213 
1.6764 


CLP CLPi-CLVLE CLP^-CLVSE CL 


CMP CMPKMVLE CMP>CMVSE CM 


CO CL««2/(PI«ARi 


0.0000 

0.0000 

.0449 

.0463 

.0895 

.0970 

.1332 

.1501 

.1759 

.2057 

.2170 

.2632 

.2564 

.3220 

.2935 

.3817 

.3282 

.4417 

.3602 

.5013 

-3892 

.5599 

.4150 

.6171 

.4375 

.6722 

.4564 

.7246 

.4717 

. 7739 

.4833 

.8195 

• 4912 

. 8610 

.4953 

.89 79 

.4958 

.9299 

.4927 

.9565 

.4861 

.9776 

.4763 

.9929 

.4633 

1.0023 

.4474 

1.0055 

.4288 

1.0027 

.40/9 

.9937 


0.0000 

0.0000 

.0452 

.0471 

.0905 

• 0981 

.1357 

.1525 

• 1802 

.2099 

.2237 

.2698 

.2658 

.3314 

• 3062 

.3944 

• 3445 

.4579 

.3804 

.5215 

.4137 

.5844 

.4440 

.6460 

.4711 

.7058 

.4948 

.7630 

.5150 

.8172 

.5315 

.8677 

.5441 

.9139 

.5530 

.9555 

.5580 

.9920 

.5591 

1.0230 

.5565 

1.0480 

.5503 

1.0669 

.5405 

1.0795 

.52 73 

1.0855 

.5110 

1.0849 

.4918 

1.0776 


0.0000 

0.0000 

-.0705 

-.0735 

-.1407 

-.1528 

-.2101 

-.2373 

— ♦ 2 786 

-.3268 

-.3457 

-.4207 

-•4111 

-.5186 

-.4745 

-.6201 

-.5356 

-.7246 

-.5940 

-.8317 

-.6496 

-.9407 

- . 702 1 

-1.0512 

-.7511 

-1.1627 

-.7964 

-1.2746 

-.8379 

-1.3863 

-.8753 

-1.4973 

-.9084 

-1.6071 

-.9371 

-1.7151 

-.9612 

-1.8209 

-.9806 

-1.9238 

-.9953 

-2.0234 

-1.0051 

-2.1192 

-1.0100 

-2.2107 

-1.0100 

-2.2976 

-1.0051 

-2.3793 

-.9953 

-2.4555 


0.0000 

0. 0000 

-.0711 

-.0741 

-.1429 

-.1550 

-.2152 

-.2424 

-.2876 

-.3358 

-.3597 

-.4348 

-.4313 

-.5388 

-.5018 

-.6474 

-.5711 

-. 7601 

-.6386 

-.8763 

-.7043 

-.9953 

-.7676 

-1.1166 

-.8283 

-1.2400 

-.8862 

-1.3643 

-.9408 

-1.4892 

-.9920 

-1.6141 

-1.0395 

-1.7383 

-1.0831 

-1.8612 

-1.1226 

-1.9822 

-1.1577 

-2.1008 

-1.1883 

-2.2164 

-1.2142 

-2.3283 

-1-2354 

-2.4361 

-1.2517 

-2.5393 

-1.2631 

-2.6372 

-1.2694 

-2.7296 


0.0000 

0.0000 

• 0016 

.0003 

• 0069 

.0012 

.0160 

• 0030 

.029 5 

.0056 

.0476 

• 0093 

.0704 

• 0140 

.0983 

• 0198 

• 1313 

.0267 

.1694 

• 0346 

.2127 

.0435 

.2610 

.0531 

.3142 

.0634 

.3722 

.0741 

.4345 

.0850 

.5009 

.0959 

.5711 

.1064 

.6445 

• 1163 

.7207 

.1253 

.7992 

.1332 

.8794 

.1398 

.9607 

.1449 

1.0424 

.1484 

1.1241 

.1500 

1.2049 

• 1498 

1.2642 

.1478 
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ALPHA 

0,0000 
2.0000 
A. 0000 
6.0000 
6. 0000 
10.0000 
12.0000 

14.0000 

16.0000 
18.0000 
20.0000 
22.0000 

24.0000 

26.0000 
28.0000 

30.0000 

32.0000 

34.0000 

36.0000 

38.0000 

40.0000 

42.0000 

44.0000 

46.0000 

48.0000 

60.0000 


PERFORMANCE CHARACTERISTICS FOR PLANFORM 2 


CN 


CLP CLP^CLVLE CLP+CLVSE CL 


CMP CHP^CMVLE CMP+CMVSE CM 


CO CL**2/<PI*ARJ 


0.0000 

0.0000 

.0685 

.0666 

.1403 

, 1327 

.2151 

.1976 

.292 5 

.2608 

.3721 

.3219 

.4536 

.3802 

• 5365 

.4353 

.6205 

.4868 

. 7051 

. 5342 

. 7899 

.5772 

.8 746 

.6155 

. 9586 

.6487 

1.0417 

.6768 

1.1233 

• 6995 

1.2031 

.7167 

1.2807 

.7284 

1 .3557 

• 7345 

1.42 78 

.7352 

1.4965 

-7306 

1.5617 

.7209 

1.6229 

.7062 

1.6798 

.6870 

1. 7322 

.6634 

I.T799 

.6359 

1.8225 

.6049 


0.0000 

0.0000 

.06 79 

• 0672 

.1377 

.1349 

.20 89 

.2026 

.2803 

.2696 

.3528 

.3355 

.4243 

-3996 

« 4945 

.4614 

.5629 

.5203 

.6289 

.5759 

• 6918 

.6277 

.7511 

.6752 

• 8063 

.7182 

• 8569 

.7562 

.9024 

.7889 

.9424 

.8162 

.9766 

.8378 

1.0048 

.8536 

1.0266 

.863 7 

1-0420 

.8679 

1.0509 

• 8663 

1.0531 

.8591 

1.0489 

.8465 

1.0382 

.8285 

1.0212 

.8057 

.9982 

.7782 


0.0000 

0.0000 

•0684 

-.1906 

.1400 

-.3803 

.2139 

-.5682 

.2896 

-.7533 

.3665 

-.9347 

.443 7 

-1.1115 

.5206 

-1.2830 

.596 5 

-1.4481 

• 6706 

-1.6063 

.742 3 

-1.7566 

.8109 

-1.8983 

.8757 

-2.0308 

.9362 

-2.1534 

.9918 

-2.2656 

1.0419 

-2.3666 

1.0861 

-2.4562 

1.1239 

-2.5338 

1.1551 

-2.5990 

1.1793 

-2.6516 

1.1963 

-2.6912 

1.2060 

-2.7178 

1.2084 

-2. 7311 

1.2033 

-2.7311 

1.1910 

-2.7178 

1.1715 

-2.6912 


0.0000 

0.0000 

-.1943 

-.1924 

-.3951 

-.3876 

-.6014 

-.5845 

-.8122 

-.7822 

-1.0264 

-.9797 

-1.2430 

-1. 1761 

-1.4610 

-1.3704 

-1.6793 

-1.5616 

-1.8968 

-1.7489 

-2.1124 

-1.9313 

-2.3252 

-2.1080 

-2.5341 

-2.2780 

-2.7380 

-2.4405 

-2.9360 

-2.5948 

-3.1271 

-2.7401 

-3.3104 

-2.8757 

-3.4850 

-3.0009 

-3.6500 

-3.1152 

-3.8046 

-3.2178 

-3.9481 

-3.3085 

-4.0798 

-3.3867 

-4.1990 

-3.4520 

-4.3052 

-3.5041 

-4.3978 

-3.5428 

-4.4764 

-3.5679 


0.0300 
1962 
-.4024 
-.6177 
-. 841 1 
-1.0714 
-1.3076 
-1.5484 
-1.7928 
-2 . 0394 
-2.2872 
-2.5349 
-2.7812 
-3.0251 
-3.2653 
-3.5006 
-3.7299 
-3.9521 
-4.1661 
-4.3709 
-4.5654 
-4.7487 
-4.9199 
-5.0782 
-5.2228 
-5.3530 


0.0000 

.0024 

.0098 

.0225 

.0407 

.0646 

.0943 

.1298 

.1710 

.2179 

.2702 

.3276 

• 3899 
.4566 
.5273 
.6015 

• 6786 
.7581 
.8392 
.9214 

1.0038 

1.0859 

1.1669 

1.2461 

1.3227 

1.3961 


0.0000 
.0006 
.002 5 
.0058 
.0107 
.0171 
.0251 
•0345 

• 0453 
.0573 
.0702 
•0837 
.0976 

• 1116 
.1252 

• 1382 
.1502 

• 1608 
.1699 
.1771 
.1822 
• 1852 
.1859 
.1844 
.1806 
.1747 


ts9 

-a 



00 


TOTAL PERFORMANCfc CHARACTERISTICS 


ALPHA 

CN 

CLP 

CLP+CLVLE 

CLP<-CLVSE 

CL 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

2.0000 

.1156 

.1115 

.1147 

.1124 

.1155 

4.0000 

.2386 

.2221 

.2347 

.2254 

.2380 

6. 0000 

.3685 

.3308 

.3 590 

.3382 

.3664 

8.0000 

.5045 

.4367 

.4365 

• 4498 

.4996 

10.0000 

.6460 

.5389 

.6160 

.5592 

.6362 

12.0000 

.79 24 

.6365 

. 7463 

• 6654 

. 7751 

14.0000 

.9429 

.7288 

.8 762 

.7675 

.9149 

16.0000 

1.0969 

.8150 

1.0046 

.8648 

1.0544 

18.0000 

1.2534 

.8944 

1 .1301 

.9563 

1.1920 

20.0000 

1.4118 

.9664 

1.2517 

1.0414 

1.326 7 

22.0000 

1.5713 

1.0305 

1.3682 

1.1192 

1.4569 

24.0000 

1.7312 

1.0862 

1.4785 

1.1893 

1.5815 

26.0000 

1.8906 

1.1332 

1.5815 

1.2510 

1.699 3 

28.0000 

2.0488 

1.1712 

1.6763 

1.3039 

1.8090 

30.0000 

2.2049 

1.2000 

1.7619 

1-3476 

1.9095 

32.0000 

2.3584 

1.2195 

1.8376 

1.3819 

2.0000 

34.0000 

2.5083 

1.2298 

1.9027 

1.4066 

2.0795 

36.0000 

2.6540 

1.2310 

1.9565 

1.4216 

2.1471 

38.0000 

2.7947 

1.2233 

1.9986 

1.42 70 

2.2022 

40.0000 

2.9298 

1.2070 

2.0285 

1.4228 

2.2443 

42.0000 

3.0586 

1-1825 

2. 0461 

1. 4094 

2.2730 

44.0000 

3.1805 

1.1502 

2.0511 

1.3869 

2.2678 

46.0000 

3.2949 

1.1108 

2.0437 

1.3559 

2.2888 

48.0000 

3.4012 

1.0647 

2.0238 

1.3167 

2.2758 

50.0000 

3.4989 

1.0128 

1.9918 

1.2700 

2.2490 


CMP 

CHP4-CMVLE 

CMP^CMVSE 

CM 

CO CL*^2/IPI*ARI 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-.2611 

-.2679 

-.2635 

-.2703 

.0040 

.0017 

-.5210 

-.5479 

-.5305 

-.5574 

.0166 

.0072 

-.7783 

-.8387 

-.7997 

-.8602 

.0385 

.0171 

-1.0318 

-1.1389 

-1.0698 

-1.1769 

.0702 

-0318 

-1.2803 

-1.4471 

-1.3395 

-1.5062 

.1122 

.0515 

-1.5226 

-1.7616 

-1.6074 

-1.8464 

.1648 

.0765 

-1.7574 

-2.0811 

-1.8722 

-2.1959 

.2281 

• 1066 

-1.9837 

-2.4039 

-2.1327 

-2.5529 

.3023 

.1415 

-2.2003 

-2.7284 

-2.3876 

-2.9157 

.3873 

• 1809 

-2.4062 

-3.0531 

-2.6356 

-3.2825 

.4829 

.2241 

-2.6004 

-3.3765 

-2.8756 

-3.6516 

• 5886 

.2703 

-2.7819 

-3.6968 

-3.1063 

-4.0212 

.7041 

.3185 

-2.9498 

*4.0126 

-3.3267 

-4.3894 

.8288 

.3676 

-3.1034 

-4.3223 

-3.5356 

-4.7545 

.9618 

.4166 

-3.2419 

-4.6245 

-3.7321 

-5.1147 

1.1025 

.4643 

-3.3646 

-4.9175 

-3.9152 

-5.4682 

1.2497 

.5093 

*3.4708 

-5.2001 

-4.0840 

-5.8133 

1.4026 

.5506 

-3.5602 

-5.4709 

-4.2377 

-6. 1484 

1.5600 

.5870 

-3.6322 

-5.7284 

-4.3755 

-6.4717 

1.7206 

.6175 

-3.6865 

-5.9715 

-4.4968 

-6. 7818 

1.8832 

.6413 

-3.7229 

“6.1990 

-4.6009 

-7.0770 

2.0466 

.6578 

-3.7411 

-6.4098 

-4.6874 

-7.3560 

2.2093 

.6664 

-3.7411 

-6.6028 

-4.7559 

-7.6175 

2.3701 

.6670 

-3. 7229 

-6.7771 

-4.8059 

-7.8601 

2.5276 

.6594 

-3.6865 

-6.9318 

-4.8373 

-8.0826 

2.6803 

.6440 


THIS CASE IS FINISHED 


APPENDIX B 



APPENDIX C 


FORTRAN PROGRAM LISTING 

This program was written in FORTRAN IV language, version 2.3, for the Control 
Data series 6000 computer systems with SCOPE 3.0 operating system and library tape. 
Minor modifications may be required prior to use with other computers. The program 
requires 53000a words of storage on the Control Data 6600 computer system and consists 
of a main program, five overlays, and five subroutines. Each program or subroutine is 
identified in columns 73 to 75 by a 3 -character identification. In addition, each of these 
parts is sequenced with a 4 -digit number in columns 76 to 79. The following table is an 
index to the program listing: 


Program or subroutine 

Identification 

Page 

WINGAL 

MAI 

30 

INFSUB 

INF 

31 

GEOMTRY 

GEO 

32 

MATXSOL 

MAT 

42 

SSLESO 

SSL 

43 

AERODYN 

AER 

44 

FTLUP 

TLU 

49 

CDICLS 

CDI 

51 

CDRAGNF 

DRA 

53 

FTLUP 

TLU 

57 

TIPSUCT 

TIP 

59 

WRTANS 

WRT 

65 


The execution time for this program is similar to that for the program presented 
in reference 6. 
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c 

c 

c 

c 

c 


10 


20 


30 

40 

so 


60 

70 

C 

c 

c 

BO 

90 

100 


OVERLAY (WINGTL»OtO) s,Ai 

program *nM6AL ( INPJT=201 *OUTPUT=100l ♦TAPES=1NPUT»TAPE6=0 UTpuT*TApEmAI 
1 10=401 ) m j 

COMMON /all/ BOT .M,BETA,PTEST,QTEST,T8lSCw ( 50) »0(200) *PN (2oO ) . P V < 2M A I 
100) » ALP (200) ,S (200) .PS I (200) .PHl (50 ) .Zh(50) MAI 

COMMON /TOTHRE/ CIP(200,2) MAI 

COMMON /THREEOR/ CC A V ( 2 . 50 ) . CL T . CLNT »NSSW. ALPO MAI 

COMMON /ONETHRE/ T V I ST ( 2 ) . CREF , SPEF. CAVE .CEDES .STRUE . AR. ARTRUE . RTCM A I 
IDHT (2) .CONFIG.NSSWSV (2) ,MSV(2) .KBOT .PL AN. TPL AN . MACH, SSWW A ( 50 ) . XL ( 2MA I 
2) .XT (2) .CLX8.CMCL.CLA (2) .BLAIR (50) .CLAMAR(?) MAI 

COMMON /MAINONE/ I CODEOF . TOT AL . A AN ( 2 ) . XS ( 2 ) .YS(2) .KFCT5 ( 2 ) . XREG ( 25m A I 
1.2) .YREG(25.2) .AREG(25.2) .DIH(25.?).MC0( 25,2) .XX (25.2) .YY(25.2) . ASmAi 


2(25.2) .TTWD( 25.2) .MMCD(25.2) .An (2) . ZZ ( 25. 2 ) . I T IPCOO 
COMMON /CCRRDD/ TSPAN. TSPANA .K8 I T 

VORTEX lattice AERODYNAMIC COMPUTATION 
NASA-LRC PROGRAM NO, A2794 


icodeof=total=o 

WINGTL=6LWINGTL 

RECALL=6HRECALL 

call, overlay (WINGTL. 1.0. recall) 
if ( ICOOEOF,GT.O) go to 70 
IF (M.GT.200) GO TO 40 
NSW=NSSWSV ( 1 ) .NSSWSV (2) 

IF (NSX,GT.50) GO TO 30 
ITSV=0 

DO 20 IT=1.IPLAN 
IF (AN(IT) ,LE,25.) GO TO 
WRITE (6.100) IT.AN(IT) 

ITSV=1 
CONTINUE 

IF (ITSV,GT,0) GO TO 60 
GO TO 50 

NSW 


20 


WRITE 
GO TO 

write 

GO TO 

CALL 

CALL 


(6.90) 

60 

(6. BO) M 
60 


Overlay (wingtl. 2 . 0 . recall) 

OVERLAY (WINGTL. 3 . 0 . RECALL) 

IF (PTEST,EQ.1.,0R,QTEST,EQ,1,) go to 60 
call overlay (WINGTL. 4.0. RECALL) 

IF (ITIPCOD.EO.I) CALL OVERLAY ( W INGTL . 5. 0 .RECALL ) 
TOTAL=TOTAL-1 , 

GO TO 10 
STOP 


( 1H1//10X. I6.93HH0RSESH0E VORTICES LAIOOUT. THIS IS MORE 
200 maximum, this CONFIGURATION IS ABORTED,) 


FORMAT 

ian the 

FORMAT (1h1//10X. I6.101H ROWS OF HORSESHOE VORTICES LAIOOUT. 
IIS MORE Than the so maximum, this configuration IS aborted,) 
format ( 1H1//10X,8HPlANFORM. I6,4H HAS.I6.74H BREAKPOINTS. THE 
IMUM dimensioned IS 25. THE CONFIGURATION IS ABORTED.) 

END 


MA I 
MA I 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
mA I 
MA I 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
MAI 
THMAI 
MAI 


THIS MAI 
MAI 

maximai 

MAI 

MAI 


10 

20 

25 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

ISO 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

soo 

510 
S20 
530 
54 0 
550 
560- 
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SUttROUTIME InFSUB ( BOT » F V I * E w I ) 



INF 

10 


COMMON /1NSUB23/ PS I I . APHI I » XXX * YYY, ZZZ 

tSNNtTOLRNC 

INF 

20 


FC=COS(PSII) 



INF 

30 


FS=SIN(PSI I) 



INF 

40 


FT=F5/FC 



INF 

50 

c 




INF 

60 

c 




INF 

70 


FPC=COS(APhI I ) 



INF 

80 


FPS=SlN(APriI 1 ) 



INF 

90 


FPT=FPS/f;PC 



INF 

10 0 


F 1=XXX*SNN*FT*FPC 



INF 

110 


F2=YYY*SNN»FPC 



INF 

120 


F3=ZZZ*5NN*FPS 



INF 

130 


F4=XXX-SNN*FT*FPC 



INF 

140 


F5=YYY-SNN»FPC 



INF 

150 


F6=ZZZ-SNN*FPS 



INF 

160 


FFA= (XXX"«2* ( YYY*F»S) **2»FPC*«2* ( (YYY*FT)**2* 

(ZZZ/FC)**2-2.*XXX» 

rr INF 

170 


1 Y*FT)-2."ZZZ*FPC* <YYY*FPS*XXX*FT*FPS)) 



INF 

IdO 


FFB=(FI*F1*F2*F2*F3*F3) »*.S 



INF 

190 


FFC=(FA*FA*F5*F5*F6*F6)«*.5 



INF 

200 


FFD=F5*F5*F6*F6 



INF 

210 


FFE=F2*F2*F3*F3 



INF 

220 


FFF= (F1«FpC*FT*F2*FPC*F3*FpS)/FF8- <FA*FPC*FT* 

Fb*FPC*F6*FPS)/FFC 

INF 

230 

c 




INF 

240 

c 




INF 

250 

c 

THE TOLEPANCE SET AT THIS POINT IN THE 

PROGRAM MAY NEED TO BE 

INF 

260 

c 

changed FOS computers other than the COC 6000 

SEKIES 

INF 

270 

c 




INF 

280 

c 




INF 

290 


IF (A8S(^FA) .LT. (B0T*15.E-5)**2) GO TO 

10 


INF 

300 


FV0NE=(XXX*FPS-ZZZ*FT»FPC»*FFF/FFA 



INF 

310 


FW0NE=(YYY*FT-XXX) *FFF/FFA*FPC 



INF 

320 


GO TO 20 



INF 

330 

10 

fvone=fmone=o. 



INF 

340 

c 




INF 

350 

20 

IF (ABS(FFO) .LT.TOlRnC) GO TO 30 



INF 

360 


FVTW0=F6* ( 1 .-F4/FFC) /FFO 



INF 

370 


FWTWO=-F»* (1 .-F4/FFC) /FFD 



INF 

380 


GO TO 40 



INF 

390 

30 

FVTWO=Fi»TwO = 0. 



INF 

400 

C 




INF 

410 

40 

IF (ABS(FFE) .LT.TOLRNC) GO TO 50 



INF 

420 


F VTHRE=-F 3* n .-FI/FFB) /FFE 



INF 

430 


FWTHRE=F2*( 1 .-F I /FFB ) /FFE 



INF 

440 


GO TO 60 



INF 

450 

50 

fvthre=f»thre=o. 



INF 

460 

C 




INF 

470 

60 

FVI=FVONt*FVTWO*FVTHRE 



INF 

480 


FWI=FWONE*FWTWO*FWTHHE 



INF 

490 


RETURN 



INF 

500 


END 



INF 

510- 


31 



ooooooooo*- ooooooooo ooooooo 


APPENDIX C 


OVERLAY(»<INGTL.1«0) GEO 

PROGRAM GEOMTRY GEO 

dimension XREF(25». yREF( 2S). SAR(25)» A(25)t RSAR(25)* X<2S). Y(2 gE 0 
15)» B0TSV(2)» GA(2). VBORD(Sl). SPY(50«2). KFX(2)* IYU50*2)* lYjtGEO 
2S0.2) GEO 

COMMON /All/ B0T.M»BETA,PTEST.QTEST,TBLSC»y(50) ,Q(200) .PN(2oO) #PV(2GE0 
100>»ALP(200),S(200).PSI(200 ).PhI(50)»2h(50) GEO 

COMMON /ONETHRE/ TNIST (2) t CREF , SREF. C A YE t CLDES . STRUE » AR , ARTRUE . RTCGEO 
10HT<2) .C0NFIG«NSSUSV (2) ,MSV(2 ) ,KBOT «PL AN » IPL AN » M ACh. 5SWW A ( 50 ) *XL(2GE0 
2) tXT(2) ,CLW8,CMCL*CLA(2) «BLAIR(50) .CLAMAR(2) GEO 

COMMON /MAINONF/ I COOEOF . TOT Al . A An ( 2) t XS ( 2 ) , YS ( 2 ) , <FCTS ( 2 ) ,XREG(2SGE0 
1 f 2) ,YREG(25,2) tAREG(25.2) .DIH(25.2) ,MC0 ( 25 . 2 ) t XX ( 25 , 2 ) » Y Y ( 25, 2 ) ,ASGE0 


2(25,2> ,TTirfO(25,2) ,MMCO<25,2 ) , An ( 2 ) ,22 ( 25,2 ) , I T IPCOD 
COMMON /CCRRDO/ TS° An , TSP AN A , KB I T 
REAL MACH 


PART ONE - GEOMETRY COMPUTATION 


SECTION ONE - INPUT OF REFERENCE WING POSITION 


IF (TOTAL. EO.O 
YTOL=l .E-10 
A2Y=1 ,E*I3 
PIT=1 .5707963 
RAD=57.2957B 
IF (TOTAL. GT.O 


set plan Equal 

VARIABLE sweep 
SET plan equal 


,0) RTCOHT ( 1 ) =RTCOHT (2)=XL (2)=XT (2)=0.0 


) GO TO 70 


TO 1 . 
WING 
TO 2. 


FOR A wInG AlONE COMPUTAION - EVEN FOR A 
FOR A WING - TAIL COMBINATION 


SET total equal TO ThE NUMBER OF SETS 
OF GROUP TWO DATA PROVIDED 

READ (5,BB0) plan, total, CREF.SREF 
IF (ENDFILE 5) 830,10 
IPLAN=PL*N 


SET AAN(IT) EQUAL TO THE MAXIMUM NUMBER OF CURVES REQUIRED TO 
DEFINE THE PLANFORM PERIMETER OF THE <IT) PLANFORM. 

SET RTCDHT(IT) EQUAL TO THE ROOT CHORD HEIGHT OF THE LIFTING 

Surface ut), whose perimeter points are being read in, with 

RESPECT TO THE WING ROOT CHORD HEIGHT 

WRITE (6*860) 

DO 60 IT=1,IPLAN 

READ (5,880) AAN(IT) ,XS(IT) ,YS(IT),RTCOHT(IT) 

N=AAN(1T1 
N1=N*1 
MAK = 0 

IF (IPLAN.EQ.1) PRTCON=10H 

IF (IPLAN.E0.2.AN0.IT.EQ.1) PRTCON=IOH FIRST 
IF (IPLAN.EQ.2.AND.IT.E0.2) PRTCON=10H SECOND 

WRITE (6*870) PRTC0N,N*RTCDHT(IT),XS(IT)*YS(IT) 

WRITE (6*990) 


GEO 

geo 

GEO 

geo 

geo 

GEO 

GEO 

geo 

GEO 

GEO 

GEO 

GEO 

geo 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

no 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 


geo 260 

GEO 270 

geo 280 

GEO 290 
GEO 300 
GEo 310 
GEO 320 
GEO 330 
geo 340 
GEO 350 
GEO 360 
geo 370 

geo 380 

geo 390 
GEO 400 
GEO 410 
geo 420 
geo 430 
GEO 440 
geo 450 
GEO 460 
GEO 470 
GEO 480 
geo 490 

geo 500 

geo 510 
GEO 520 


GEO 

GEO 

GEO 

GEO 

GEO 

GEO 

geo 

geo 


530 

5A0 

550 

560 

570 

580 

590 

600 
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20 


30 


40 


50 


60 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

70 


BO 


00 50 I=1*N1 

READ (5*680) XREG ( I ♦ I T ) ♦ YREG ( I * I T ) *DlH (I * I T ) ♦ AMCO 

MC0(I*1T)=AMCD 

IF (I.EO.l) GO TO 50 

IF (MAK,NE.O.OR.MCD( I-l * IT) .NE.2) GO TO 20 
MAK=I-1 

IF <A3S(TREG(I-1,IT)-YREG(I.IT) ) .LT.YTOL) GO TO 30 

AREG(I-1*IT)=(XREG(I-1,IT)-XREG(I*IT) )/(YHEG<I-l.IT)-YkEG(I.in ) 

ASWP=ATA4IAREG( I-1*IT) )*RA0 
GO TO 40 

YREG(I*IT)=YREG(I-1*IT) 

AREG(I-1*IT)=AZY 

ASnP=90. 

J=I-I 


geo 610 
geo 620 

geo 630 
geo 640 
geo 650 

geo 660 

geo 670 
GEO 650 
GEO 690 


•(RITE PLANFORM perimeter POINTS AND ANGLES 

WRITE (6*960) J , XREG ( J . I T > ♦ YREG ( J . I T ) . ASWP , 0 IH ( J, I T ) *MCD ( J ♦ I T ) 
0IH(J*IT)=TAN(0IH(J*IT)/RAD) 

CONTINUE 
KFCTS( IT)=MAK 

write (6*960) Nl *XREG(N1 * IT) *YREG(N1« IT) 

CONTINUE 


READ GROUP 2 


PART I - section 2 

DATA and compute DESIRED WING POSITION 


SET SA(1)*SA(2) EQUAL TO THE SWEEP ANGLE* IN DEGREES, FOR T HE 
CURVE(S) THAT CAN CHANGE SWEEP FOR EACH PLANFORM 


IF A PARIICULAR 
TO BE COMPUTED, 
SET CLDES EQUAL 


IF 

IF 

DO 


PTEST IS SET 
QTEST IS SET 
NOT SET both 


SET 

FOR 


TWIST(l) OR 
A PLANFORM 


TWIST(2) equal to 0. 
THAT HAS TWIST AND/OR 


geo 700 
GEO 710 
GEO 720 
GEO 730 
geo 740 
GEO 750 
geo 760 
GEO 770 
geo 750 
GEO 790 

geo 800 
geo 810 

GEO 820 
GEO 830 
GEO 840 
geo 850 

geo 860 

GEO 870 
GEO 850 
FIRSTGEO 890 
geo 900 
GEO 910 
ARE geo 920 
geo 930 
TO I.OGEO 940 
geo 950 

COMPUTE ClP geo 960 

COMPUTE CMQ AND CLQGEO 970 
SINGLE configuration GEO 980 

geo 990 

FOR A FLAT PLANFORM AnD TO 

camber 


value of cl is desired at which the loadings 

SET CLDES EQUAL TO THIS VALUE 
TO 11. FOR A DRAG POLAR AT CL VALUES OF-.l 


equal to one the program will 
equal to one the program will 

PTEST AND QTEST TO ONE FOR A 


SET ATPCOD TO ONE IF THE CONTRIBUTIONS TO LIFT, DRAG AND MOmENI 
FROM SEPERATEO FLOW AROUND THE LEADING ANO/OR SIDE EDGES IS 
DESIRED. OTHERWISE SET ATPCOD TO ZERO. 


1 .GEOIOOO 
GEOlOlO 
GE01020 
GEOI030 
GE01040 
GE01050 
GE01060 
GE01070 

READ (5,950) CONE I G, SCW , V I C *’M ACH, CLOES.PTEST ,QTEST , TW I ST ( 1 ) , SA ( I ) *^^01050 
1TWIST(2) *SA(2) .ATPCOD 
ITIPCOD=ATPCOO 
IF (ITIPCOD.NE.l) GO TO 110 
DO 100 IT=1*IPLAN 
N8B6=AAN(IT) 

DO 90 IB8G=2,NBBG 

( YREG( IBBG, IT ) .EQ.YREG ( IBBG* 1 * IT) ) GO TO 80 
TO 90 

( YREG( I8HG*2, IT) .LT . YREG( I80G*1 , IT) ) GO TO 90 
(YREG(IBBG-1*IT).LT.YREG(I88G,IT)) go TO 90 
XL(IT)=XREG(I8BG,IT) 

XT(IT)=XREG(I8BG*1*IT) 


IF 

GO 

IF 

IF 


GE01090 
GEOI 100 
GEOlllO 
GEOI 120 
GEOI 130 
GEOI 140 
GEOH50 
GE01160 
GE01170 
GEOI 180 
GEO1190 
GE01200 
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GO TO 100 


GE01210 

90 

CONTINUE 


GE01220 


Xl(IT)=0.0 


GEO1230 


XT(IT)=0.0 


GE012‘*0 

100 

CONTINUE 


GE01250 

110 

CONTINUE 


GE01260 


WRITE (6*890) CONFIG 


GE01270 


IF (ENDFILE 5) 830fl20 


GEO1280 

120 

IF (PTEST.NE.O. .AND.0TE5T.NE.0. ) GO 

TO 850 

GE01290 


IF (SCW.tQ.O.) GO TO UO 


GE01300 


00 130 I=l»50 


GEO1310 

130 

T8LSCW(I)=SCw 


GE01320 


GO TO 150 


GEO 1330 

140 

READ (5.880) STA 


GEOl 340 


NSTA=STA 


GEO13S0 


READ (5.880) (TBLSCW ( I ) . TBLSCW ( I ♦ 1 ) 

.TBLSCW( 1*2) .T8LSCW( I*J) .THLSCWGEOl3<iO 


1 (!♦<.) .TBLSCW (1*5) .TBlSCw (1*6) .TBLSCW(I«7) . 1=1 .NSTA.8) 

GE01370 

150 

DO 410 IT=1.IPlAN 


GE01380 


N=AAN(IT) 


GEOl 3'90 


N1=N*1 


GEOUOO 


DO 160 1=1. N 


GE01410 


XREF(I)=XRtG(I.IT) 


GE01420 


YREF(I)=YR£G(1,IT) 


GEO1430 


A(n=AREG( I.IT) 


GE01A40 


RSAR(I)=ATAN(A(I) ) 


GEOIaSO 


IF (A(I) .EO.AZY) R5AR(I)=PIT 


GEO l4b0 

160 

CONTINUE 


6E01470 


XREF(N1)=XREG(N1.IT) 


GEOU80 


YREF(N1)=YREG(N1.IT) 


GEOU90 


IF (KFCT5 ( IT) .GT.O) GO TO 170 


GE01500 


K = 1 


GEO1510 


SA(IT)=RSAR(l)*RAO 


GEO1520 


GO TO 180 


GEO1S30 

170 

K=KFCTS(IT) 


GEO1540 

IBO 

WRITE (6.920) K.SA(IT).IT 


GEO1S50 


SB=SA( IT)/RAD 


GE01S60 


IF (ABS(b8-RSAR(K) ) .GT. (. 1/RAO) ) GO 

TO 210 

GE01570 

C 

REFERENCE PLANFORM COORDINATES ARE STORED UNCrtANGEO FOR WiNGS 

GE01S80 

c 

WITHOUT CHANGE IN SWEEP 


GE01S90 


DO 200 1=1. N 


GE01600 


X(I)=XREF(I) 


GE01610 


Y(I)=YREF(I) 


GEO1620 


IF (RSAR(I).EQ.PIT) GO TO 190 


GE01630 


A(I)=TAN(RSAR(I) ) 


GE01640 


GO TO 200 


GE01650 

190 

A(I)=AZY 


GE01660 

200 

SAR ( I ) =RSAR ( I ) 


GEO 167 0 


X(N1)=XREF(N1) 


GE01680 


Y(N1)=YREF(N1 ) 


GE01690 


GO TO 390 


GEOUOO 

C 



GEO1710 

c 

changes In wing sweep are made here 


GE01720 

c 



GEOU30 

210 

IF (HCD(X.IT) .NE.2) GO TO 840 


GEO 1740 


KA=K-1 


GE01750 


DO 220 1=1. KA 


GEOU60 


X(I)=XREF(I) 


GE01770 


Y(I)=YREF (I) 


GEOU80 

220 

SAR(I)=RSAR(I) 


GE01790 

C 

DETERMINE LEADING EDGE INTERSECTION 

between fixed and variable 

GE01800 
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c 


c 

d30 


^40 

C 

c 

dbO 

2b0 


270 


2ti0 

290 

300 

310 

320 

330 

3A0 


C 

c 

3b0 


360 

370 

C 

c 


380 

390 


SWEEP WING SECTIONS 0^18^0 

SAR(K)=5C rFnlftSO 

A(K)=TANISB) 6EO1840 

MK*n=XSnTM(XREF(K*l)-XS(IT) )*COS(SAl)MYRtF ( K * 1 ) -YS ( I T ) ) *S IN ( SGEO 1 850 

GtOlQbU 

Y(K*1)=YS(IT)MYREF(K*1)-YS(IT) )*COS(SAI)-(XRtF(K»l)-XS(lT) )*5IN(SGE01870 


lAI ) 

IF ( A8S (SB-SAR (K-1 )) ,LT .(. 1/PAD) ) GO TO 230 

y(K)=X(K*l)-X(K-l)-A(K)»Y(K*l)*A(K-l)“Y(K-l) 

X(K)=A(K)*x*k-1 ) -A*k- 1 )*X (K*l )*A(K-1)«A(K)*<Y(K*1)-Y(K-1)) 
X(K)=X(K)/(A(K)-AK-1) ) 

GO TO 240 

ELIMINATE EXTRANEOJS BREAKPOINTS 
X(K)=XREF (K-1) 

Y(K)=YREF (K-I) 

SAR(K)=SAR<K-I) 

K ^ K 1 

SWEEP the breakpoints ON THE VARIABLE SWEEP PANEL 

(IT ALSO KEEPS sweep ANGLES IN FIRST OR FOURTH QUADRANTS) 

K=K*1 

SAR(K-1)=SAI*RSAR(K-1) 

IF (SAR(K-1 ) .LE.PIT) go TO 270 
SAH (K-1 )=SAR (K-1 ) -3. 1415927 
60 TO 260 

IF (SAR(K-l) .GE. l-PIT) ) GO TO 280 
SAR(K-1)=SAR(K-1) *3.1415927 
60 TO 270 

IF ( (SAR (K-1 ) ) .LT. .0 ) GO TO 290 

IF (SAR (K-1 ) -P IT ) 320*300*300 

IF (SAR(K-1 ) *PIT) 310*310*320 

A(K-1>=A/Y 

GO TO 330 

A(K-1 )=-AZY 

GO TO 330 

A(K-1)=TAN(SAR(K-1) ) 

KK=MCD(K*IT) 

Y°Kl = YSnT) *1 yREF (K) -Y5 ( IT ) ) *C0S ( S A I ) - ( XREF (K)-XS( IT))*SIN(SAI) 

X (K ) =XS ( IT ) ♦ ( XREF ( K ) -XS (IT ) ) *C0S (SAI ) ♦ ( YREF ( K ) -YS ( I T ) ) *S IN ( SAI ) 
60 TO 250 

determine the trailing edge intersection 

between FIXED AND VARIABLE SWEEP WING SECTIONS 
IF (ABS(»SAR(K)-5AR(K-1) ) .LT. {. 1/RAO) ) GO TO 360 
Y(K)=XREF (K*1)-X(K-1)-A(K)*YREF(<*1)*A(K-1)»Y(K-1) 

Y (K) =Y(K) / ( A (K-1 ) -A (K) ) tv ivi 

X(K)=A(K)*X(K-1)-A(K-1)*XREF(K*1)*A(K-1)*A(K)*(YREF(K*1)-Y(K-1)> 

X(K)=X(K)/(A(K)-A(K-1)) 

GO TO 370 
X(K)=XREF(K*1) 

Y(K)=YREF (K*1 ) 

K ^ 1 

STORE reference planform cooroinates on inboard fixed trailing 
EDGE 

DO 380 I=K*N1 
X(I)=XREF (I) 

Y( I)=YREF ( I ) 

SAR(I-1)=RSAR(I-1) 

DO 400 I=1*N 


GEO108O 

GE01890 

GE01900 

GEO1910 

GEO1920 

GE01930 

GE01940 

GE01950 

GE01960 

GE01970 

GE01980 

GE01990 

GE02000 

GEO2010 

GEO2020 

GE02030 

GE02040 

GE02050 

GE02060 

GE02070 

GE02080 

GE02090 

GE02100 

6E021 10 

6E02120 

GE02130 

GEO2140 

GE02150 

GE02160 

GE02170 

GE02180 

6E02190 

GE02200 

GE02210 

GEO2220 

GE02230 

GE02240 

GEO2250 

GE02260 

GEO2270 

GE02280 

GEO2290 

GEO2300 

GEO2310 

GEO2320 

GE02330 

GE02340 

GE02350 

6E02360 

GE02370 

GEO2380 

GE02390 

GE02400 
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XX{I.IT)=X<I) 

yv(i.iT)=y(i) 

HMCD(I*II )=MC0(I*IT) 

TTi«(0(I»ir)=0IH(I*IT) 

<*00 A5(I»IT)=A(I) 

XX (N1 . IT) =X (,N1 ) 
yy(Ni.iT)=Y(Ni) 

AN( IT) =AAN (IT) 

<*10 COMTINUE 
C 

c line up breakpoints among planeorms 

c 

BOTSV ( 1 ) =B0TSV (2) =0. 

WRITE (6*9B0) 

DO 530 IT=1,iplAN 

NIT=AN(IT)*1 

DO <*70 ITT = 1,IPLAN 

IF (ITT.EO.IT) go to <*70 

NITT=AN(ITT) *1 

DO <*60 I = 1,NITT 

JPSV=0 

DO <*20 JP=1,NIT 

IF ( yy ( jpf IT ) .Eo.yy ( I * ITT ) ) go to <*60 

<*20 CONTINUE 

DO <*30 JP=1,nIT 

IF <yy ( jp* IT ) .LT.yr ( I , ITT ) ) go to <*<*o 
<♦30 CONTINUE 
GO TO <*60 
<*<♦0 JPSV=JP 

IND=NIT- ( JPSV-l ) 

DO <*50 JP=1*IN0 

K2=NIT-JP»2 

K1=NIT-JP*1 

XX (K2* IT) =XX (K1 .IT) 

yy(K2»iT)=yy(Ki,iT) 

MHCD (K2. IT ) =HMCD (K 1 . IT ) 

AS(K2,IT)=AS(K1,IT) 

<*50 TT«(DtK2,IT)=TTW0(Ki,iT) 
yy(JPSV,IT)=YY(I,ITT) 

AS( JPSV. IT) = AS ( JPS(/-1 , IT) 

TTWD ( JPSV. I T ) =TTWD ( JPSV-1 . I T ) 

^ XX ( JPSV. IT ) = (yy(JPSV.IT)-YV( JPSV-I . it ) ) *AS ( JPSV- l * 

MMCD ( JPSV. IT ) =MMCD ( JPSV- 1 . I T ) 

AN(IT)=AN(IT).l. 

Nlf=NIT*I 
<*60 CONTINUE 
<*70 CONTINUE 
C 

C SEQUENCE WING COORDINATES FROM TIP TO ROOT 

N1=AN( IT) ♦! . 

DO <*80 1 = 1. N1 
<*80 Q(I)=YY(I.IT) 

DO 520 J=1.N1 
HIGH=1 . 

DO A90 I=1.N1 

IF ( (Q(D-HIGH) .GE.O.) GO TO <*90 

HIGH=0(I) 

m=i 


GEO2<*10 

GE02<*20 

GE02<*30 

GE02490 

GE02<*50 

GE02<*60 

GE02<*70 

GE02<*80 

GE02490 

GE02500 

GE02510 

GEO2520 

GE02530 

GE025<*0 

GEO25S0 

GE02560 

GE02570 

GE02580 

GE02590 

GE02600 

GEO26I0 

GE02620 

GE026JO 

GE026<*0 

GE02650 

GE02660 

GE02670 

GE02680 

GE02690 

GE02700 

GEO27I0 

GE02720 

GE02730 

GE02740 

GE02750 

GE02760 

GE02770 

GEO2780 

GE02790 

GE02800 

GE02810 

T)*XX(JPSV-1.ITgE02820 

GE02830 

GE02840 

GE02850 

GE02860 

GE02870 

GE02880 

GE02890 

GE02900 

GE02910 

GE02920 

GE02930 

GE029<*0 

GE02950 

GE02960 

GE02970 

GE02980 

GE02990 

GE03000 
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490 


500 


510 

5^0 

530 

C 

c 

c 

c 

c 


c 

c 

c 


540 

550 


DbO 

570 


560 

590 


600 

610 

C 

C 

c 

c 


CONTINUE 

IF (J.NIE.l) GO TO 500 

80TSV<IT)=HIGH 

KFX ( IT)=IH 

0(IH)=1. 

5PY ( Jf IT)=HIGH 

IF ( IH.Gl .KFX ( IT ) ) GO TO 510 

IYL(J.IT)=1 

IYT(JtIT)=0 

GO TO 520 

IYL(J»IT)=0 

IYT(J.IT)=1 

CONTINUE 

CONTINUE 

SELECT MAXIMUM B/2 AS THE WING SEMISPAN. IF bOTH FIRST AnD 
SECOND PLANFORMS HAVE SAME SEMISPAN THEN THE SECOND PLANFOrM IS 
taken TO BE THE KING. 

KB0T=1 

IF (BOTSV (1 ) .GE.BOTSV (2) ) KB0T=2 
BOT=BOTSV(KBOT) 

COMPUTE NOMINAL HORSESHOE VORTEX WIDTH ALONG WING SURFACE 
TSPAN=0 

iSAVE=KFX(KBOT)-l 
I=KFX (KBOT ) -2 
IF (I.EO.O) GO TO 550 

IF (TTWD<I<KBOT).Ea.TTWO(ISAVE.KaoT)) GO TO 560 
CTWD=C0S«ATAN( TTWDI ISAVE.KBOT ) ) ) 

TLGTH= ( YY ( 1SAVF*1 *KBOT ) -YY ( I ♦ 1 ,KBOT ) ) /CTWO 

tspan=tspan*tlgth 
IF (I.EO.O) GO TO 570 
ISAVE=I 
1 = 1-1 
GO TO 5A0 
VI=TSPAN/VIC 
VST0L=VI/2 
TSPANA=0. 

KBIT=2 

IF (IPLAN.EQ.l) GO TO 610 
IF (KbOT.EO.2) KBIT=1 
ISAVEA=KFX (KBIT) -1 
IA=KFX(KctlT)-2 
IF (lA.EO.O) GO TO 590 

IF (TTWO( I A.KBIT ) .EO.TTWO ( ISAVEA.KBIT) ) GO TO 600 
CTwOA=COS(ATaN(TTwO( ISAVEA.KBIT) ) ) 

TLGTHA=(YY( ISAVEA*! .KBIT)-YY(IA»1 .KBIT) )/CTWDA 

TSPANA=TSPANA.TLGTHA 

IF (lA.EO.O) GO TO 610 

ISAVEA=IA 

IA=1A-1 

GO TO 580 
CONTINUE 

ELIMINATE PLANFORM BREAKPOINTS WHICH ARE WITHIN (B/2)/2000 UNITS 

laterallt 

DO 630 ir=l.lPLAN 


GE03010 

GE03020 

GE03030 

GE030AO 

GE03050 

GEO3060 

GE03070 

GE03080 

GE03090 

GE03100 

GE031 10 

GE03120 

GEO3130 

GEO31A0 

GEO3150 

GE03160 

GEO3170 

GE03180 

GE03190 

GE03200 

GEO3210 

GE03220 

GEO3230 

GEO32A0 

GE03250 

GE03260 

GE03270 

GE03280 

GEO3290 

GE03300 

GEO3310 

GEO3320 

GE03330 

GE033AO 

GE03350 

GE03360 

GE03370 

GE03380 

GE03390 

GE03AOO 

GE03A10 

GE03A20 

GEO3A30 

GE034AO 

GE03450 

GE03460 

GE03470 

GE03480 

GE03490 

GEO3500 

GE03510 

GE03520 

GE03530 

GE035A0 

GEO35S0 

GEO3560 

GE03570 

GE03580 

GEO3590 

GE03600 
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b20 

630 

C 

c 

c 


640 

6b0 


660 


670 

C 

C 

C 


6b0 

690 

C 

c 

c 


c 

700 


N=AN(IT) 

M=N*1 

00 630 J=1,N 

AA=ABS(SPY(J.IT)-S3Y(J*1,IT) ) 

IF (AA.EQ,O..Oi5.AA.GT.ABS(TSPAN/aOOO.) ) GO TO 630 
IF (AA.GT.YT0L) WRITE (6*1010) SPY ( J* 1 . I T ) , SP Y IJ. I T ) 
00 620 1=1 *N1 

IF (YY(I.IT) .NE.SPY(J*1,IT) ) GO TO 620 
YY(I»IT)=SPY( 3, IT) 

CONTINUE 

SPY(J*1.IT)=SPY(J*IT) 

CONTINUE 


COMPUTE 2 COORDINATES 


00 670 IT=1,IPlAN 
JM=NI=ANHT)*1. 

00 640 J2=I*n1 
2Z(JZ.IT)=RTC0HT(IT) 

JZ = I 
JZ=JZ*I 

IF (JZ.GT.KFX(IT) ) GO TO 660 

ZZ ( JZ* IT) =ZZ ( JZ-1 * IT) ♦ ( YY ( JZ* IT) -YY ( JZ-1 . IT) ) *TTWO ( JZ- 1 ♦ I T ) 

GO TO 650 

JM=JM-1 

IF (JM.EQ.KFX ( IT) ) GO TO 670 

ZZ( JM* IT)=ZZ UM*I. IT) ♦ (YY ( JM. IT)-YY( JM*l* IT) ) *TTWD ( JM. I T ) 

GO TO 660 
CONTINUE 


WRITE PLANFORM PERIMETER POINTS ACTUALLY USED IN THE COMPUTATIONS 


WHITE (6*900) 

DO 690 n = l , IPlAN 
N=AN(IT) 

NI=N*I 

IF (IT.EQ.2) write (6*1000) 
DO 680 KK = 1 ,N 


T0UT=ATAN(TTWD(KK. IT) )*RAD 
A0UT=ATAN(AS(KK*IT) )oRAO 
IF (AS1KK.IT) .EQ.AZY) AOUT=90. 
WRITE (6*910) KK*XX (KK* IT) *YY (kK 
IT) 

CONTINUE 

WRITE (6*910) N1 *XX (Ml * IT) *YY (Nl 
CONTINUE 


IT) *ZZ(KK. IT) *AOUT*TOUT*MmCO(KK, 
IT) *ZZ(N1,IT) 


PART ONE - SECTION THREE - LAY OJT YAWED HORSESHOE VORTICES 
STRUE=0. 

NSSWSV ( I ) =NSSWSV (2) =MSV ( I ) =MSV (2) =0 
DO 780 II=1*IPLAN 
Nl=AN(IT) *1 . 

1 = 0 
J=I 

YIN=B0TSV(IT) 

ILE*ITE=KFX(it) 

determine spanwise Borders of horseshoe vortices 

iXL=IXT=0 
i = i*i 


GE03610 
GEO3620 
GE03630 
GE03640 
GEO3650 
GE03660 
GE036 7 0 
GE03680 
GE03690 
6E03700 
GE03710 
GE03720 
GE03730 
GE03740 
GE03750 
GE03760 
GE03770 
GEO3780 
GE03790 
GEO3800 
GEO3810 
GE03820 
GEO3830 
GE03840 
GE03850 
GE03860 
GE03870 
GEO3880 
GE03890 
GE03900 
GEOJrIO 
GEO3920 
GEO3930 
GE03940 
GE03950 
GE039t>0 
GE03970 
GEO3980 
GE03990 
GE04000 
GE040 1 0 
IGEO4020 
GE04030 
GEO404O 
GE04050 
GE04060 
GE04070 
GE04080 
GE04090 
GE04100 
GE041 10 
GE04120 
GE041 30 
GE041<*0 
GEO4150 
GE04160 
GE041 70 
GE04180 
GE04190 
GE04200 
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10 


c 

7dQ 


730 


740 


750 


C 

760 


C 

c 

c 

c 


770 

c 

c 

c 


OF BWtAKpOlNl 
BORDER 


IF (YIN.0E.(SPY(J»IT)^VST0D) go to 710 
bOPOER It> within vortex SPACING TOLERANCE (VSTOL) 
therefore use the NiEXT BREAKPOINT INBOARD FOR THE 
VBORD(I)=YIN ^ 

GO TO 740 

USE NOMINAL VORTEX SPACING TO DETERMINE THE BORDER 
VBORD( I ) =SPY ( IT) 

COMPUTE SUBSCRIPTS ILE AND ITE TO INDICATE WHICH 

BREAKPOINTS ARE ADJACENT AND WHETHER THEY 4RE ON THE WING LEADING 

EDGE OR THE TRAILING EDGE 
IF (J.GE*N1) GO TO 730 

IF (5PY(Jf IT) •NE.SPY(J^lf IT) ) GO TO 730 
ixl=ixloyl ( Jf IT) 

IXT=IXT^1YT(J*IT) 
j= j* 1 
00 TO 720 
YIn=SPY( J t IT) 

1XL=IXL*IYL ( J. IT) 

1XT=IXT*1YT ( J» IT) 

J=J*1 

CPHI=COStATAN(TTWD( I lE. IT) ) ) 

IPHI=ILE-IXL 
IF (J.GE.Nl) IPHI=1 

YIN = YIN-VI*C0S ( ATAV) ( TTWO ( IPHI ♦ IT) ) ) 

IF (I.NIE.l) GO TO 760 
ILE=ILE-IXL 
ITE=ITE»IXT 
GO TO 700 

COMPUTE COORDI^JATES FOR CHORDWISE RO«f OF HORSESHOE VORTICES 
YQ= (VBOROI I-l ) ♦VBORO( I ) ) /2. 

HW= ( VB0R0< I ) -VBORD ( I-l ) ) /2. 

IM1=I-1*NSSWSV(1) 

ZH(IM1)=2Z(ILE»IT)*{YQ-YY(ILE.IT))*TTW0(ILE«IT) 

PHI ( IMl ) =TTW0( ILE* IT) 

SSrtWAdMJ )=AS(ILE*IT) 

XeE=XX(1lE*IT)*AS(ILE.IT)*(YQ-YY(ILE.IT)) 

XTE=XX(irE.IT) ♦AS(ITE.IT)*(Ya-YY(ITE.IT)) 

XLOCAL=(XLE-XTE)/T3LSCli) ( IMl) 

COMPUTE <<ING AREA RROJECTEO TO THE X - Y PLANE 

STRUE = STRUE*XL0CAL*T8LSCii( (IM1)»(HW*2.)*2. 

NSCW = TBlSCW ( IMl) 

DO 770 jC«=ltNSCW 
AJCW=JCW-1 

xlel=xle-ajcw*xlocal 
NTS=JCV*MSV(1) *MSV<2) 

PN(NT5)=XLEL-.25*XL0CAL 
PV (NTS) =XL el- . 75*X LOCAL 

PS I (NTS) = ( (XlE-PN(NTS) ) * AS ( I T£ , I T ) ♦ (PN (NTS ) -XTE ) * AS ( ILE ♦ I T ) ) / ( XLE 
IXTE) 

S(nTS) =H«*/CPHI 
Q(NTS)=YU 

continue 

MSV(IT)=MSV(IT)*NSCW 


TEST TO DETERMINE <(HEN WING 
IE (VBORD(I) .LT.YREG(1*IT) ) 


root is REACHED 
GO TO 750 


GEOA210 

GEOA220 

GEOA230 

GEO42A0 

GE0A2S0 

GEOA260 

GE0^270 

GEO42a0 

GE04290 

GE04300 

GE04310 

GE04320 

6E043JO 

GE04340 

GE043SO 

GE04360 

GEO4370 

GE04380 

6E043RO 

GE04400 

GE04410 

GE04420 

GE04430 

GE04440 

GE04450 

GE04^*60 

GE04470 

GE044HO 

GE044VO 

GE04500 

GE04510 

GE04520 

GE04530 

GE04540 

GE04550 

GE04560 

GE04570 

GE045BO 

GE04590 

GE04600 

GE04610 

GE04620 

GEO4630 

GE04640 

GE04650 

GE046bO 

GE04670 

GE04680 

GEO4690 

GE04700 

-GE04710 

GE04720 

GE04730 

GE04740 

GE047SO 

GE04760 

GE04770 

GEO4780 

GEO4790 

GE04800 
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NSSWSVtlT)=I-l 

GE04B10 

780 

continue 

GEO4820 


M=HSV(1)*MSV(2) 

GE04830 

C 


GE04840 

C 

COMPUTE ASPECT RATIO AND AVERAGf CHORD 

GEO405O 

c 


GE04860 


BOT=-BOT 

GE04870 


AR=4.*B01»B0T/SREE 

GEO480O 


ARTRUE=4.*BOT»0OT/STRUE 

GE04890 


CAVE=STRUE/(2.*B0T) 

GE04900 


8ETA= ( 1 ,-MACH*MACH) *«.5 

GE04910 


NVTWO=0 

GE04920 


DO 810 lT=ltIPLAN 

GE04930 


NVONE=l*(IT-n*MSV(l) 

GEO4940 


NVTWO=NVTWO*MSV (IT) 

GE049SO 


IF (TWIST ( IT ) .LE.O. ) GO TO 7R0 

GEO4960 


READ (5.OB0) (ALP(NV) »ALP(NV*1) .ALP(NV*2) ,alP(NV* 3) t 

ALP (NV.4) . ALP (GE04970 


1NV*5) »ALP(NV*6) .AL3(nV*7) . NV=NV0NE *N VT WO. 8 ) 

GEO490O 


GO TO 810 

GE04990 

790 

DO 800 NV=NVONF .NVTWO 

GE05000 

800 

ALP(NV)=0. 

GE05010 

810 

CONTINUE 

GE05020 


WRITE (6.1040) M 

GE05030 


WRITE (6.1050) ( IT,MSV ( IT) ,N5SWSV (IT) . IT=l , IPLAN) 

GE05040 


IF (SCW.NE.O.) WRITE (6fl020) SCw 

GEOSoSO 


IF (SCW.EO.O.) WRITE (6tl030) ( TBLSCW ( I ) . 1=1 ,NSTA) 

GE05060 

C 


GE05070 

c 

apply PRANDTl-GLAUERT CORRECTION 

GEO5O0O 

c 


GE05090 


DO 820 NV=1,m 

GE05100 


PSI (NV)=ATAN(PSI (NV)/BETA) 

GE051 10 


PN(NV)=PN(NV)/BETA 

GE05120 

820 

PV(NV)=PV(NV)/0ETA 

GE05130 


RETURN 

GE05140 

830 

ICOOEOF=l 

GE051SO 


WRITE (6.930) CONFIG 

GE05160 


RETURN 

GEO5170 

840 

IC00E0F=2 

GE05180 


WRITE (6.940) K.IT 

GE05190 


RETURN 

GE05200 

850 

IC00E0F=3 

GE05210 


WRITE (6.970) PTEST.QTEST 

GE05220 


RETURN 

GE05230 

C 


GE05240 

c 


GE05250 

c 


GE05260 

860 

FORMAT (1H1//63X,13HGE0METRY DATA) 

GE05270 

870 

FORMAT (///45X.A10.22HREFERENCE =>LANFORM HAS.I3.7h CURVES// 1 2X , 19H6E052BO 


IROOT Chord height =.F12.5.4X,29HVARIABlE SWEEP PIVOT 

POSIT ION. 4X.6GE05290 


2HX(S) =.F12.5,5X,6HY (S) = . F 1 2 . S//46X .40HBREAK POINTS 

FOR The REFEHGE05300 


3ENCE PLANFORM /) 

GE05310 

880 

FORMAT (6F10.4) 

GE05320 

890 

FORMAT ( 1H1//47X* 1 7HCONFIGURAT ION NO..F8.0/) 

GE05330 

900 

FORMAT (22Xt5HP0INT t6X. 1 HX 1 11 X 1 1 HY . 1 1 X t IHZ . 1 OX .5HSWEEP . 7X . 8H0I HEDRGE05340 


1 AL. 4X»4HM0VE/68X. 5HANGL E.8X.5H angle. 6X,4HC00E/) 

GE05350 

910 

FORMAT (20X. I5.3F12.5.2F14.5. 16) 

GE05360 

920 

FORMAT (/40X,5HCURVE, I3*9H IS SWEPT ,F12.5, 20H DEGREES ON PL ANFORM t GE05370 


113) 

GEO530O 

930 

FORMAT UHI///41X»43hEND OF FILE ENCOUNTERED AFTER CONFIGURATION, FGE05390 


1 /.O) 

GE05400 
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940 FORMAT < iHl /// I 8X . 45 hThE FIRST VARIABLE SWEEP CURVE SPECIFIED (K =GE05410 
1*I3.44h ) DOES not HAVE AN M CODE OF 2 FOR PLANF0RM»I4) GE05420 

9S0 FORMAT (8Fb.I*F10.A,F5.1.F10.4,F5.I) GEO5430 

960 FORMAT ( 26X , 1 5 , 2F 1 2 .5 ♦ 2F 1 6 . 5 *4X * 14) GEOS440 

970 FORMAT ( 1 HI //30X , 36HERR0R - PROGRAM CANNOT PROCESS PTEST =,F5. 1 1 12GEOS450 
IH AND QTEST =,FS.l) GE05460 

9B0 FORMAT ( // 48 X . 3SHBRE AK POINTS FOR THIS CONFIGURATION//) GE05470 

990 FORMAT ( 28X » 5HP0 1 NT ♦ 6X * IHX* IIX* IHY. I IX » 5HSWEEP » 1 OX t BHOIHED rAL t 7X ,4GE05480 
1HM0VE/38X,3HREF»9X,3HREF» lOX.SHANGLEf UX«5HANGLEt9X.4HC00L/) GE05490 

1000 FORMAT ( /52X , 28HSEC0ND PLANFOHM BREAK POINTS/) GE05SOO 

1010 FORMAT (////25X.3AHTHE BREAKPOINT LOCATED SPANWISE AT .F U . 5 * 3* * 20HGE0551 0 
IHAS BEEN ADJUSTED T0.F9.5////) GE05520 

1020 FORMAT ( /43X , F5 . 0 » A 1 H HORSESHOE VORTICES In EACH CHOROWISE ROW) GEO5530 
1030 FORMAT ( /23X , 98HT ABLE OF HORSESHOE VORTICES IN EACH CHORDWISE ROW GEOSS^O 
KFROM TIP TO ROOT BEGINNING WITH FIRST PL ANFORM ) //25F5. 0/2SF5. 0 ) GEOSSbO 
1040 FORMAT (///33X I5»62H HORSESHOE VORTICES USED ON THE LEFT HALF OF ToEOSSbO 
IHE CONFIGURAT10N//50X»36HPLANFORM TOTAL SPANWjSE/) GE05570 

1050 FORMAT ( 52X . 14 , 1 OX , 1 3 , 1 1 X . 14 ) GEOS580 

END GE05S90- 
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OVERLAY (<iINGTL»2*0) 

mat 

10 


program matxsol 

MAT 

20 


DIMENSION YY(2)* Ev(2)* FVI(2). EVN(200) 

MAT 

10 


COMMON /All/ BOT.M*BETA,PTEST*OTEST.TBlSCw(50) »Q(;?00) »PN<200) 

♦PV (2mAT 

40 


100) ♦ ALP (200) ,S (200) »PSI (200) .PHl (50) *ZH(50) 

MAT 

SO 


COMMON /TOTHrE/ CIR(200f2) 

mat 

60 


common /INSUB23/ APS I ♦ APHI . XX . YYY , 2Z f SNN. TOLC 

mAt 

70 

c 


MAT 

80 

c 


mat 

90 

c 

PART 2 - COMPUTE CIRCULATION TERMS 

MAT 

100 

c 


MAT 

110 

c 


MAT 

120 

c 

THE tolerance set AT THIS POINT IN THE PROGRAM MAY NEED TO BE 

mat 

no 

c 

CHANGED FOR COMPUTERS OTHER THAN THE COC 6000 SERIES 

MAT 

140 

c 


MAT 

ISO 

c 


mat 

160 


TOLC=(BoT*15.E-05) **2 

MAT 

170 


DO 10 Nv=l*M 

mat 

no 


CIR(NV. 1 ) =12. 5663704* ALP (NV) 

mat 

190 


CIR (NV»2) =12.5663704 

MAT 

?00 


IF (PTEST.NE.O. ) CIR(NV*2)=-1.0R64l55*a(Nv)/BOT 

mAt 

210 


IF (QTEST.NE.O. ) CIR(NV,2)=-1.0B64155*oV(NV)*BETA 

mat 

220 

10 

CONTINUE 

mAt 

230 


IZZ = 1 

mat 

240 


NNV=TBLSCW ( IZZ) 

mat 

250 


REWIND lO 

mat 

260 


DO 70 NV=1»M 

mat 

270 


DO 20 I=1.M 

mat 

280 

20 

FVN(I)=0. 

mAT 

290 


IZ = 1 

mat 

100 


NNN=TBLSCW(IZ) 

mat 

310 


Do 60 NN=1»M 

MAT 

320 


APHI=ATAN(PHI ( IZ) ) 

MAT 

130 


APSI=PSI (NN) 

mat 

340 


XX=PV (NV) -PN (NN) 

MAT 

ISO 


YY(1)=Q(NV)-0(NN) 

mat 

360 


YY(2)=0(NV) *0(NN) 

MAT 

370 


ZZ=ZH(I2Z)-ZH(IZ) 

MAT 

380 


SNN=S(NN) 

mat 

390 


DO 30 1=1,2 

MAT 

400 


YYY=YY ( I ) 

mat 

410 


CALL INFSUB ( BOT , F V ( I ) , FW ( I ) ) 

MAT 

420 


APH1=-APHI 

mAt 

430 


APS1=-AP5I 

MAT 

440 

30 

CONTINUE 

MAT 

450 


IF (PTEST.NE.O.) GO TO 40 

MAT 

460 


FVN(NN) =Fw ( 1 ) *FW(2) -(FV ( 1 ) *FV(2) )*PHl ( IZZ) 

mat 

470 


GO TO 50 

MAT 

480 

i*0 

FVN(NN)=FW(1)-FW(2)-(FV(1)-FV(2) )*PHI(IZZ) 

mat 

490 

50 

IF (NN.LT.NNN.OR.NN.EQ.M) go to 60 

mat 

SOO 


IZ=IZ*1 

mat 

sio 


NNN=NNN*T8LSCW ( IZ) 

mat 

520 

60 

CONTINUE 

mat 

S30 


DUM8=-CIR (NV, 1 ) 

mat 

S40 


DUMY=-CI« (NV,2) 

mat 

S50 


write (10) (FVN( I ) , 1=1 ,M) ,DUMB,0UMY 

mat 

S60 


IF (NV.lT.NNV.OR.NV.EQ.M) go to 70 

mat 

570 


IZZ=IZZ*1 

mat 

500 


NNV=NNV*T8LSCW( IZZ) 

MAT 

S90 

70 

CONTINUE 

mat 

600 


CALL SSLESO (M,2) 

mat 

610 


RETURN 

mat 

620 


END 

MAT 

630- 
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SUBfiOUTINE SSLESO (NTtNCFLG) 

COMMON /TOTHRE/ CH<200,2) 

DIMENSION RV(205)* CV(205)« R(205)« V(10350) 

REWIND 10 

N1=NT*NCFLG 

J=N1-1 

READ <10) (R(I)»I=l*Nl) 

DO 10 I=1.J 

10 V ( I)=-R (1*1 )/R ( 1 ) 

IN = 1 

20 READ (10) (R ( I ) * 1 = 1 *N1 ) 

12 = 0 

DO ^0 I=1»J 
RV(I)=0. 

DO 30 11 = 1 » IN 
12=12*1 

30 RV < I ) =RV ( I ) *R ( I I ) *'/ ( 12) 

N2=I*IN 

40 RV ( I ) =RV ( I ) *R <N2) 

12=IN*1 

NN=J*IN*1 

KK=J*I2 

J=J-1 

DO 60 1=1 »J 
DO 50 II=1*IN 
NN=NN-1 
KK=KK-1 

50 V(KK)=V(NN) 

e>0 KK=KK-1 

DO 70 I=1*IN 
70 R(I>=V(I> 

K = 0 

DO 90 I=1»J 
CC=-RV ( 1*1 )/RV ( 1 ) 

DO 80 1 1=1 » IN 
CV ( 1 1 ) =CC«R (II) 

NN=K*I I 
12=12*1 

BO V(NN)=CV(1I)*V(I2) 

K=NN*1 
12=12*1 
90 V(K»=CC 

IN=IN*1 

IF (J.EO«NCFLG) GO TO 100 
GO TO 20 
100 K=l 

DO 110 J=1.NCFLG 
DO 110 I=1»NT 
CIR ( I * J) =V (K) 

K=K*1 

110 CONTINUE 

RETURN 
END 


SSt 

10 

SSL 

^0 

SSL 

30 

SSL 

40 

SSL 

50 

SSL 

60 

SSL 

70 

SSL 

80 

SSL 

90 

SSL 

100 

SSL 

no 

SSL 

no 

SSL 

130 

SSL 

\^Q 

SSL 

ISO 

SSL 

160 

SSL 

170 

SSL 

ISO 

SSL 

190 

SSL 

?00 

SSL 

210 

SSL 

220 

SSL 

230 

SSL 

240 

SSL 

250 

SSL 

260 

SSL 

270 

SSL 

280 

SSL 

290 

SSL 

300 

SSL 

310 

SSL 

320 

SSL 

330 

SSL 

340 


SSL 350 
SSL 360 
SSL 370 
SSL 380 
SSL 390 
SSL 400 
SSL 410 
SSL 420 
SSL 430 
SSL 440 
SSL 450 
SSL 460 
SSL 470 
SSL 480 
SSL 490 . 
SSL 500 
SSL 510 
SSL 520 
SSL 530- 


43 


APPENDIX C 


10 


20 


30 


<*0 

c 

c 

C 


OVERLAY(wI!^GTL.3»0) 

program aerooyn 


DIMENSION) YCP(2). CLCC(200f2). 
1. CP<200). P(200)* SMOAD(2,50) 

COMMON /All/ got.m»beta,ptest* 

100). ALP (200) ,5(200) .PSI (200),P 
COMMON /TOTHPE/ CIR(200,2) 
COMMON /THREFOR/ CCAv (2,50) ,CL 
COMMON /ONETHRF/ T ,( I ST ( 2 ) , CREF 
1DHT(2) ,C0nFIG,nSSW5V (2) ,MSV(2) 
2) ,XT (2) .CLWB.CMCL'CLA (2) .BLAIR 
COMMON /THRECOI/ 5LOAD(3,50) 
Common /InsuB23/ a=si , aphi ,xx. 


CH(2,50). SIIM(2), AC(2), 
5LDT(50), SMLO(2,50) 
OTEST,TBlSCw( 50) ,0(200) ,P) 
Hi (SO) ,ZH(50) 

T.CLNT.NSSW, ALPD 
.SREF.CAvE.CLOES.STRUE.AR, 
,KBOT,PLAN, IPLAN.MACH.SSW) 
(50) ,CLAMAR(2) 

YYY./Z.SNN.TOLCSO 


PART 3 - COMPUTE OUTPUT TERMS 


RAD=57. 29578 
TWST=TWIST ( 1) ,TWIST (?) 
ALREf=l 
QINF=1 . 

NSSW=NSS»(SV(1 ) ,NSSi>(SV(2) 


PART 3 - SECTION 1 

COMPUTE LIFT AND PITCHING MOMENT HERE 


10 


1Z = 1 

NNN=TBLSCW ( IZ ) 

DO 10 I=1,M 

P(I)=S(I)»C0S(ATAN(PHKIZ) ) ) 

IE ( I ,LT .NNN.OP. I . EO.M) GO TO 
IZ=IZ»1 

NNN = NNN*TBLSCW ( IZ) 

CONTINUE 
DO 20 NV=1,2 
SUM(NV)=0 
DO 20 I=1,M 

SUM (NV ) =SUM (NV) *C IR ( I ,NV) *P ( I ) 

IF (NV.EG. 1 . AND. I .EO.MSV ( 1 ) ) 

IF (NV.EQ.2.AND.I.E0.MSV(1) ) 

CONTINUE 
CLT=8.*SUM ( 1 ) /SREF 
CLNT=8,*SuM(2)/SREF 
IF (KBOT.EO.I) GO TO 30 
CLi«(NGT = CLT-CLWNGT 

clwing=clnt-clwing 

CRL=0. 

DO AO I=1,M 

CRL=CRL* (0(I)«CIR(I,2)*2.*P(I))*2. 
CLCC (I,1)=CIR(I,1) *2*P ( I ) / (CAVE*S ( I ) ) 
CLCC (I,2)=CIR(T,2) »2*P ( i ) / (CAVERS ( I ) ) 

COMPUTE CLP 

CLP=CRL/ (SREF*90T*0. 08725) 

CLA(2)=CLNT 
DO 120 IXX=1,2 


CLWNOT=SUM ( 1 ) 08 ./SREF 
ClWING=SUM(2)*8./SREF 


AER 

10 

AER 

20 

AER 

30 

AER 

40 

1 ) AER 

50 

AER 

60 

2AER 

70 

AER 

80 

AER 

90 

AER 

100 

CAEP 

110 

2AER 

120 

AER 

130 

AER 

140 

AER 

150 

AER 

160 

AER 

170 

AER 

160 

AER 

190 

AER 

?00 

AER 

210 

AER 

220 

AER 

230 

AER 

240 

AER 

250 

AER 

260 

AER 

270 

AER 

280 

AER 

290 

AER 

300 

AER 

310 

AER 

320 

AFR 

330 

AER 

340 

AER 

350 

AER 

360 

AER 

370 

AER 

380 

AER 

390 

AER 

400 

AER 

410 

AER 

420 

AER 

430 

AER 

440 

AER 

450 

AER 

460 

AER 

470 

AER 

480 

AER 

490 

AER 

500 

AER 

510 

AER 

520 

AER 

530 

AER 

540 

AER 

550 

AER 

560 

AER 

570 

AER 

580 
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30 


60 

70 


dO 


'iO 

100 


110 

120 


c 

c 

c 

c 


GO TO 60 


6A=S0=SC=O. 

1=0 

J8=NSSWSV(1) 

JA=1 

CONTINUE 

DO 70 JSi>W = JA.JB 
SD=SE=0. 

SLOADI IXX. JSSW) =0 
NSCW = TBtiC*( ( JS5W) 

DO 70 JSC«=1.NSCW 
IE (TkST.EO.O. . AND. IXX.EO. 1 ) 

I = I*1 

SA=SA*CH 1 1 . IXX) *P ( I ) 
b8=SB*CH ( I ♦ IXX ) ( I ) *P ( I ) 

SC = SC*CI><( I . IXX) *PN( I ) *P ( I)“BETA 
SLOAD ( IXX, JSSW ) =SLOAO ( IXX * jSSiK ) ♦ (BOT*CIH ( I . D 
1 (IXX) ) 

SD=SD*CIf< ( I ♦ IXX) 

SE = SE*CIP ( I I IXX ) *PN ( I ) "BETA 
IF ( JSCW.NE.NSCW) GO TO 70 
SMOADI IXX, JSSW) =SE 
SMLD( IXX* JSSW) =SD 
GO TO 70 

SLOAD ( 1 , JSSW ) =SMOAD (1 , JSSW ) =SMLO ( 1 , JSSW) =0 . 

CONTINUE 

IF ( JSSW.GE.NSSW) GO TO 80 

JA=NSSWS'/(1) .1 

JB=:NSSW 

IF (IXX.EQ.l) GO TO 50 

SC2=SC 

SA2=SA 

CL AM AW ( 1 ) =SC/ ( SA"CPEF ) 

GO TO 50 
CONTINUE 

IF (IXX.tO.l) GO TO 100 
IF (IPLAN.EO.l) GO To 00 
SC3=SC-SC2 
SA3=SA-SA2 

CLAMAH (2) =SC3/ <SA3»CREF ) 

GO TO 100 

CLAMAH(D=SC/(SA*CREF) 

CONTINUE 

IF (TWST.EO.O. .AND. IXX.EQ.l ) 
YCP(IXX)=SB/(SA*BOT) 

AC ( IXX) =SC/ (SA"CREF) 

GO TO 120 

YCP(1)=AC(1)=0. 

CONTINUE 
CMCL=AC(2) 

CMO=(AC(l)-AC(2) )"CLT 

part 3 - SECTION 2 
compute other- and PRINT ALL FINAL 

DO lAO IXX=1,2 
JN = 0 

DO lAO JSSW=1,NSSW 
CH(IXX,JSSW)=0 
NSCW=TBLi>CW { JSSW) 

00 130 JSCW=1,NSCW 


GO To 110 


- OUTPUT data here 


AER 

610 

AER 

620 

AER 

630 

AER 

6^0 

AER 

650 

AER 

660 

AER 

670 

AER 

6ao 

AER 

690 

AER 

700 

AER 

710 

aER 

720 

AER 

730 

AER 

7^0 

AER 

7S0 

1AER 

760 

AER 

770 

AER 

780 

AER 

790 

AER 

000 

AER 

010 

AER 

020 

AER 

030 

AER 

040 

AER 

050 

AER 

060 

AER 

870 

AER 

080 

AER 

090 

AER 

900 

AER 

910 

AER 

920 

AER 

930 

AER 

940 

AER 

950 


AER 960 
AER 970 
AER 980 
AER 990 
AERIOOO 
AERIOIO 
AER1020 
aER 1 030 
aerioho 
aerioso 

AERIOoO 

AER1070 

aerioso 

AER1090 
AERl 100 
AERlllO 
AER1120 
AER1130 
AERl 190 
AERl 150 
AERI 160 
AER1170 
AERl 180 
AERl 190 
AER1200 
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130 


l<*o 


IbO 


160 


170 

c 

c 


JN=JN*1 

CH( IXX. JS5W) = (-2,0) * (PV ( JN) -PN( JN) )*BE:TA*CH( IXX, JSSW) 
CONTINUE 

CCAV ( IXX* JSSW) =CH( IXX* JSSW) /CAVE 
CLCL(IXX»JSSW)=SL0AD(IXX, JSSK)/CCAV(1XX* JSSw) 

CONTINUE 

CLD=CLDES 

IE (CLDE5.EQ.11) CLO=l. 

DO 150 1 = 1. M 

CP(I)=(CLCC(I.1)*CLCC(I.2)*(CL0-CLT)/CLNT)*CAVE/(2.*(PN(I)- 

IBETA) 

CONTINUE 

WRITE (6.2A0) CONFIG 
IF (PTEST .NE.O. ) WRITE (6.350) 

IF (QTeST.NE.O, ) WRITE (6.330) 

IF (PTEST. EO.O.. AND. qTEST.EQ.O.) WRITE (6,3aO) 

WRITE (6*360) CLD 
HEA0=8hDESIRED 
IF (CL0E5.EQ.11.) HEAD=8H 
IEND=11 

IF (CLDE5.NE.11.) IEND=1 
DO 190 IUTK=1.IEND 

IF (IENO.EQ.il) CLDES=(FLOAT(IJT<)-1.)/10. 

IF (CLDES.EO.O. ) CLOES=-.l 
NR = 0 

DO 160 NV=1.nS5W 
NSCW=T8LSCW(NV) 

NP=NR*1 

NR=NR*N5CW 

PH I PR= AT an (PHI (NV) ) *RAD 
SLOADO.NV) =0. 

IF (NV.EQ. (NSSwSV ( 1 ) ♦! ) . AND. lENO.EQ. 1 ) WRITE ( 6 . 230 ) 

DO 160 I=NP.NR 
IF (IUTK.GT.D GO TO 160 
PNPR=PN( 1 ) *8eTA 
PVPR=PV ( I ) *8ETA 

PSIPH=ATAN(BETA*TAN(PSI ( I ) ) )*RAD 

WRITE (6.370) PNPR.PVPR.Q ( I ) . ZH (NV) .S ( I ) .PSIPH.PHIPR. ALP ( I ) 
SLOADO.NV) =SLOAO( 3. NV) *CLCC ( I .2 ) *CLOES/CLNT *CLCC ( 1 . 1 )-CLCc 
ILT/CLNT 

IF (IUTK.GT.D GO TO 170 
WRITE (6.270) 

WRITE (6*280) CREF.CAVE*STRUE.5REF.BOT,AR.aRThue*MACH 
CONTINUE 


IF (PTEST. NE. 0. ) WRITE (6*380) CLP 
IF (PTEST. NE.O.) GO TO 220 

CCMPUTE CMQ.CLQ 

CNQ=2.0*CMCL*CLNT/ (0.08725«CREF) 

CLO=2.0*CLNT/ (0.08725*CREF) 

IF (QTEST.NE.O.) WRITE (6.390) CHQ.CLQ 
IF (OTESi.NE.O.) GO TO 220 

COMPUTE INDUCED DRAG FOR FLAT WING-BOOy WITH NO DIHEDRAL 


AER1210 
AER1220 
AER1230 
AER12A0 
AER1250 
AER1260 
AER1270 
AER1280 
AER1290 
•PV ( I ) )«AER1300 
AER1310 
AER1320 
AER1330 
AER1340 
AER1350 
AER1360 
AER1370 
AERl 380 
AER1390 
AER1400 
AERUIO 
AERU20 
AERU30 
AERlA^O 
AERU50 
AERUGO 
AERU70 
AERUaO 
AERU90 
AER1500 
AER1510 
AER1S20 
AER1530 
AERlS-tO 
AER15S0 
AERIBGO 
AER1S70 
.CPU) AER1580 
(I.Sl^CAERlSRO 
AER1600 
AER1610 
AER1620 
AER1630 
AER1640 
AER1650 
AER1660 
AER1670 
AER1680 
AER1690 
AER1700 
AER1710 
AER1720 
AER1730 
AERl 7A0 
AER1750 
AER1760 
AER1770 


NSV=NSSW5V ( 1 ) *1 
MTOT=MSV(l)*l 


AER1780 
AERl 790 
AERiflOO 
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180 


190 


200 


210 

220 

C 

C 

c 

230 

240 

250 


260 

270 


280 

290 


300 


310 


IF (KBOT.tQ.l) GO TO 180 
NSV=NSV*^SSHSV (2) 

MT0T=MT0l*MSV(2) 

Call COICLS (AQtAHTRUE.NSSWSVO<80T).MT0T.NSV»C0I»C0lT) 

CLAPD=CLA(2) /S7.29578 

ALP0=-<CLT/CLA(2) >*57.29578 

ALPO=ClOES/CLAPO*AlPO 

ALPW=1 ./CLAPD 

CLrtB=CLWlN6*ALPD/57.29570*CLrtNGT 
CDIWB = CDI/ (CLWB*CL*/B) 

IF (lUTK.EQ.l) WRITE (6*250) HEAO*COIT 

write (6*260) CLDES* ALPD*CLWB*CDI*C0IW3 

WRITE (6*290) CL A ( 2 ) * CL APO * CLT * ALPO * YC^ ( 2 ) . CMCL *CMO 

WRITE (6*300) CLT 

NR=J=0 

00 210 NV=1 *NSSW 
bclcc=badlae=basld=o. 

NSCW=TBLSCW (NV) 

NP=NR*I 
NR=NR*NSCW 
DO 200 I=NP*NR 
ADLAE=CLCC ( I *2) *CLT/CLNT 
BSLO=CLCC(I*l )-ADLAE 
BCLCC=BCLCC*CLCC (1*1) 

baolae=badlae*adlae 
baslo=basld*bsld 
continue 

SLOT (NV) = (SMOAD ( 1 *NV) ♦5M0AD(2*NV)*(CL0ES- 
1LD(2*NV)*(CL0ES-CLT)/CLNT) 

J=J*NSCW ’ 

YQ=0( J)/80T 

IF (NV.E(J. (NSSWSV ( 1 ) ♦ 1 ) ) WRITE (6*310) 


aERIHIO 
AERI820 
AERI830 
AER1840 
AER1s50 
AER1860 
AER1870 
AER1380 
AER1890 
AER1900 
AER1910 
AERI920 
AER1930 
AER1940 
AER1950 
AER1960 
AER1970 
AER1980 
AER1990 
AER2000 
AER2010 
AER2020 
AER2030 
AER2040 
AER2050 
AER2060 
AER20 7 0 

CLT)/CLND/(SNL0(1 *NV) *SMAER2080 

AER2090 
AER2100 
AER21 10 
AER2120 


WRITE (6*320 ) NV*YO*SLOAD (2*NV ) * CLCL (2*NV) ,CCAV (2*NV) * BCLCC * B AOL A£ AER2 1 30 


I *0ASLO*SlOAD(3*NV) *SlDT (NV) 

continue 

RETURN 


FORMAT 

FORMAT 

17.0//) 

FORMAT 


(/12X,45HSEC0ND PLANFORM HORSESHOE VORTEX DESCRIPTIONS/) 
( 1 Hi ///58X* 16H AERODYNAMIC OAT A///54X * 1 7 hCONF I GUR AT ION 

( IHl *18X*22HC0MPLETE CONFIGURAT ION* 31 X * 25HW ING-BOOY 


AER2140 
AER2150 
AER2160 
AER2170 
AER2180 
AER2190 
AER2200 
NO. *FAER2210 
AER2220 
CHARACTAER2230 


1ERISTICS/64X,4HLIFT*9X*33HINDUCED DRAG (FAR FIELD SOLUT I ON ) // 1 6X A8 AER2240 


2*21h cl computed ALPHA * 1 9X * 6HCL ( W8 >* 7X * 1 3HCDI 
301/ (CL (Wt))»*2) /80X* 12H( 1/ (PI*AR) =F3.5*2H )) 

FORMAT (11X*2F15.5*15X*3F15.5) 

FORMAT (////4X.11H REF. CH0RD*6X*25HC AVERAGE 
14HREFERENCE ARE A * 9X * 3HB/2 * 0X * 7h«EF . AR*8X*7HTRUE 
2BER/) 

FORMAT (BF15.5) 

FORMAT (///47X.38HC0MPLETE CONFIGURATION 
1 ALPHA*8X*53hCL (TWIST) ALPHA AT CL=0 


AT CL(wB) *aX*15hCaER2250 
AER2260 
AER2270 

TRUE area *2X,1aER2280 
AR*4X* 1 IHmACH NUMAER2290 
AER2300 
AER2310 

CHARACTER1STICS//36X*8HCLAER2320 
Y CP CM/CL CMOAER2330 


2/27X*23HPER RADIAN PER DEGREE/24X * 7F 12.5 ) AER2340 

FORMAT (//25X* 10HADOITIONAL L0ADING/24X*24HWITh CL BASED On S ( TRUE aER2 350 
1)71X*11H-AT CL 0ES-/67X.34HL0A0 DUE ADO. LOAD AT BASIC LOA03X * 27 AER2360 
2HSPAN LOAD AT X LOCATON 0F/8H STATI0N6X*5H 2Y/B9X.9H SL COEF *4XAER2370 
3*8HCL RAH04X.7HC R AT 1 0 , 7X * 1 4HT0 TWIST CL = *F9.5 * 3X * THAT CL = 05X , 2 AER2380 


46HDESIREO CL LOCAL CENT PR/) 
FORMAT (/47X.61HC0NTRIBUTI0N OF 


The SECOND planform to span load 


AER2390 

OAER2400 
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320 

330 

340 

3S0 

3b0 


370 

380 

390 


ARE 


AER2410 

AER2420 

AER2430 

CaER2440 

AER2450 

AER2460 


1 I5TRIBUTION/) 

format l4,F12.5»5X»3F12.5f3X.3F12.5»3X,2Fl?.5) 

FORMAT (/54X,2<*HCM3 AND CLQ ARE COMPUTED//) 

FORMAT (/38X,57HSTaTic LONGITUDINAL AERODYNAMIC COEFFICIENTS 
IOMPuTED//) 

FORMAT </59X,15HCL^ IS COMPUTED//) 

FORMAT (/20X,lHXfnx»lHXtllXflHY»lIXtlHZfl2Xt 1hS*5X»9HC/4 SWEE°t AX AER 2470 
1 tSHDIHEDRAL *2X. 1 IHLOCAL ALPHA t 2X » 1 9HDElT A CP AT DES IRED/ I 9X * 3HC/4 . AER2aB0 
29X »4H3C/^*A2X»5H ANGLE* 7X * 5H ANGLE ♦ AX » lOHiN R AO I ANS * AX t AHCL =*F10.5/AER2a90 

AER2500 

(I2X,RF12.5) AER2S10 

(^/////////36X*AHCLP=»F9.5////> AER2520 

(//////////A2X. AHCM0=,F9.5, 10X*AHCLO=*F9.5////) AER2S30 

AER25A0- 


3) 
FORMAT 
FORMAT 
FORMAT 
END 


48 



APPENDIX C 



SUBROUTINE FTLUP ( X . Y * M ♦ N * V AR I , V ARO) 



TLU 

10 

c 

•“•document date 09-12-69 SUBROUTINE REVISED 

07-07-69 


20 

c 

MODIFICATION OF LIBRARY INTERPOLATION SUBROUTINE FTLUH 

TLU 

30 


DIMENSION VARKD* VARD(1)» V(3t» YY(2) 



TLU 

40 


DIMENSION 11(43) 



TLU 

SO 

c 




TLU 

60 

c 

Initialize all interval pointers to -i.o for 

monotonicity checktlu 

70 


DATA ( 1 1 ( J) . J=1 *43)/43“-l/ 



TLU 

80 


MA=IA8S(M) 



TLU 

90 

c 




TLU 

100 

c 

ASSIGN interval POINTER FOR GIVEN VARI TABLE 


TLU 

110 

c 

THE same pointer WIlL BE USED ON A GIVEN VARI 

TABLE EVERY TIME 

TLU 

120 


LI=M0D(L0CF (VARI ( 1 ) ) *43) *1 



TLU 

IJO 


1=11 (LI) 



TLU 

140 


IF (I.GE.O) GO TO 60 



TLU 

150 


IF (N.LT.2) go TO 60 



TLU 

160 

c 




TLU 

170 

c 

MONOTONICITY CHECK 



TLU 

180 


IF (VARI (2)-VARI ( 1 ) ) 20*20*40 



TLU 

190 

c 

ERROR IN HONOTONICITY 



TLU 

200 

10 

K=L0CF(VAWI (1) ) 



TLU 

210 


print 170* J,K* ( VARI ( J) * J=1 *N) * ( VARD( J) * J=1 *N) 



TLU 

220 


STOP 



TLU 

230 

c 

MONOTONIC DECREASING 



TLU 

240 

20 

DO 30 J=2*N 



TLU 

250 


IF (VARI (J)-VARI ( J-1 ) ) 30*10*10 



TLU 

260 

30 

CONTINUE 



TLU 

270 


GO TO 60 



TLU 

280 

C 

MONOTONIC increasing 



TLU 

290 

40 

DO 50 J=2*N 



TLU 

300 


IF (VARI (J)-VARI (J-1 ) ) 10*10*50 



TLU 

310 

dO 

continue 



TLU 

320 

C 




TLU 

330 

c 

INTERPOLATKON 



TLU 

340 

60 

IF (I.LE.O) 1=1 



TLU 

350 


IF (I.GE.N) I=N-1 



TLU 

360 


IF (N.LE.I) go to 70 



TLU 

370 


IF (MA.Nt.O) GO TO 80 



TLU 

380 

c 

ZERO ORDER 



TLU 

390 

70 

Y=VARD(1 ) 



TLU 

400 


GO TO 160 



TLU 

4l0 

c 

LOCATE I INTERVAL ( X ( I ) . lE . X .LT . X ( I * 1 ) ) 



TLU 

420 

80 

IF ((VARI(I)-X)*(VARKI*1)-X) ) 110*110,90 



TLU 

430 

C 

IN GIVES DIRECTION FOR SEARCH OF INTERVALS 



TLU 

440 

90 

IN=SIGN( 1 .0* (VARI ( 1*1 ) -VARI ( I ) ) *(X-VARI ( I) ) ) 



TLU 

450 

C 

IF X outside ENDPOINTS* EXTRAPOLATE FROM END INTERVAL 


TLU 

460 

100 

IF ( (I*IN) .LE.O) GO TO 110 , 



TLU 

470 


IF ( ( 1»IN) .GE.N) GO TO 110 



TLU 

480 


I = I*IN 



TLU 

490 


IF ( (VARI ( I )-X)* ( VARI ( 1*1 ) -X) ) 110*110*100 



TLU 

500 

110 

IF (MA.EO.Z) 60 TO 120 



TLU 

sio 

c 




TLU 

520 

c 

first order 



TLU 

530 


Y=(VARD( I )• (VARI ( 1 ♦ 1 ) -X ) -V ARD ( I ♦ 1 ) * ( VAR I ( I )-X) )/( VARI (I* 

1 )-VARI ( I ) TLU 

540 


1 ) 



TLU 

550 


00 TO 160 



TLU 

560 

c 




TLU 

570 

c 

SECOND 0*^0ER 



TLU 

580 

120 

IF (N*EQ.2) GO TO 10 



TLU 

590 


IF (I.EQ. <N-1) ) GO TO UO 



TLU 

600 
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IF JI.EO.l) GO TO 130 

TLU 

610 

c 

PICK thi»<d point 

TLU 

620 


Sk=VARI (1*1) -VARI ( I ) 

TLU 

630 


IF ( (SK*(X-\/ARI ( I-l ) ) ) .LT. (SK* (VARI ( I*2)-X) ) ) GO TO 140 

TLU 

6^0 

130 

L = I 

TLU 

650 


GO TO 150 

TLU 

660 

140 

L=I-1 

TLU 

670 

150 

V(1)=VAHI (L)-X 

TLU 

680 


V(2)=VARl (L*l )-X 

TLU 

690 


V(3)=VARl (L*2)-X 

TLU 

700 


YY(1)=(VARD(L)*V(2)-VARD(L*1)*V( 1 ) )/(VARI (L*1)-VARI (L) ) 

TLU 

710 


YY(2)= (VARD(L*1 )*V (3)-VARO(L*2)*V(2) )/(VARI (L*2)-VARI (L*l ) ) 

TLU 

720 


Y=(YY(l)»V(3)-rY(2)*V(l) ) / ( VAH i ( L*2) -V AR 1 (L) ) 

TLU 

730 

160 

II(Ll) = r 

TLU 

740 


RETURN 

TLU 

750 

C 


TLU 

760 

c 


TLU 

770 

170 

FORMAT (iHlfBOH TABLE BELOW OUT OF ORDER FOR FTLUP AT POSITION 

♦ TLU 

780 


1I5,/31H X TABLE IS STORED IN LOCATION ,06,// (8G15.8) ) 

TLU 

790 


END 

TLU 

800- 
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SOBROUTI'^E CDICLS ( AP ♦ ARTRUE ♦ I SE'<5P»MT0T »NSV «CDI tCOI T ) 
dimension ETAN(51). GAMPR(51.1)t ETA(41). GAMMA(41), Vt(41), t 

^COMMOn''/ALL/ 80T,M,BETA,PTE5T,QTeST,TBlSCW(50).0(200) .PN(2oO) 

100) ♦ ALP (200) ,S (200) tPSI (200) *PHl (50) »2H(50) 

COMMON /IHRECOI/ SlOAD(3»50) 

00 10 1=1.41 

DO 10 J=1.41 

FVM(I.J)=0 

SPAN=2.*B0T 

CAVB=5PAN/ARTRUE 

PI=.314159265E*01 

NST=lSEMSP.l 

NN=MT0T 

DO 20 N=1»ISEMSP 

NM=NSV-N 

NSCW=TBLSCK(NM) 

NN=NN-NSCW 

ETAN(N) =ASIN (-0 (NN) *2./SPAN) 

GAMPR (N. 1 ) =SLOAD (3.NM) *CAVB/ ( 2.*SPAN) 

CONTINUE 

ETAN(NST)=PI/2. 

GAMPR (NST ♦ 1 ) =0 
DO 30 NP=1.41 
ANP=NP 

ETA(NP) = ( ANP-21 . ) “RI/42. 


DO 40 JK=21.41 

CALL FTLUP ( eta ( JK )» gamma ( JK ). 1 .'1ST »ET an. GAMP f( ) 

CONTINUE 

DO 50 NY=22.4l 

ETA(NY)=SIN(ETA(NY) ) 

NR=42-NY 

ETA(NR)=-ETA(NY) 

GAMMA(NR)=GAMMA(NY) 

DO 90 NU=21.41 
ANU^NU 

DO 80 M=l.41 
AN=N 

NNUD=IAB5(N-NU) 

VE (N) =C05 ( ( ( AN-21 . ) *Pl ) /42. ) 

IF (NNUD.NE.O) GO TO 60 

B(N)=(42.)/(4.0*COS( ( (ANU-21. )*PI»/42,) ) 

GO TO 80 

IF (MOD(NNU0.2) .EO.O) GO TO 70 

fa(N)=VE(N)/((42.)*(ETA(N)-ETA(stU))**2) 

GO TO 80 
B(N)=0.0 
CONTINUE 
DO 90 NP=21.41 
NUST=I ArtS(NU-21 ) 

IF (NUST.EQ.O) GO TO 90 
IF (MOD(NUST,2) .EO.O) GO TO 90 
NPST=IABS (NP-20) 

IF (MOD(NPST.2) .EO.O) GO TO 90 
NPNUD=IAB5(NP-NU) 

IF (NPMUO.EQ.O) GO TO 90 
IF (M0D(NPNUD.2) .EO.O) GO TO 90 

F VM (NU .NP )=2.0*8(NP)/21.*COS((ANU-21.)*PI/42.) 

IT=42-NU 
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90 

C 

100 

no 


ITT=42-niP 

CDI 

610 

FVN(NU»ITT)=2.0*B(ITT)/21.*COS< <AnU-21.)»P1/42.) 

COI 

620 

FVN( IT*NP) =rvN (NU* ITT ) 

CDI 

630 

FVNMIT.lTT)=rVN(NU.NP) 

CDI 

6^0 

continue 

COI 

650 


CDI 

660 

CCC=0,0 

CDI 

670 

DO 100 N=lt41 

CDI 

680 

CCC=CCC* (GAHhA(N)«GAmMA(N) ) 

CDI 

690 

DO no NUP = 1,41 

CDI 

710 

DO 110 N=lt41 

COI 

720 

CCD=CCD-2.0*FVN (NUP*N) * ( GAMMA ( NUP) *GAMMA (N) ) 

CDI 

730 

Continue 

CDI 

740 

CDI=PI*AP/4.* (CCC*CCO) 

COI 

750 

CDIT=1./(PI*AR) 

CDI 

760 

RETURN 

COI 

770 

END 

CDI 

780- 
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I 


c 

c 


10 

20 


OVERLAY (ml I NGTL*^«0) 

PROGRAM CDRAGNE 

dimension GAM(IOOO)* XCAOOOO)* YO(IOOO)* CCRI20 
1XCC(20)« CCC(200)» CRR(200)* Ya<50)« CRI(51)» NM 
2HD(50)» XC^A(50)» YY(2)* PPHI(50t« ZZH(SO)* Z(10 
3A(50)t SSA(IOOO)t ALOP(200). AlL°(50)* ALPPOUOO 
A), YQO(50) 

COMMON /All/ 0OT.M.BETA,PTESTOTEST.TBlSCW( 5O> ,0 
100)*ALP(200),S(200)»PSI(200)*PHI(50)»ZH(50) 
COMMON /ONETHPE/ T W I ST ( 2 ) • CREE * S^EF » C A VE» CLDES » S 
IDHT (2) *CONF IG«NSSW5V(2) .MSV(2) fKBOT »PL ANt IPLANtM 
2) ,XT<2) *CLWB.CMCL*CLA(2) . BL A IP ( 50 ) . CL AMAR < 2 ) 
COMMON /TOTHRE/ CIP(200*2) 

COMMON /INSUB23/ A=S I . APH I . XX . Y Y Y . ZZ t SNN* TOLCSO 
COMMON /THREFOR/ CCAv ( 2 ♦ 50 ) . CLT « CLNT *NS5W , ALPO 
COMMON /CCRROD/ TSP AM * TS® AN A » KB I T 


)♦ FW(2)» FV(2)* XDRA 
A(2)» XCCA(200)* CORA 
00>« PHIIdOOO)* SORA 
0)* Al 0(20»* rC(51DRA 
DRA 

(200) tPN(200) .PVt?ORA 
DRA 

TRuf*AR* ARTR uEtRTCoRA 
ACH»SSWWA (50) .XL (20RA 
DRA 
DRA 
DRA 
DRA 
DRA 
DRA 
DRA 


WRITE (6.250) 

APSI=TOLC5Q=TBLS=0. 

PI=A.*ATAN( 1 . ) 
fpi=a.»pi 

eOTL=ABS (TSPAN) 

BOL=ABS(TSPAnA) 

SNN=BOTL/ (2.*N55WSV (KBOT) ) 

DELTYB=2.*SNN 
NMA(KBOT)=0OTL/OELTYB 
NMA(KBIT)=B0L/0ELTYB 
NMAX=NMA ( 1 ) *NMA (2) 

DO 10 I=1.M 

CRRm=ClR(I.l) ♦CIR(I.2)»(CL0 ES-CLT)/ClNT 

CONTINUE 

5CWMIN=20. 

DO 20 I=1.N5SW 

SCWMIN=AMIN1 (SCWMIN.T8LSCW( I ) ) 

NSCWMIN=SCWMIN 

MM=NSCWMlN»NMAX 

DELTX0C=1 ./SCWMIN 

DO 100 LA=1 .NSSW 

CHO(LA)=CCAV (2.LA)*CAVE/BETA 

DELTXX=1 ./TBlSCW(LA) 

XC=-.75*DELTXX 
ITBL=TBlSCW (lA) 

DO 30 LB=1.ITBL 
XC=XC*DELTXX 
XXCC (LB) =XC 

lc=lb.tbls 

AL0(LB)=ALP(lC) 

XlE=PN<LC) *Ch0(LA) * (1 ,-.75/TBLSCW(LA) ) 

XOC=-.75»DELTXOC 

KCODE=LB=0 

DO 90 K=1.NSCWMIN 

j=K*(LA-1)*NSCWMIN 

X0C=X0C*DELTX0C 

XCCA(J)=-XOC*ChO(LA) *xle 

call FTLOP (XOC.ALOP(J) ..I.ITHL.XXCC.AlO) 
AXMN=K*D£LTX0C 

cat=o. 

IF (KCOOE.EQ.2) CAT=CCR(l6)-CUT 
KC0DE=0 


ORA 190 
DRA 200 
DRA 210 
DRA 220 
DRA 230 
DRA 240 
DRA 250 
DRA 260 
DRA 270 
DRA 280 
ORA 290 
DRA 300 
DRA 310 
DRA 320 
DRA 330 
DRA 340 
ORA 350 
DRA 360 
DRA 370 
DRA 380 
ORA 390 
ORA 400 
DRA 4l0 
DRA 420 
ORA 430 
DRA 440 
ORA 450 
DRA 460 
ORA 470 
DRA 480 
ORA 490 
DRA 500 
DRA 5l0 
DRA 520 
DRA 530 
ORA 540 
DRA 550 
DRA 560 
ORA 570 
ORA 580 
DRA 590 
ORA 600 
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40 

50 

60 

70 

80 

90 

100 


no 

120 

130 


L0=LB*1 

ORA 610 

UC=La*T8LS 

ORA 620 

CCR(L8)=CHR(LC) 

DRA 630 

AXITBL=LB*DELTXX 

ORA 690 

IF (AXHN-AXITBL) 50t60.70 

DRA 650 

CUT=CCR (LB)* ( AXMN-(LB-l) *DEL TXX ) /OELTXX 

ORA 660 

KC0DE=2 

ORA 670 

GO TO 00 

ORA 680 

KCODE=l 

ORA 690 

CUT=CCR (LB) 

ORA 700 

CAT=CAT*CUT 

DRA 710 

IF (KCOOE.GE. 1 ) GO TO 90 

ORA 720 

IF (LB.LI.ITBL) GO TO ^*0 

ORA 730 

CCC(J)=CAT 

DRA 790 

TBLS=T8LS*T8LSCW(LA) 

ORA 750 

CONTINUE 

DRA 760 

11 = 1 

ORA 770 

DO 150 I=1,IPLAN 

ORA 780 

BOTT=BOTl 

ORA 790 

IF (I.EQ.KBIT) B0TT=B0L 

ORA 800 

IUZ=NSSW5V ( I ) 

ORA 010 

IUX=IUZ*1 

DRA 020 

IC=MSV(lj*(I-l)*MSV(2) 

ORA 830 

ID=iC*l 

ORA 890 

IZ=NSSWSV( 1 ) ♦ ( I-l ) *NSSWSV(2) 

ORA 850 

VCAT=0. 

ORA 860 

IAMM=NMA(I) 

ORA 870 

DO lAO LA=1*NSCWMIN 

ORA 880 

YC(l)=-PI/2. 

ORA 890 

CRI (1 )=0. 

ORA 900 

DO 120 J=1.IUZ 

ORA 910 

L = J»1 

ORA 920 

LU=LA» ( J-1 ♦ ( I-l ) *N5SWSV ( 1 ) ) *NISCWNIN 

ORA 930 

ALLP(J)=AL0P(LU) 

ORA 960 

XCAA(J)=XCC'»(LU) 

ORA 950 

CRI (L)=CCC(LU) 

ORA 960 

IF (LA.NE.l) 60 TO 120 

ORA 970 

JJ=J*(I-1)»NSSWSV(1) 

DRA 980 

ZZH( J)=Zri( JJ) 

ORA 990 

SA(J)=SS»IWA( jj) 

DRAlOOO 

PPHI ( J) =PHI ( JJ) 

DRAlOlO 

YQQ(J)=Q(II) 

ORA1020 

II = II*T8LSC«MJJ) 

ORA1030 

IE=IUZ-J*1 

ORA109O 

ITl=T8LSC»i ( IZ) 

ORAIOSO 

IO=ID-ITL 

DRAl 060 

IA=ID*1TL 

ORA1070 

IF (IA.gT.IC) YCAT=YcAT-S(IO) 

0RA1080 

IF (IA.gT.IC) go to 110 

ORA1090 

YCAT=YCAf-S(IO)-S(lA) 

DRAl 100 

IZ=IZ-1 

DRAl 110 

YB(IE)=YCAT 

ORAl 120 

CONTINUE 

DRAl 130 

DO 130 JP=1,IUZ 

DRA1190 

JZ=JP*1 

ORAl 150 

YC(JZ)=A5IN(YB( JP) /BOTT) 

ORAl 160 

CONTINUE 

DRAl 170 

YOB=-NMA( I )*2.*SNN-SNN 

ORAl 180 

DO lAO K=1.IAMM 

ORAl 190 

KP=CA* (K-1* ( I-l )*NHA(1 ) ) *NSC«MIN 

ORA1200 
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Y0B=Y08*i)ELTr8 

Y0C=ASIN(y0B/80TT) 

CALL TTLUP (Y0B*Y(J(KP) .♦1*IUZ*YB*Y00) 
call FTLJP (Y0R»AL3PD(KP) .♦1.IUZ.Y0»ALlP) 
call FTLBP (YOBtSSA(KP) .»1.IUZ,YB,5A) 

CALL FTLOP (YORtXCA (KP) . *1 » lUZ* YB.XC^^) 

CALL FTLiJP ( YOB»Z (<P) » *1 f IUZ» Y8*ZZH) 

CALL FTLBP (YOfl.PHII (KP) ,*1.IUZ»YB»PPHI) 

CALL FTLUP ( YOC * GA8 ( KP ) * ♦ 1 » I UX , YC * CR I ) 

IF <Y0B.GT,YB( lUZ) ) GAM (KP) =CWI ( lUX) 

I'fO CONTINUE 
IGO CONTINUE 

CDRAG=CTHRUST=CSUCT=0. 

C0NST=16.*SNN»B0T/5REF 
DO 190 LI=1*NMAX 
LA=(L1-D*NSCWMIN*1 
LB=LI*NSCWMIN 
CDRAGIT=CTT=0. 

DO 100 NV=LA,LB 
CPT=COS(ATAN(PhI I (NV) ) ) 

velin=o. 

DO 170 NN=1.MM 
XX=XC4(NY>-XCA(NN) 

YY<1)=Y0<NV)-Y0(NN) 

YY(2)=YQ(NV) ♦YO(NN) 

2Z=2(NV)-Z«NN) 

APHI=ATAN (Pril I (NN) ) 

DO 160 I=l»2 
YYY = YY ( I ) 

call INFBUB (BOTtFY ( I ) *FW ( I ) ) 

APHI=-APHI 

160 CONTINUE 

VFLIN={ (FWU ) ♦FW (2) ) -(FV ( 1 ) *FV (2) ) *PHI I (NV) ) *GAM (NN) /FPI ♦ VfL IN 
170 CONTINUE 

CTT=CTT*GAM(NV) * ( AlP0/57.2957«*ALPP0(NV) )*CPT/(2.*B0T) 

100 CDRAG1T=CDRA61T*VELIN*GAM(NV)*CPT/(2.*B0T ) 

CTT=CTT-C0RAG1T 
SWLE=ATAN(SSA(LA) ) 

CST = CTT/COS(Sv^LE) 

CCC(LI)=C0RAGIT 

CRR(LI)=CTT 

XCC4(LI)=C5T 

CDRAG=CDRAG*C0RAGIT*C0NST 
CThRUST=CTHRUST»CTT*CONST 
CSUCT=CSUCT*CST "CONST 

190 continue 

T8LE=II=0 
LI = 0 
LBLR=0 

DO 220 1=1.2 
IAMH=NHA( I ) 

00 200 J=1.IAMM 
JJ=J*(I-1)*NMA(1) 

LA=1*( J-1* ( I-l )*NMA( 1 ) ) "NSCWMIN 
GAM( J)=CCC(JJ) 

XC4(J)=CRRUJ) 

2(J)=XCCA(JJ) 

PHII (J)=YQ(LA) 

200 CONTINUE 

IU2=NSSW5Y(I) 


DRA1210 

DRA1220 

ORA1230 

ORA1290 

ORA1250 

DRA1260 

ORA1270 

ORA120O 

ORA1290 

DRAI 300 

ORA1310 

DRA1320 

DRA1330 

DRA1390 

DRA13S0 

DRA1360 

DRaI 370 

0RA1380 

DRA1390 

ORAl AOO 

nRAl4lO 

DRA1420 

0RA1430 

ORA1440 

DRA1450 

DRA1460 

DHA1470 

DRA140O 

ORA1490 

ORA1500 

DRA1510 

ORA1520 

DRA15J0 

ORA1540 

DRA15B0 

ORA1560 

OR A 1 570 

DRA1580 

ORA1590 

DRA1600 

DRA1610 

DRA1620 

ORA1630 

DRA1640 

DRA1650 

DRA1660 

DRA1670 

DRA16B0 

DRA1690 

DRA1700 

DRA1710 

ORA1720 

ORA1730 

ORA1740 

DRA17G0 

0RA1760 

ORA1770 

DRAI780 

ORA1790 

ORA1800 


55 



APPENDIX C 


^10 

^^0 

230 


C 

C 

c 

Z>*0 

250 


260 

270 


DO 210 LBLAIR=1 ♦ lUZ 

oRAieio 

LI=LI^1 

DRA1820 

LU = UT8LE 

DRA1830 

LBLR=L8L^*1 

DRAI840 

YBB=Q(LU) 

DRA1850 


0RA1860 

TBLE=TBLE^TBLSCW(II) 

ORA1870 

YOOB=Y08/BOT 

DRAI880 

CALL FTLOP { YBB»COf?AGlT, ♦ lAMHfPHi I fGAM) 

DRA1890 

CALL FTL*JP (YB0»CTT»*l,IAMMfPHII»XC4) 

DRA1900 

call FTLJP (YBB»C5T,^1,IAHM»PHII,Z) 

DRAI910 

LL=L0LA I ( I -1 ) *NS5WSV ( 1 ) 

DRA1920 

SWALE=ATAN<SSWWA(LL) )*57. 29578 

DRA1930 

IF ni.EQ. (NSSWSVU ) ♦!) ) WRITE <6f240) 

DRA1940 

WRITE (6»260) L I » YDOB f SW ALE f CDRAGI T ♦ CT T » CST 

DRA1950 


BLAIR(LBI-R)=CST 


DRA1960 


GO TO 230 


CONTINUE 

IF (NSSWSV (2) .EO.O) 

CONTINUE 

CDCL2=C0R AG/CL DES»»2 

NRITE (6*270) COCL2.CTHRUST.CSUCT 


DRA1970 
ORA1980 
DRA1990 
ORA2000 
DRA2010 
ORA2020 
0RA2030 
ORA20^0 
ORA2050 

OF THE SECOND PLANFORH TO THE CHORD OORA206O 

ORA2070 

LEADING EDGE THRUST AND SuCTION CDRA2080 
IOEFFICIENT CHARACTERISTICS/40X,53 hCOMPJTED at THE DESIRED CL FROM DRA2090 
2A near field SOLUT ION//58X,20HSECTION COEFF IC I ENTS/48X . 1 1 Hl . E. SWDRA2100 
3EEP/25X, /HSTaTI0N*9X,5H 2Y/B . 5X * SHANGLE . 7x , 9HCDI I C/2B*5X*7HCT C/2DRA2110 
4B*5X*7HCS C/2B) 0RA2120 

FORMAT (20X,I10*5X,5F12.5) DRA2130 

FORMAT (///57X , 18HTOTAL COEFF I C I ENTS//36X , 1 2HC0I I/CL"«2 =F 1 0 .5* SX * QRA2 1 40 
13HCT=*F1 0.5, 5X,3HCS=,F 10.5) ORA2150 

end ORA2160- 


RETURN 


FORMAT 
IR DRAG 
FORMAT 


(/37X,62HC0NTRIBUTI0N 
FORCE/) 

(////30X,73HINDUCEO DRAG* 
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c 

c 


c 

c 


c 

c 

c 


c 

c 

c 

10 


c 

20 

30 

c 

40 

50 

C 

C 

60 


C 

70 

C 

BO 

C 

90 

C 

100 


110 

c 

c 


c 

c 

120 


SU0ROUTI^<E TTLUP (X* Y*M*S*VARI .VARD) 'LU 

•••DOCUMENT DATE 09-12-69 SUBROUTINE REVISED 07-07-69 »«***o**«tLU 
modification of library INTER°OLAT10N SUBROUTINE FTlUP TLU 

DIMENSION VARI(1)» VARD(I)* V(3)* YY(2) TLU 

dimension 11(43) TLU 


10 

20 

30 

40 

so 

60 


initialise all interval pointers to -1.0 

DATA ( I I (J) » J=1 *43) /43*-1/ 

MA=IABS(M) 


FOR MONOTONICItT CHECKTLU 70 

TLU 00 
TLU 90 
TLU 100 


ASSIGN INTERVAL POINTER FOR GIVEN VARI TABLE 
THE SAME pointer MILL BE USED ON A GIVEN VARI TABLE EVERY TIME 
LI=MOD(LOCF (VARI ( 1 ) ) *43) *1 
I = IKLI) 

IF (I.GE.O) GO TO 60 
IF (N.LT.2) GO TO 60 


TLU 110 
TLU 120 
TLU 130 
TLU 140 
TLU 150 
TLU 160 
TLU 170 


MONOTONICITY CHECK 
IF { VARI (2) -VARI ( 1 ) ) 20.20*40 
ERROR IN MONOTONICITY 
KsLOCF ( VARI ( 1 ) ) 

PRINT 170* J,K * ( VARI ( J) * J=1 *N) * ( VARDI J) * J=1 *N) 
STOP 

MONOTONIC DECREASING 
DO 30 J=2*N 

IF (VARI ( J)-VARI ( J-1 ) ) 30*10*10 

CONTINUE 

GO TO 60 

MONOTONIC INCREASING 
DO 50 J=2*N 

IF (VARI (J)-VARI (J-1 ) ) 10*10*50 

continue 


TLU 100 
TLU 190 
TLU 200 
TLU 210 
TLU 220 
TLU 230 
TLU 240 
TLU 250 
TLU 260 
TLU 270 
TLU 200 
TLU 290 
TLU 300 
TLU 310 
TLU 320 
TLU 330 


interpolation 

IF (I.LE.O) 1=1 
IF (I.GE.N) I=N-1 
IF (N.LE.l) GO TO 70 
IF (MA.nE.O) go to 00 
SERO ORDER 
Y=VARD( 1 ) 

60 TO 160 

LOCATE 1 interval ( X ( I ) . LE . X . L T . X ( I ♦ 1 ) ) 

IF { (VARI(I)-X)*(VARI (I*D-X) ) 110*110.90 

IN GIVES DIRECTION FOR SEARCH OF INTERVALS 
IN=SI6N( I .0* ( VARI ( I»1 ) -VARI ( I ) )* (X-VARI ( I ) ) ) 

IF X OUTSIDE ENDPOINTS* EXTRAPOLATE FROM END INTERVAL 
IF < (I.IN) .LE.O) GO TO 110 - 
IF ( ( I» IN) .GE.N) GO TO 110 
I=1.IN 

IF ( ( VARI ( I ) -X) • ( V ARI ( I ♦ 1 ) -X) ) 110.110*100 
IF (MA.EO.2) go to 120 


TLU 340 
TLU 350 
TLU 360 
TLU 370 
TLU 300 
TLU 390 
TLU 400 
TLU 4l0 
TLU 420 
TLU 430 
TLU 440 
TLU 450 
TLU 460 
TLU 470 
TLU 400 
TLU 490 
TLU 500 
TLU SlO 
TLU 520 


FIRST ORDER 

Y=(VAR0( I ) * ( VARI ( I ♦! ) -X) -VARD( 1*1 ) 


•(VARI <I)-X) )/(VARI (I 


1 ) 


GO TO 160 


TLU 530 
♦ D-VARI (I)TLU 540 
TLU 550 
TLU 560 
TLU 570 


SECOND ORDER 

IF (N.EQ.2) GO TO 10 

IF (I.EQ.tN-D) 60 To 140 


TLU 500 
TLU 590 
TLU 600 
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ir (I.EQ.n GO TO 130 
C PICK ThIPD point 

SK=VARI(I*1)-VARKI) 

IF ( (SK* (X-VAP I n-1 ) ) ) .LT. (SK» ( VARI < 1*2) -X) ) ) GO TO lAO 
130 L=I 

GO TO 150 
140 L=I-1 

150 V(1)=VARI (U-X 

V(2)=VARI (L*l)-X 
V(3)=VARI (L*2)-X 

YY ( 1 ) =(VAR0 (L) *V (2) -VARD(L*1 ) *V(l))/<VARI<L*l)-VAHI(L)) 
YY<2)=(VARD(L*1)*V(3)-VAR0(L*2)*V(2 ) )/(VAri (L* 2)-VARI (L*l) ) 
Y=<YY(1>*V(3)-YY(2)*V(1))/<VARI(L*2)-VARI(L)) 

160 II(LI)=r 
RETURN 

70 FORMAT (iHl.SOH TABLE BELOW OUT OF ORDER FOR FTlUP AT POSITION 
115. /31H X TABLE IS STORED IN LOCATION .06.// (8G15.8) ) 

END 


TLU 610 
TLU 620 
TLU 630 
TLU 640 
TLU 650 
TLU 660 
TLU 670 
TLU 680 
TLU 690 
TLU 700 
TLU 710 
TLU 720 
TLU 730 
TLU 740 
TLU 750 
TLU 760 
TLU 770 
.TLU 780 
TLU 790 
TLU 000- 
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I 



10 


20 


OVERLAY (KiflNGTLt5» 0) 

dimension YY(2)* WVOU(<^0)* FV(?)* rw(2)f XTLEG(60)f CIRSUH(50)i 
1EGSV<50)* ZlEGSV(50) _ 

COMMON /All/ bot,m, beta, pTEST.QT5ST,TBlSCw<50). 0(200). PN(^00).► 
100) , alp (200 ) ,S (200) ,DSI (200) .PHI (50) ,ZH(50) 

COMMON /TOTHOE/ CIP(200.2) .oto,.p 

COMMON /OMETHRE/ TK I ST ( 2 ) . CREF , SPFF , C A vE , CLDES. STRUE , AR , APtPiJE ; 

lOHT (2) .COnF IG.NS5*lSV(2).MSV(2).KBOT.PLAN.IPLAN.MACH,5S.(WA(bO).. 

2) .XT (2) .CLWB.CmCL » CLA ( 2) .BlAIW (50 ) .CLAMAR (2) 

COMMON /THREFOR/ CC A V ( 2 , 50 ) . CL T . CLNT .NSSirf. ALPO 
COMMON /INSUB23/ APS I . ApH I , XX , Y Y Y , 22 .SNN . TOLCbO 
DIMENSION XKVSEW(2). CENTR(2) 

XKVSEW ( 1 ) =XKVSE*( (2) =CENTR ( 1 ) =CENTR ( 2 ) =0 . 0 

IF ( IPLAN.EO. 1 . AND.Xl ( 1 ) .EO.XT ( 1 ) ) GO TO 540 

IF (XL ( 1 ) .EQ.XT ( 1 ) . AnD.XL (2) .EQ.XT (2) ) GO TO 540 

BLAMAR=1 ./BETA 

XT(1)=XT(1)»BLAMAR 

XT (2) =XT (2) »BLAHAR 

XL(1)=XL‘1 )*BLAMAR 

XL(2)=XL(2)*BLAMAR 

THE TOLERANCE SET AT THIS POINT IN THE PROGRAM MAY NEED TO BE 
changed FOR COMPUTERS OTHER THAN THE COC 6000 SERIES 

TOLC=.0100«BOT 

T0LC50=TDLC*T0LC 


TIPSU=PITCH=0. 
NS5W=NSSMSV ( 1 ) .NSSRSV (2) 


GEOMETRY FOR TIP TRAILING LEGS 


ITT = 1 

1M=IHM=NSSw1=0 

CCIRS=0. 

NSSW2=NSSW3 = NS5«SV ( 1 ) 

L=1 

NSCW=MSV(1)/NSSWSV(1) 

GO TO 20 
NSSMl=NSSl«(SV ( 1 ) 
NSSW2=NS5WSNSS«3=NSSii(SV (2) 
L=NSSWSV(1) *1 
NSCW=MSV (2) /NSSMSV (2) 

IF (XL (2) .EQ.XT (2) ) GO TO 500 

I=IMM.l 

J=IMM»2 

1UU=2 

APHI=ATAN(PHI ( IM*1) ) 
SA=SIN(APHI ) 

CA=COS(APHI) 

TLX1=PN(I)-5(I)*TAN(PSI (I) )*CA 

TLX2=PN(J)-S( J)*TAN(PSI ( J) )*CA 

CLFTLG=T(-Xl-TLX2 

XTLEG( 1 ) =TLX1/2.*TLX2/2. 

VLEG=Q(D-S( I)*CA 

IE (NSSWl.EQ.O) YLEGSV ( I ) =YLEG 

2LEG=2H ( IM* 1 ) -S ( I ) *SA 


TIP 360 
TIP 370 
TIP 330 
TIP 3R0 
TIP 400 
TIP 410 
TIP 420 
TIP 430 
TIP 440 
TIP 450 
TIP 460 
TIP 470 
TIP 480 
TIP 490 
TIP 500 
TIP 510 
TIP 520 
TIP 530 
TIP 540 
TIP 550 
TIP 560 
TIP 570 
TIP 580 
TIP 590 
TIP 600 
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IF (NSSWUEO.O) ZLEGSV( 1)=ZLEG 
IF (XL(ltT) .EQ.XT(ITT) ) GO TO 100 
DO 30 NV=2.NSCW 

nvt=mv-i 

30 XTLEG(NV)=XTlEG(NVT) -CLFTLG 

NCTL=0 
NA = l 
NB=NSCW 

40 DO 70 NV=NA»nB 

c 

c 

c The ratio of w/u is initialized to -i because in The term 

C -U*ALPHA/UfUSED IN ThIS SUMMATlONf ALPHA IS SET TO 1 RADIAN 

C SO THAT the RESULTING TIP SUCTION CAN BE USED DIRECTLY TO FIND 

C KV SIDE EDGE 

C 

c 

WVOU(NV) . 

IZ = I 

NNN=T9LSLW(IZ) 

DO 60 NN=1»M 
APHI=ATAN(PHI ( IZ) ) 

APSI=PSI (NN) 

XX=XTLEG<NV)-PN(NN) 

YYU)=YLEG-Q(NN) 

YY(2)=YLEG^Q(NN) 

ZZ=ZLEG-^H( IZ) 

Snn=S (NNJ 
DO 50 I = U2 
YYY = YY ( I ) 

CALL INFSUB { QOT » F Y ( I ) f Fl^ ( I ) ) 
aphi=-aphisapsi=-apsi 
SO CONTINUE 

WVOUCNV)=WVOU(NV) ♦ (FW( 1 ) ♦FWCZ) ) *C I R C NN ♦ 2 ) / I 2 • 5663704 
IF (NN.LT.NNN.OR.NN.EQ.M) GO TO 60 
IZ = IZM 

NNN=NNN^TBLSCW (IZ) 

60 CONTINUE 

70 continue 

NCTL=NCTL>1 

IF (NCTL-2) Q0»100»150 
C 

C GEOMETRY FOR SPANWISE BOUND VORTICES 

C 

bO NA=NSCw^i 

NB=2<^NSC«^ 

JA=IMM4l 

YLEG=Q(JA) 

ZLEG=ZH( IM^l) 

DO 90 J=i*NSCW 

JK=IMM^J 

NV=J^NSC^ 

90 XTLEG(NV)=PN( JK) 

60 TO 40 
C 

C geometry along RIGHT TRAILING LEGS 

C _ 

100 NA=2*NSC«f^l 

N8=3*NSCrt 
CCIRS=0* 


TIP 610 
TIP 620 
TIP 630 
TIP 640 
TIP 650 
TIP 660 
TIP 670 
TIP 680 
TIP 690 
TIP 700 
TIP 710 
TIP 720 
TIP 730 
TIP 740 
TIP 750 
TIP 760 
TIP 770 
TIP 780 
TIP 790 
TIP 800 
TIP 810 
TIP 820 
TIP 830 
TIP 840 
TIP 850 
TIP 860 
TIP 870 
TIP 880 
TIP 890 
TIP 900 
TIP 910 
TIP 920 
TIP 930 
TIP 940 
TIP 950 
TIP 960 
TIP 970 
TIP 980 
TIP 990 
TIPIOOO 
TiPlOlO 
TIP1020 
TIP1030 
TIP1040 
TIP1050 
TIP1060 
TIP1070 
TIP1080 
TIP1090 
TIPIIOO 
TIPlllO 
TIP1120 
TIP1130 
TIP1140 
TIP1150 
TIPU60 
TIP1170 
TIP1180 
TIP1190 
TIP1200 
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no 

120 

130 

uo 

c 

c 

IbO 

160 


170 

180 

190 

200 

210 


220 


JK=1MM*1 

APhI=ATAM(PHI ( IM*1 ) ) 

SA=SIN<APHI) 

CA=C0S(APHI ) 

YLEG=Q(jm ♦St JK)*CA 

ir (NSSWI.EO.O) YLE6SV(IUU)=YLEG 

ZLEG=ZH(IM*1)^S(JK)»SA 

IF (NSSWl.EQ.O) ZLEGSV ( lUU) =ZLEG 

IF (XLnfT).EQ.XT(ITT) ) GO TO 150 

TLX1=PN < JK ) ♦$ ( JK ) *T AN (PSI ( JK ) ) *CA 

JK=JK*l 

TUX2=PN( JK) ♦$( JK) *TAN(PSI ( JK) ) *CA 

CRTTLG=TLX1-TLX2 

XTLEGINA)=TLX1/2.*TLX2/2. 

NAA=NA^1 

IF INSSWl .EO.NSSWSV ( 1 ) ) GO TO 110 
GO TO 130 
DO 120 11=2. L 

IF ((A0S<YLEGSV<IT)-YLEG)»UT,TOLC).AND.(A8S(ZLEGSV(IT) 

ILC)) CCIPS=CIRSUH ( IQ) 

IF (CCIRS.NE.O. ) GO TO 130 
CONTINUE 

DO 140 nV=NAA.NB 
NVT=NV-1 

XTlEG(NV)=XTLEG(NVT)-CRTTLG 

GO TO 40 


CONTINUE 

IF (CCIRS.NE.O. ) GO TO 160 

GO TO 270 

IJ=2*NSC«*1 

XLT=XTLE6tl ) ♦CLFTLG/2. 

XRT=XTLE&n J) *CRTTLG/2. 
XLL=XLT.CLFTlG/4. 

XRL=XRT^CRTTLG/4. 

IF (XLL.GE.XLt ITT) .AnD.XLT.LE.XT (ITT) ) 

IF (XLL.LE.XL(ITT) .ANO.XLT.GE.XT (ITT) ) 
IF (XLL .GT .XL ( I TT ) . ANO.XLT .GE . XL ( I TT ) ) 
IF (XLL.LE.XT ( ITT) ) GO TO 200 
IF (XLL.GT.XLUTT) .and. XLT.lt. XL (ITT)) 
C0N4=(XT(ITT)-XLL)/(XLT-XLL) 

GO TO 210 

C0N4= (XL tlTT)-XT(ITT))/ (XLL-XLT ) 

GO TO 210 

C0N4= (XL UTT)-XLT)/ (XLL-XLT.) 

GO TO 210 
C0N4=1 . 

GO TO 210 

CON4=0.0 

CONTINUE 

IF (XRL.GE.XL ( ITT ) . AND.XRT .LE.XT ( ITT) ) 
IF (XRL.LE.XL ( ITT) .AND.XRT .GE.XT(ITT) ) 
IF (XRL.6T.XL (ITT) . AND . XRT . GE . XL ( ITT) ) 
IF (XRL.LE.XT ( ITT) ) GO TO 250 
IF (XRL.GT.XL(ITT) . and. XRT.lt. XL ( ITT) ) 
C0N5=(XT (ITT)-XRL) /(XRT-XRL) 

GO TO 260 

C0N5= (XL (ITT)-XT(ITT))/t XRL-XRT ) 


GO 

TO 

170 

GO 

TO 

190 

GO 

TO 

^00 

GO 

TO 

IBO 


GO 

TO 

220 

GO 

TO 

240 

GO 

TO 

250 

GO 

TO 

230 


-ZLEG) .LT 


TIP1210 
TIP1220 
TIP1230 
TIP1240 
TIP1250 
TIP1260 
TIP1270 
TIP1280 
T1P1290 
TIP1300 
TIP1310 
TIP1320 
TIP1330 
TIP1340 
T IPl 350 
TIP1360 
TIP1370 
TIP1380 
TIP13R0 
TOTIP1400 
TIP1410 
TIP1420 
TIPU30 
TIP1440 
TIP1450 
TIP1460 
T IP1470 
T IP1480 
TIP1490 
TIP1500 
TIP1510 
TIP1520 
TIPI530 
TIP1540 
T IP1S50 
TIP1560 
TIP1570 
TIP15B0 
TIP1590 
TIP1600 
TlPlGlO 
TIP1620 
TIP1630 
TIP1640 
TIP1650 
TIP1660 
TIP1670 
TIP1680 
TIP16V0 
TIP1700 
TIP1710 
TIP1720 
TIP1730 
TIP1740 
TIP1750 
TIP1760 
TIP1770 
TIP1780 
TIP1790 
TIP1800 
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GO TO 260 




TIP1810 

^30 

C0N5=(XL ( ITT)-X«T) /URL-XRT) 




TIP1820 


GO TO 260 




TIP1830 

^40 

C0N5=l. 




TIP1840 


GO TO 260 




TIP1850 

^bO 

CON5=0,0 




TIP1860 

^60 

CONTINUE 




TIP1870 


T IPSU=TIHSU*CCIRS*0.25* (CON4»*(VO J< 1 )*CLFTlg- 

C0N5*WV0U( IJ)"CRTTLG) «TIP1880 


l2./SREr*dETA 




TIP1890 


P1TCH=PI TCH*CC IRS* 0.25* ( -C0N4*rfV3U ( 1 ) * 

clftlg 

«XTLEG( 1 ) 

♦ CONb^NVOU ( I JT IP 1900 


1 )*CRTTLG*XTLEG( IJ) ) *2./ (5REF*CREF) *0ETA** 

2 


TIP1910 

270 

C1RCUS=CC1RS 




TIP1920 


DO A60 nP05=1*NSCW 




TIP1930 


JK=IMM»NP05 




TIP1940 


JN=2*NSCW*NP0S 




TIP1950 


NPIS = NSC«K*NPOS 




TIP1960 


CIRCUS=CI«CUS»CIR UK,2) 




TIP1970 


IF (XL ( I IT ) .EQ.XT ( ITT) ) GO TO <*60 




TIP1980 


XLLEG=XTLEG(NP0S) 




TIP1990 


XRLEG=XTLE6( JN) 




TIP2000 


XLL = XTLEG(.MP0S) ♦CLFTLG/2. 




TIP2010 


XLT=XTLEG(MP0S)-CLFTlG/2. 




TIP2020 


XRL=XTLEG( JN) *CRTTLG/2. 




TIP2030 


XRT=XTLEG( JN)-CRTTLG/2. 




TIP2040 


IF (XLL.GE.XL( ITT) .AniO.XLT.LE.XT ( ITT) ) 

GO 

TO 

280 

TIP2050 


IF (XLL.LE.XL ( ITT) . ANO.XLT.GE.XT ( ITT) ) 

GO 

TO 

300 

TIP2060 


IF (XLL.GT.XLI ITT) .AND.XlT.GE.XL(ITT) ) 

GO 

TO 

310 

TIP2O70 


IF (XLL.LE.XT(ITT) ) GO TO 310 




TIP2080 


IF (XLL.GT.Xl (ITT) .and. XLT.lt. XL ( ITT) ) 

GO 

TO 

290 

TIP2090 


CONl=(XT ( ITT)-XLL) /(XLT-XLL) 




TIP2100 


XLLEG=XT (ITT) *C0N1*CLFTLG/2. 




TIP21 10 


GO TO 320 




TIP2120 

280 

CONl=(XL ( ITT) -XT (ITT) )/ (XLL-XLT) 




TIP2130 


XLLEG=(XL(ITT)*XT(ITT) )/2. 




TIP2140 


GO TO 320 




TIP2150 

290 

C0N1= (XL ( ITT) -XLT) / (XLL-XLT) 




TIP2160 


XLLEG=XL T* CON l*CLFT LG/2. 




TIP2170 


GO TO 320 




TIP2180 

300 

C0N1=1. 




TIP2190 


GO TO 320 




TIP2200 

310 

CON1=0.0 




TIP2210 

320 

CONTINUE 




TIP2220 


IF (NP0S.EQ.N5CW. AND. CONI .EO. 1 . ) GO TO 

360 


TIP2230 


IF (XRL.GE.XL ( ITT) . AnD . XR T .LE . XT ( ITT) ) 

GO 

TO 

330 

TIP2240 


IF (XRL.lE.XL (ITT) .AnD.XRT.GE.XT ( ITT) ) 

GO 

TO 

3b0 

TIP2250 


IF (XRL.GT.XL (ITT) . AnO . XRT . GE . XL ( ITT) ) 

GO 

TO 

370 

TIP2260 


IF (XRL.LE.XT ( ITT) ) 60 TO 370 




TIP2270 


IF (XRL.GT.XL (ITT) . ANO . XRT . LT , XL ( ITT) ) 

GO 

TO 

340 

TIP2280 


C0N2=(XT( ITT) -XRL)/ (XRT-XRL) 




TIP2290 


XRLEG=XT ( ITT) ♦C0N2*CRTTLG/2. 




TIP2300 


GO TO 380 




TIP2310 

330 

C0N2=(XL(ITT)-XT (ITT) )/ (XRL-XRT) 




TIP2320 


XRLeG= (XL ( ITT) «XT ( ITT) )/2. 




TIP2330 


GO TO 380 




TIP2340 

340 

C0N2= (XL (ITT)-XRT) / (XRL-XRT) 




TIP2350 


XRLE6=XPI*C0N2*CRTTLG/2. 




TIP2360 


GO TO 380 




TIP2370 

3b0 

C0N2=l . 




TIP23B0 


GO TO 380 




TIP2390 

360 

C0N1=.75 




TIP2400 
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CON2=.75 


TIP2410 


GO TO 380 


T IP2420 

370 

CON2=0.0 


TIP2430 

380 

IF (XRL.GT.XLL) GO TO 390 


TIP2440 


XSIGN=-1.0 


T IP2450 


XBLL=XLL 


TIP2460 


XBLT=X«L 


T IP2470 


GO TO 400 


TIP2480 

390 

X9LL=XRL 


TIP2490 


XBLT=XLL 


TIP2500 


X5IGN=1. 


TIP2510 

400 

BVDLG=ABG(XBLL-X8LT) 


TIP2S20 


IF (X8LT.GE.XL ( ITT) ) GO TO 440 


TIP2530 


IF (XBLL.LE.XT(ITT) ) GO TO 440 


TIP2540 


IF (XBLL.GE.XL(ITT) .AND.XBLT.LE.XT(ITT)) GO 

TO 

430 TIP2550 


IF (XBLL.LE.XL(ITT) .AND.XBLT.GE.XT(ITT) ) GO 

TO 

420 TIP2560 


IF (XBLL.GT.XL(ITT) .ANO.XBLT.GE.XT(ITT) ) GO 

TO 

410 TIP2570 


CON3=(XT ( ITT) -X3LL) /(XBLT-XBLL) 


TIP2580 


XTLEG(NPIS)=XT (ITT) *c0N3*BV0LG/2 . 


TIP2590 


C0n3=C0N3*XSIGN 


T IP2600 


GO TO 450 


TIP2610 

410 

CON3=(XL ( ITT)-XBLT) / (XBLL-XBLT ) 


TIP2620 


XTLEG(NPIS)=XBLT*C0M3*BV0LG/2. 


TIP2630 


C0N3=C0N3*XSIGN 


TIP2640 


GO TO 450 


TIP2650 

420 

CON3=l .*XSIGN 


TIP2660 


GO TO 450 


TIP2670 

430 

C0N3=(XL(ITT)-XT( ITT) ) / ( XBLL-XBL T ) 


TIP2660 


XTLEG(NPI5) = (XL ( ITT) .XT ( ITT) )/2. 


TIP2690 


C0N3=C0N3*XSIGN 


T IP2700 


GO TO 450 


TIP2710 

440 

COI^3 = 0.0 


TIP2720 

450 

TlPSU=TlPSU* (CIRCUS* {WV00(NP0S)*CLFTLG*C0 n1 

-C0N2*WV0U( JN)»CPTILG) -TIP2730 


1CIR( JK.2) * (*(VOU (NPIS ) *C0N3*BV0LG) ) *2./5REF* 

beta TIP2740 


PITCH=PI ICH* (CIRCUS* (-WVOU(NPOS) *CLFTLG*C0NI*XLLEG*WV0U( JN) »CUN2*CTIP27S0 


1HTTLG*XRlEG)-CIR ( JX*2)* (*(V0U(NPIS) *CON3*BVOLG* 

XTLEG(NPIS) ))*2./(SdTIP2760 


2EF*CREF) *BETa«»2 


TIP2770 

460 

CONTINUE ' 


TIP2760 


IM=IM*1 


TIP2790 


IMM=IMM*TBLSCrt ( IM) 


TIP2P00 


IF (NSSwi.EO.O) CIRSUM ( IM) =CIRCUS 


TIP2810 


IF (NSSWl.EO.O) IUJ=IM*2 


TIP2820 


IF ( IM.EQ.NSSWSV ( 1 ) ) GO TO 470 


TIP2830 


IF (XL ( ITT) .EO.XT ( ITT) ) GO TO 100 


TIP2840 


GO TO 480 


T1P2850 

470 

CTSW=TIPSU 


TIP2860 


cmw=pitch 


TIP2870 


IF (NSSW2.EQ.N5S»() GO TO 520 


TIP2880 


ITT=2 


TIP2890 


GO TO 10 


TIP2900 

480 

IF (IM.EQ.NSSW) go To 500 


TIP2910 


NCTL=1 


TIP2920 


DO 490 NV=1,NSCW 


T IP2930 


clftlg=c«ttlg 


TIP2940 


NY=NV*2*NSCW 


TIP2950 


XTLEG(NV)=XTLEG(NY) 


TIP2960 

490 

WVOO(NV)=WYOU(NY) 


TIP2970 


GO TO 80 


TIP2980 

500 

XKVSEH(2)=2.*ABS(TIPSU-CTSW) 


TIP2990 


IF (XKVSEW(2) .LT. 0.000001) GO TO 510 


TIP3000 
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CENTR(2) =(PITCH-CM»rf) «C^EF/A0S(TIPSU*CTSW) 

TIP3010 


60 TO 520 

TIP3020 

310 

CENTR(2)=0.0 

TIP3030 

5^0 

XKVSEW ( 1 ) =2,«ABS (CTSW) 

TIP3040 


IF (XKVSEW ( 1 ) .LT ,0.000001) GO TO 530 

TIP30S0 


CENTR( l)=CM^<>CkEF/A&S(CT5W) 

TIP3060 


GO TO 540 

TIP3070 

530 

CENTR(1)=0.0 

TIP3080 

540 

call rtRTANS (XKVSEiirftCENTR) 

TIPOO'^O 


RETURN 

TIP3100 


END 

TIP31 10- 
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I 


10 

^Q 


30 

40 


50 

60 


70 

BO 


SUBROUTINE WRTAN5 ( XK VSEii# ♦ CENTR ) 

COMMON /all/ BOTfM*BETA»PTESTf QTESTfTBLSC*^ (SO) »Q(200) t PN ( 200 ) » PV ( 2wRT 
100) f ALP(^OO) ♦S(200) *PSI (200) fPHi (50) *2H(50) WRT 

COMMON /OnETHPE/ T^IsT (2) fCPEFfSREFtCAVEf CLOESfSTRUEfARf AkTKUE fWTCWRT 
1DHT(2) *C0NFIG*NSSWSV (2) tM5V(2) t K30T t PL AN* IPL AN ♦ MACH» SSWWA ( 50 ) *XL(2wRT 
2) *XT (2) ♦CLW8*CMCL*CLA(2) t BLAIR (50) *CLAMAR(2) wRT 

COMMON /fHREFOR/ CCAV ( 2 * 50 ) » CLT *CLNT*N5Snf * ALPO WRT 

DIMENSION XKV(2)* XKP(2)f YCHLO(2)t YCHHI(2)* CENTR(2)» CEnT(2)* XwRT 
lKV5Ert(2) 

LCH=0 WRT 

LAMAR=NSbWSV ( 1 ) WRT 

SUMY=SUM=0*0 WRT 

C0nV=3. 1415926536/180. WRT 

CINV=1 ./(3.1415926S36*^AR) WRT 


DElTA=2.«C0NV 
C0NST=16.*B0T/SREF 
ALPHA = ALP0<»C0NV 
522 = ALPhA»<>2 
KBT = 1 

IF (KBOT.EQ.l) KBT=2 

IF (IPLAN.EQ.l) GO TO 10 

XKP(KBOT) =CLWB/ ALPHA 

XKP (KBT) = (CLDES-CLWB) /ALPHA 

GO TO 20 

XKP(1)=CLA(2) 

READ (5f 180) YCHLO < 1 ) *YCHHI ( 1 ) *YChL0 (2) ♦YCHhI (2) 
NCH=1 

NC1=N5SWSV(1) 

DO 40 J=lfNCl 

IF (Q(NCH).GT.YChLO(I) ) GO TO 30 
IF (Q (NCH) .LT. YCHHI ( 1 ) ) GO TO 30 


i^RT 

WRT 

WRT 

wRT 

WRT 

WRT 

WRT 

WRT 

WRT 

WRT 

WRT 

wRT 

WRT 

wRT 

wRT 

wRT 

WRT 


SUM: 

IT 


:SUM* (PN (NCH) ♦ ( °N (NCH) -PN (NCH^l ) ) /4* ) *BET A^BLAIR ( J) ( NCH ) ^CONSwRT 

wRT 


SUMY = 5UMY4BLAIR ( J) *^S(NCH)*CONST 
IF (J.EQ.NCl) GO TO 40 
NCH=NCH*TBLSCW ( J) 

CONTINUE 

XKV(l)=SUMY/S22 

IF (XKV ( 1 ) .LT. 0.000001 ) GO TO 50 

CENT(1)=^UM/SUMY 

GO TO 60 

CENT(1)=0,0 

CONTINUE 

IF (IPLAN*EQ.l) GO TO 100 
5UMY=SUM=0.0 
NCH=MSV(1)^1 
NC2=NSSWSV(2) 

DO 80 J=1*NC2 

IF (O(NCH) .GT.YCHL0(2) ) GO TO 70 
IF (CLINCH) ,LT. YCHHI (2) ) 60 TO 70 


WRT 

WRT 

WRT 

WRT 

wRT 

wRT 

WRT 

wRT 

wRT 

WRT 

WRT 

WRT 

WRT 

WRT 

wRT 

WRT 

WRT 


SUM = SUM^ (PN(NCH) ♦ ( :>N (NCH) -PN (NCH^ 1 ) ) /4. ) *BET A*BLA IR ( NCI ♦ J ) (NCH) »WRT 


I CONST WRT 

SUMY = SUMY^BLAIR (NCUJ)<^S(NCH) <»CONST wRT 

IF (J.EQ.NC2) GO TO 80 wRT 

NCH=NCH*TBLSCW (J ♦LAMAR) WRT 

CONTINUE WRT 

XKV(2) =SUMY/S22 WRT 

IF (XKV(2) .LT. 0.000001 ) GO TO 90 WRT 

CENT (2) =i>UM/SUMY WRT 

GO TO 100 WRT 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

4l0 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 
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90 

100 


110 

120 


130 


140 


IbO 


160 


CENT(2)=0.0 

*^RT 

610 

CONTINUE 

wRT 

620 

WRITE (6*190) 

wRT 

630 

DO 110 I<=1*IPLAN 

wRT 

640 

CEnTPM = ClAMAR ( IK ) *CREE 

wRT 

650 

WRITE (6*200) IK 

wRT 

660 

WRITE (6*210) XKP(IK) .CENTPM 

wRT 

670 

WRITE (6*220) XKV(IK) tCENT(IK) 

WRT 

680 

WRITE (6*230) XKVSEW(IK) ,CENTk(I<) 

WRT 

690 

CONTINUE 

wRT 

700 

DO 160 IK=1,IPLAN 

wRT 

710 

IF (LCH.EQ.l) GO TO 130 

WRT 

720 

WRITE (6*2A0) IK 

wRT 

730 

WRITE (6*250) 

WRT 

7^0 

ALPHA=0.0 

wRT 

750 

DO 160 J=l*26 

wRT 

760 

V=SIN(ALPHA) 

wRT 

770 

C=C0S(ALPHA) 

wRT 

780 

C2=C**2 

wRT 

790 

S2=V**2 

wRT 

800 

CLP=XKP(1K)*V*C2 

wRT 

810 

CLVL=CLP*XKO ( IK) *S2*C 

WRT 

820 

CLSL=CLP*XKVSEW ( IK ) «C*S2 

wRT 

830 

ClT0T=CLVL*XKVSEW(IK)*C*S2 

WRT 

8^0 

IF (LCH.EO.O) GO TO 140 

WRT 

850 

CMP=V*C* (XKP (2) *CLAMAR (2) ♦ (XKP( 1 ) -XKP (2) )*CLAMAR ( 1 ) ) 

WRT 

860 

CMPL = CMP*S2/C«EF“ (CENT (2)*XKV (2) ♦ (XKV( 1 )-XKV(2) ) *^CEnT ( 1 ) ) 

WRT 

870 

CMPS=CMP*S2/C«EF* (CENTR (2) *XKVSEw (2) ♦ ( XK VSEW ( 1 ) -XKVSEW ( 2 ) ) 

<^CENTR(iwRT 

880 


\ ) ) \ftRJ &90 

CHT0T=CMPL^CMPS-CMP 900 

GO TO IbO 910 

CMP=CLAMAR ( IK ) *XKP ( IK) *V*C 920 

CMPL=Cr4P^CENT ( IK) <^XKV ( I K ) <^S2/C9EF »^)PT 9 30 

CMPS=CHP^XKVSEW ( IK) *cENT9 ( IK) *S2/C«EF 940 

CMTOT = CMPL^XkVSEW( 1K)^*CENTP( IK) *52/CR£F W«T 9b0 

C0I=CLT0T*TAN (ALPHA) WRT 9b0 

CDII=(CLtOT*tt2)*ClNV WWT 970 

ALPHl=AL^HA/COMV WRT 980 

CNTT-CLTOT/C '"'RT ^'^0 

«^RITE (6t260) ALPHl*CNTT*CLP»CLVLfCLSLfCLTOTtC4P»CMPL»CMPS,CMT0T,CWRTl000 


170 


C 

c 

c 

ibo 

190 

200 


iDlfCDI I 

alpha=alrha^delta 

IF (IPLAN.EO.l) GO TO 170 

IPLAN=1 

LCH=1 

WRITE (6*280) 
XKP(1)=XKP(1)4XKP(2) 
XKVU)=XKV(1) ♦XKV(2) 

XK VSEW ( 1 ) =XKVSEW ( I ) ♦ XK V5EW ( 2 ) 
GO TO 120 
write (6*270) 

IPLAN=PLAN 

return 


format (4F10.5) 

FORMAT (lHl*////f31X,60HKP * KV AND RESPECTIVE CHOROWISE CENTROI 
1 FOR EACH PLANFORM) 

FORMAT (///*52X»12HPlANF0RM NO. *12) 


wRTiOlO 
WRT1020 
WRT1030 
WRT 1040 
wRTlObO 
WRTlObO 
WRT1070 
WRT1080 
WRT1090 
WRTllOO 
WRTUIO 
WRT1120 
WRTl 1 30 
WRTl 140 
wRTllbO 
WRTU60 
WRTU70 
OSWRTl 180 
WRT1190 
WRT1200 
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210 FORMAT (40X * 3HKP= ♦ F 1 0 . 5 1 1 OX . 1 1 hCENTROI D AT.F10.5) WRT1210 

220 FORMAT (37X.6HKV LE= ,F 1 0 ,5. 1 OX* 1 IHCENTROID AT.F10.5) WRT1220 

230 FORMAT {37X.6HKV SE= ♦ F 1 0 . 5 ♦ 1 OX * 1 1 HCENTRO lO ATtFlO.S) WRT1230 

2<*0 FORMAT ( 1 HI ♦////* A3X , <*0HPERF0RHANCE CHARACTERISTICS FOR PL ANFORM , I wRT 1 2^0 

12) WRT12S0 

2b0 FORMAT ( // 7X , 5HALPHA , 6X , 2HCN , BX . 3HCLP * <*X* 9 hCLP«CL VLE ♦ 1 * * 9HCLP*CL VSWRT 1 260 

1E*AX*2HCU*8X,3HCMP*AX.9HCMP*CMVLE*1X*9HCMP*CM\/SE»9X.2hCM.BX»2hCO.3wRTI270 
2X,13HCL**2/(PI»AR)/) WRT1280 

260 FORMAT (3x*12F10.4) WRT1290 

270 format (///.50X.21HTHIS CASE IS FINISHED) RRT1300 

280 FORMAT ( 1 HI ,////• 48X , 33HT0TAL PERFORMANCE CHARACTERISTICS) wRT1310 

END WRT1320- 


i 


i 

1 

1 
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TABLE I.- Ky se AND ITS CHORDWISE CENTROID OBTAINED 
BY FOUR METHODS FOR THREE PLANFORMS 


[m = 0; for method 2, Ng = 6 and Ng = 25; for method 3, Ng = 3 
and Ng = 25; for method 4, Ng = 11 and Ng = llj 


Method 

Cropped diamond 
planform 

Cropped arrow 
planform 

Cropped delta 
planform 

KvjSe 

Ax/ct 

*^v,se 

Ax/ct 

^v,se 

Ax/ct 

1 

1.200 

0.5367 

1.693 

0.5320 

1.397 

0.5373 

2 

1.232 

.5207 

1.726 

.5098 

1.456 

.5182 

3 

1.200 

.5451 

1.688 

.5350 

1.412 

.5441 

4 

1.300 


1.742 


1.60 



TABLE n.- EFFECT OF BOUND VORTEX SWEEP ANGLE 
ON Ky^se and Ax/ct 


[a = 3.5; \ = 1.0; M = 0.30; for method 2, Nq = 6 and Ns = 3^ 


A, deg 

Method 1 

Method 2 

^v,se 

Ax/ct 

^v,se 

Ax/ct 

0 

1.0968 

0.6004 

1.1037 

0.5956 

20 

1.3050 

.5777 

1.3869 

.5325 

40 

1.4149 

.5486 

1.3630 

.5270 

50 

1.3631 

.5418 

1.3076 

.5164 

60 

1.2431 

.5464 

1.1698 

.5066 

70 

.9446 

.5956 

.9243 

.4980 

75 

.8335 

.5729 

.7556 

.4916 
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TABLE in. - POTENTIAL AND VORTEX LIFT FACTORS OBTAINED 
FROM THE PRESENT METHOD AND METHOD 1 

[m = 0; for method 2, N^ = 6 and Ng = 2^ 


0.05 

.10 

.20 

.30 

.40 

1.00 


Type 


□ 


1.00 

G 

^873 

c 

^1.069 

<a 

^.738 

o 


Present method (method 2) 

Method 1 

Kp 

^v,fe 

*^v,se 

Kp 

^v,fe 

K^VjSe 

0.0798 

0.0399 

2.9816 

0.07844 

0.0393 

3.1799 

.1596 

.0798 

2.9477 

.1571 

.0785 

3.0188 

.3188 

.1597 

2.8533 

.3138 

.1571 

2.7913 

.4769 

.2395 

2.7497 

.4693 

.2356 

2.7208 

.6329 

.3194 

2.6467 

.6227 

.3141 

2.6341 

1.4862 

.7969 

2.1157 

1.4614 

.7816 

2.1255 

1.4475 

.7923 

1 

2.3581 

1.4335 

.7787 

2.3863 

1.3064 

1.5345 

1.4563 

1.2789 

1.5041 

1.3967 

1.5049 

1.8575 

1.7256 

1.4868 

1.8241 

1.6929 

1.1298 

1.3000 

1.2321 

1.1066 

1.2744 

1.2001 


^Same wings as presented in figure 2. 
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Vortex filaments 
Wing planform 
Elemental panel perimeter 


X 



End view of forces 


Figure 1, - General layout of axis systems, elemental panels 
and horseshoe vortices for a typical wing. 
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Figure 3,- Variatioi 
for simp! 
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O Experiment (ref. 4) 
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a, deg 

Figure 6.- Theoretical and experimental results for 0,4 -aspect -ratio rectangular flat 
wing with sharp leading and side edges at M = 0.20 with N^. = 6 and Ng = 25. 
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a. deg 

Figure 7.- Theoretical and experimental results for 1,0 -aspect -ratio rectangular flat 
wing with sharp leading and side edges at M = 0,20 with Nc = 6 and Ns = 25, 
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Figure 10.- Theoretical and experimental results for 1.069-aspect-ratio cropped arrow 
wing at M = 0.20 with Nc = 6 and N. = 25. 




Figure 11.- Theoretical and experimental results for 0.873 -aspect -ratio cropped delta 
wing at M = 0.20 with TTc = 6 and Hg = 25. 
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Figure 12.- Theoretical and experimental results on 44° swept wii^ in presence of canard in wing chord plane 
M = 0.30 with Nc = 6 and Ng = 12 for canard and Ng = 17 for wing. 
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Figure 14.- Theoretical and experimental results on 44° swept wing in presence of canard with 18.6° anhedral 
above wing chord plane. M = 0.30 with Nc = 6 and Ng = 12 for canard and Ng = 17 for wing. 




